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Atmospheric  Transmittance  from  0.25  to  26.5  ;/m 

Computer  Code  LOWTRAN  3 


i in  nusin  cri'  N 

Th<-  not’d  for  predicting  tin:  tr  ansmittanre  of  the  atmosphere  at  *..- !•  .«  broad 
spectral  interval  at  low  resolution  is  not  a new  one,  and  mum  methods  have  been 
proposed  to  do  this.  A major  problem  with  most  of  tile  techniques  if  that  the'  ar>- 
difficult  to  apple,  in  order  to  iMevi.it*  this  situ  ition  and  to  provide  a fuir!\  accu- 
rate, simple,  and  rapid  «.n  of  estimating  atmospheric  transmittance  in  th*’  0.2:i 
to  2H.fi  pm  rcginii,  an  ontpiric.il  graphical  prediction  scheme  was  <j  cased  using 
some  techniques  originally  suggested  hv  Altshuler.  The  prediction  scheme  is 
based  mainly  on  recent  laborato-y  tranimittance  measurements  -'omplt  rnented  hv 

using  available  the  theoretical  molecular  lint  constants  in  ’,ne-oy  - line  transmit- 

) 

tancc  calcu’ations,  and  is  presented  hv  Mel  latrhev  et  al, 

Hecause  of  the-  large  amount  of  interest  shown  in  this  work,  it  was  decided  !o 

computed 7*.-  the  prediction  scheme  and  to  digitize  the  spectral  curves,  tians- 

o 

mittancc  functions,  and  model  atmospheres  contained  m Met  latchey  et  al,"  which 
forms  the  basis  of  l.O'AT  HAN  i. 

( Received  for  punl  iration  li  Mai  l’MSi 

1.  Altshuler.  T.I..  ilDGP  Infrared  Transmission  arid  Background  Radiation  by 

Clear  Atmospheres.  vTl"  Report  til  Si  1 5 l'fffl,  A If  - 4 0 1 n 2 T. 

2.  Mcilatchey.  K.  A.  , I enn.  R.  W . . Selhv,  i.  1:.  A.  , Volz.  1.1...  and  daring  l.b 

< 1 P7  2 > Optical  Properties  of  the  Atmosphere  (Third  l.rtition',  Al  CHI  -72-0-J 


The  I ortr.in  computer  code  I .*  >V\  T H\\  3 is  designed  to  calculate  the  trans- 
mitt  nice  (averaged  over  .a  20  cm  ' interval)  for  a given  ntmospl  one  p..th  at 
Steps  of  5 cm  ' from  350  to  40,000  cm  ' <0.25  to  2h.fi  pint).  \ choice  of  six 
mode!  atmospheres  is  given  v.  itti  in  option  for  a seventh  model  which  can  lie  in- 
serted as  a set  of  radiosonde  data.  Aerosol  ittenuation  is  calculated  for  a given 
visual  range  based  on  an  interpolat ion  extrapolation  sch.eme  using  two  aerosol 
models  (See  Sect. on  3.1'. 

The  computer  code  1 i )U  l'UW  < supersedes  two  earlier  versions  of  the  pro- 
jr  ,1.,,  name  I’  I.DWTH  \A  1 * ami  I t >\\  T It  \ N 2.**  I ( 'U  i If  \ \ 2 is  a mod ifical ion 
of  the  1 .<  H\  T H iN  2 < omoutcr  code,  nid  provides  an  updating  of  the  original  data 
us  well  as  giving  more  flexibility  to  the  user.  The  differences  between  the  two 
programs  v.all  he  described  in  detol  -n  the  following  sections. 

f or  horizontal  path  transmittance  calculations  under  nonstandard  conditions, 
the  user  can  specif'.  Ins  own  meteorological  conditions.  2 he  amount  of  water 
Vapor  in  the  path  is  calcu’ated  m I a >\\  TltW  2 using  either  dew  point  temperature 
or  ambient  temperature  and  relative  humidity,  whichever  the  user  specifies  (set. 
Section  2.1). 

The  euo  sequence  numbering  system  used  m I .<  '\\  T !!  \ N 2 has  be  cn  preserved 
so  that  workers  who  are  1 1 ready  using  It'UTItW  2 can  update  tb.eir  card  decks 
with  o minimum  of  effort.  \ 11  changes  and  additions  to  the  latter  program  hoc 
been  indicated  by  a symbol,  (for  example',  . \.  H,  C,  etcetera)  against  an  ori- 
ginal sequence  number  isec  Appendix  M. 

We  will  first  briefly  describe  'lie  theory  and  input  data  used  in  the  program. 
General  instructions  for  using  l.<H\TI?W  2 are  given  in  S-ction  a.  A senes  of 
examples  illust  r iti'  g the  input  data  net  ess. rev  lor  making  n variety  of  typical 
atmospheric  transmittance  calculations  is  given  in  Section  6.  A listing  of  the 
computer  code  and  data  is  giver  in  Appendix  A.  supplemented  by  u How  chart 
(Appendix  lb  and  , detmition  of  symbols  (Appendix  I".  An  iterative  refraction 
Scheme  use!  for  one  particular  application  of  tne  program  (see  Section  f> . t> ' is 
described  m Appendix  ('.  Hxamples  of  atmospheric  transmittance  spectra  ob- 
tained from  I .<  >U  i It  A N 2 together  with  comparisons  with  laboratory  aid  field 
measurements  are  given  iri  Section  7. 

If  anv  discrepancies  are  encountered  m the  program,  we  would  npprec'ate 
notification  in  writing. 

3.  Marilt’v.  i ).  1’.  , Smith,  11.  i.  I’..  Treve,  f.  M.,  < a'penter,  I.W..  ilegges.  T.  ( 

I loan,  i . It.  (1^71)  nl'Tlli  li,  A I GUI  -7  1-0528  (Vol.  2 & 3* 

( 1 2 7 3 ) OITO  HI.  At  i HI  -Tlt-73- 02  17  and  04 0 1 
(1474'  OPTi  It  1IIH.  At  GUI  -Tit- 74-0319. 

4.  Selbv.  .1.11.  A.,  and  Met  Tatchcy.  H.A.  il‘'72'  At  niosplie  ri  c Trans  ini  ttanee 

from  0.25  to  28.  ' um-  Computer  c ode  1 ( H\  fTTN  ? . A TCR t .-1'  2 - 0T-f!5. 


2.  MODEL  ATMOSPHERES 


The  altitude,  pressure,  tempera*ure,  water  vapor  density,  and  ozone  density 

for  the  U.S.  Standard  Atmosphere  and  five  seasonal  model  atmospheres,  as  well 

3 

as  the  number  of  particles  per  cm  for  two  haze  models  — corresponding  to  sea 
level  visual  ranges  of  5 and  23  km  - are  provided  as  basic  input  data  for  LOW- 
TRAN  3.  The  model  atmospheres  correspond  to  the  1962  U.S.  Standard  Atmo- 

5 

sphere'  and  the  five  supplementary  models:  that  is.  Tropical  (15°N),  Midlatitude 

Summer  (a5“N,  July),  Midlatitude  Winter  (45°.  January).  Subarctic  Summer  (60°N, 

July),  and  Subarctic  Winter  (60°N,  January).  The  different  models  are  digitized 

in  1 km  steps  from  0 to  25  km,  5 km  steps  from  25  to  50  km.  then  at  70  km  and 

2 

100  km  directly  as  given  by  McClatchey  et  al. 

The  water  vapor  and  ozone  altitude  profiles  added  to  the  1962  l '.  S.  Standard 

g 

Atmosphere  by  McClatchey  et  al  were  obtained  from  Sissenwine  et  al  and 
Herring  et  al1  respectively,  and  correspond  to  mean  annual  values.  The  u/ater 
vapor  densities  for  the  1962  1J.S.  Standard  Atmosphere  correspond  to  relative 
humidities  of  approximately  50  percent  for  altitudes  up  to  10  km,  whereas  the 
relative  humidity  values  for  the  other  supplementary  models  tend  to  decrease 
with  altitude  from  approximately  8ft  percent  at  sea  level  to  approximately  30  per- 
cent at  10  km  altitude. 

In  addition  to  the  model  atmospheres  provided  in  this  report,  the  user  has  the 
option  of  inserting  his  own  model  atmosphere  (specifically  designed  for  direct  in- 
sertion of  radiosonde  data  (see  Section  6.9>),  or  of  building  another  model  by 
combining  various  parts  of  the  six  standard  models  (see  Section  6.10). 

One  major  difference  between  LOWTKAN  3 and  LOWTRAN  2 is  that  the 
reader  no  longer  has  to  look  up  the  saturation  vapor  density  of  water  (from  Table  1 

4 

in  Selby  and  McClatchey  ) when  using  meteorological  data  ns  input  f»  tiie  program. 
In  ...OWTRAN  3 the  actual  water  vapor  density  is  calculated  for  a given  ambient 
temperature  and  relative  humidity  or  dew  point  temperature  using  the  following 
empirical  expression  for  the  saturation  vapor  density: 


F(t)  = A exp  (18.9766  - 14.9595A  - 2.4388A 


2 1 -3 

! gm  in 


5.  Valley,  S.  L.  , Ed.  (1965)  Handbook  of  Geophysics  and  Space  Environments, 

AFCRL,  

6.  Sissenwine,  N, , Grantham,  D.  D, , Salmela,  H,  A,  (1968;  Humidity  Up  to  the 

Menopause,  AFCRL-68-0550. 


7.  Herring,  W.  S. , and  Borden,  T.  R.  (1964)  Ozone  Observations  Over  North 
America,  AFCKL-64-30,  Vol.  2.  " 
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v lere 


A -=  27  3. 1 5/(27  3.15  + t) 

and 

t is  given  in  °C  . 

The  above  expression  was  found  to  give  a good  fit  to  published  values  of  sat- 
uration water  vapor  density  measured  over  water  to  better  than  1 percent  for 
temperatures  between  -50°C  and  5t)“C. 

If  t is  the  dew  point  temperature  then  F(t)  gives  the  actual  water  vapor  den- 
sity directly  at  the  corresponding  ambient  temperature.  If  * is  the  ambient  tem- 
perature then  the  water  vapor  density  is  given  by  F(t)  X III]/ 100  where  Rli  is  the 
percent  relative  humidity. 

Thus  the  user  has  a choice  of  meteorological  parameters  necessary  to 
specify  the  amount  cf  water  vapor,  that  is,  ambient  temperature  and  relative 
humidity,  or  dew  point  temperature,  or  water  vapor  density.  The  procedure  for 
inserting  radiosonde  data  into  the  program  and  the  necessary  formats  are  des- 
cribed in  Sections  5.3  and  6.8. 

a.  ATMOSPHERIC  CONSTITUENTS 


3,1  Atmojiplieric  (,j.n 

It  is  assumed  in  this  report  that  mixing  ratios  of  the  gases,  COg  . N?0 , , 

CO,  N., , and  00  remain  constant  at  all  altitudes  at  the  following  values:  330, 

i - 5 >; 

0.28,  1.6,  0.0,5,  7.505  V 10',  and  2.005  > io'  parts  per  million  respectively. 
These  gases  as  a whole,  w ith  the  exception  of  nitrogen,  w ill  be  inferred  to  as  the 
unifoimlv  mixed  gases. 

Absorption  coefficients  for  water  vapor,  ozone,  and  the  cor r.V  .,ied  effects  of 
the  uniformly  mixed  gases  were  digitized  from  the  spectra'  .urves  (Figures  16- 
26)2  Mrt']n(cf,ev  et  a]  and  3re  included  as  data  for  l.OWTKAN  3.  The  trans- 
mittance spectra  from  which  the  coefficients  were  derived  were  first  degraded 

in  resolution  to  20  cm  ' and  the  data  points  were  digitized  at  steps  of  5 cm  . 

2 

l or  the  ultraviolet  and  visible  ozone  bards  (see  McClatchey  et  al.  Figure  26), 
the  absorption  coefficients  were  digitized  at  200  cm  1 and  500  cm  ' intervals 
respectively. 


8.  List,  H.  I.  (1568)  Smithsonian  Meteorological  'Tables  (6th  revised  edition', 
Smithsonian  Institute  Press,  Washington. 
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* 

There  has  been  one  modification  to  the  spectra!  data  for  water  vapor  in  the 
2.7  pm  region  (since  LOWTRAN  2 was  published)  based  on  a more  recent  review 
of  available  experimental  measurements.  The  above  modification  leads  to  slightly 
higher  atmospheric  absorption  by  water  vapor  in  the  2.  2 to  3,4  pm  spectral  region 
than  given  by  LOWTRAN  2. 

3.2  Continuum  Abn^rptMii 

Absorption  coefficients  for  the  water  vapor  continuum  near  ’0  pm  and  4 pm 

in  LOWTRA  N 3 are  based  on  measurements  of  Burch  et  al,  McCoy  and  Rensch, 
9-13 

and  Bignell.  The  effect  of  absorption  by  the  water  vapor  continuum  between 

14  p m and  28.5  pm  has  not  been  included  at  this  time  because  there  ie  insufficient 
data  in  this  region. 

The  continuum  due  to  collision  induced  absorption  by  nitrogen  in  the  4 pm 

14 

region  is  included  in  LOWTRAN  3 based  on  the  measurements  of  Reddy  and  Cho 

15  2 

and  Shapiro  and  Gush  ' (see  also  McClatchey  et  al  ). 

In  all  cases  the  transmittance  due  to  continuum  absorption  is  assumed  to 

follow  a simple  exponential  Law. 

3.3  ftioieiuLr  .5eatlenng 

The  absorption  coefficient  due  to  molecular  scattering,  C6.  is  introduced 
into  LOWTRAN  3 via  the  following  expression: 

C6  =■•  9,87  X 10"20 

where  v is  in  wavenumbers  (cm  *). 

The  above  expression  was  obtained  as  a best  fit  to  molecular  scattering  no- 
efficients  published  by  Penndorf1”  and  is  shown  in  Figure  1 together  with  the 
aerosol  extinction  coefficient. 

3 4 Vroial  ModeU 

Two  aerosol  models  are  incorporated  into  LOWTRAN  3 and  correspond  to 
visual  ranges  of  approximately  5 km  and  23  km  at  sea  level.  However,  an  aerosol 
attenuation  for  any  visual  range  is  calculated  by  LjOWTRAN  3 using  an  interpola- 
tion/extrapolation procedure  (described  in  Section  3.  5)  which  utilizes  these  two 
models.  As  these  aerosol  models  are  based  on  measurements  of  continental 
aerosol  under  moderate  visibility  conditions,  they  may  rot  be  valid  for  very  low 

No-  included  in  LOWTRAN  2. 

NOTE:  For  references  9-16,  see  list  of  References  on  page  63. 


1 


WAV£NUMBEA  Ion  1] 


SO. 000  25.000  10.000  5.OJ0  2,600  1.000  500  300 


001 


100 


Figure  1.  Attenuation  Coefficients  for  Aerosol  Transmittance 
(Absorption  and  Total  Extinction) 


visibility  cond'tions  less  than  2 km,  Reported  low  visibility  conditions  less  than 
2 km  are  probably  representative  of  fog  conditions.  In  this  case,  the  LOU.  TRAN  3 
results  will  tend  to  underestimate  the  attenuation  (overestimate  the  transmittance) 
in  the  infrared  and  overestimate  the  attenuation  in  the  ultraviolet.  A more  realis- 
tic result  can  be  obtained  by  assuming  the  attenuation  to  be  wavelength  independent 
(see  McClatchey.  et  al,  1972,  p.  7 9)^  and  assuming  that  the  aerosol  attenuation 
provided  by  IjOWTRAN  3 at  5500  J is  also  valid  throughout  the  infrared  and  m ar 
ultraviolet.  The  application  of  this  "fog  model"  should  only  be  applied  to  the  low- 
est few  hundred  meters  above  the  surface.  In  IXJWTRAN  3,  a message  will  be 
printed  out.  to  warn  the  user  in  the  event  that  results  are  required  for  sea  level 
visual  ranges  less  than  2 km.  The  two  aerosol  models  are  based  on  the  following 
assumptions. 

17  18 

(1)  A particle  size  distribution  similar  to  Z)eirmeniji;.n's  Raze  Model  C.  * 
but  where  the  large  particle  radius  cutoff  has  been  extend-d  to  100  ;i m (compared 

17.  Deirmendjian,  D.  (1964)  Appl.  Opt.  JL187. 

18.  Deirmendjian.  D.  (1969)  Electromagnetic  Scntierin^  c.~*  Spherical  8\Mv dis- 

persions, American  Elsevier  Co.',  H’.T . " 


U 


to  5 g m in  Peirmendjian1  ~ and  10  (im  in  McClatche.v  et  al^  and  Selby  and 
McClatchey^ ). 

(21  The  particle  size  distribution  is  assumed  to  remain  constant  with  altitude. 

<31  The  variation  of  aerosol  number  density  with  altitude  is  assumed  to  be 

2 

the  same  as  previously  given  by  Mci'latchey  et  al  for  the  23  km  visual  range 
model.  The  latter  aerosol  number  densities  were  adjusted  to  give  extinction  co- 
efficients at  a wavelength  of  0.55 /im  that  corresponded  to  those  obtained  by 
19  20 

Elterman  ' ’ at  each  altitude. 

(4)  The  variation  of  aerosol  refractive  index  with  wavelength  has  been  ob- 

2 1 22 

tained  from  measurements  by  Volz  (see  also  McClatche.v  and  Selby  1,  who 
has  found  that  aerosols  are  composed  of  water-soluble  substances  as  well  as 
dust  like  material, 

Aerosol  extinction  (07!  and  aerosol  absorption  (07 A!  values  -.ere  calculated 
based  on  single  scattering  Mie  theory  using  the  above  aerosol  size  distribution 
and  refractive  index  values  (assuming  the  aerosols  to  be  composed  of  70  percent 
water-soluble  substance  and  30  percent  dust-like  substance,  which  appears  to  be 
representative  of  continental  aerosol),  Figure  1 shows  the  variation  of  the  calcu- 
lated aerosol  extinction  and  absorption  coefficients  with  wavelength.  In  I.OW- 
TRAN  3,  (C7!  and  (C7A)  were  digitized  directly  from  Figure  1 at  discrete  wave- 
lengths (see  Appendix  A and  Table  A 21. 

The  above  aerosol  model  replaces  the  empirical  function  previously  used  in 

I.OWTRAN  2,  f igure  1 can  be  used  to  calculate  aerosol  extinction  and  absorption 

2 

in  the  same  way  as  described  in  McCIatchey  et  al  and  replaces  Figure  22  in  the 
latter  reference. 


3.5  Vniiwl  Intnpolitioii/bxlrcpolition  Scheme 


The  total  extinction  coefficient  0„.  at  0.55  um  is  inverselv  oroDortiona!  to 

1 23 

visual  range,  VIS,  and  can  be  w,  itten  as  follows  (Middleton  ); 


m 


3.91 
VIS  * 


19.  Elterman,  L.  (1968)  UV,  Visible  a.-.d  111  Attenuation  for  Altitudes  up  to 
50  km,  AFCHL-ee^UTST: 


20, 

21. 


22. 


23. 


Elterman,  I,.  (1970)  Vertical  Attenuation  Model  with  Eight  Surface  Meteoro- 
logical Ranges  2 to  13  km,  ArCRL-7t)-t)Z0b. 

Vole,  F.  E.  (1972)  Appl.Opt.  lj:7r5. 

McCIatchey,  H.  A,,  and  Selby,  .I.E.A.  (1974)  Atmospheric  Attenuation  of 
Laser  Radi ation  from  0.76  to  31.25  pm,  A id?  T.-TR - 7 4 

Middleton,  W.  E.  K.  (1952)  Vision  Through  the  Atmosphere,  Univ.  of  Toronto 
Press. 
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assuming  a 2 percent  contrast  threshold  where  the  suffixes  a and  m refer  to  the 
aerosol  and  molecular  components  respectively.  The  aerosol  extinction  coefficient 
can  thus  be  written  as 


a 


a 


3.91 

VIS 


V 


m 


Since  the  aerosol  extinction  coefficient  <7  is  directly  proportional  to  the 
aerosol  number  density  N(z),  we  can  write 

N(zl  = — - biz)  . 

where  a(z)  and  b(z>  are  constants  for  a given  altitude  z.  It  will  be  noted  that 
b(z*  is  proportional  to  the  molecular  scattering  coefficient  at  0.55  pm  at  altitude 
z,  where  molecular  absorption  has  been  assumed  negligible  at  A = 0.55pm. 

The  above  equation  forms  the  oasis  for  the  interpolation/cxtrnpolation  pro- 
cedure used  in  I.OWTRAN  3 to  determine  aerosol  attenuation  at  any  given  \ isual 
range. 

The  coefficients  a and  b are  determined  from  the  nhove  equation  at  each 
altitude  using  vhc  two  aer r .ol  moqels.  that  is. 

a(z'  » |Ng<r)  - *.23(z>j  / | 1 / 5 - 1/23  | 

biz)  = j N23(z)/5  - Ns<z>/23j  J jl/5  - 1/23  j 

where  Nj.  and  N,,^  refer  to  the  number  densities  for  the  5-ktn  and  23-km  sea 
level  visual  ranges.  Note  that  the  above  procedure  is  used  only  in  the  lower  5 kn 
of  the  atmosphere  since  the  two  aerosol  models  arc  identical  above  5 km  altitude. 


4.  THKOR1 

4.1  B«uc  Awm^Kns 

The  computer  program  LOWTfiAN  3 follows  almost  exactly  the  procedures 

2 

outlined  by  McClatchey  et  al.  The  main  assumptions  made  are  that  the  ntmo- 
spher  can  be  represented  by  a 33-layer  model,  and  that  the  average  transmit- 
tance T over  a 20  cm  * interval  (due  to  molecular  absorption)  can  be  represented 
by  a single  parameter  model  of  the  foi  m 

T = f (C„W*)  (1> 
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where  0^  is  a wavelength  (or  wavenumber)  dependent  absorption  coefficient  and 
a)  is  an  "equivalent  absorber  amount"  for  the  atmospheric  path,  which  is  defined 
in  terms  of  the  pressure  P(z).  temperature  T(z),  concentration  of  absorber  AL, 
and  an  empirical  constant  n as  follows: 


<u 


AI. 


I ft  7.1 

i Po 


(2) 


If  Eq.(2)  is  substituted  in  Eq.  (1)  and  n is  set  equal  to  zero  and  unity,  respec- 
tively, Eq.(l)  reverts  to  the  well  known  weak  line  and  Strong  line  approximations 
common  to  most  band  models. 

The  form  of  the  function  f and  parameter  n was  determined  empirically 
using  both  laboratory  transmittance  data  and  available  molecular  line  constants. 

In  both  cases,  the  transmittance  was  degraded  in  resolution  to  20  cm  1 through- 
out the  entire  spectral  range  covered  here.  It  was  found  that  the  functions  f for 
11,0  and  the  combined  contributions  c.f  the  uniformly  mixed  gases  were  essen- 
tially identical,  although  the  parameter  n differed  in  the  two  cases.  Mean  values 
of  n were  determined  to  he  0.9  for  11,0.  0.75  for  the  uniformly  mixed  gases. 

0.4  for  ozons. 


4.2  F»rt!i  (jirwturr  and  Rrlrjrlion 

In  general,  earth  curvature  has  a greater  influence  on  the  path  length  (and 
hence  on  the  transmittance)  than  atmospheric  refraction.  For  long  slant  paths 
with  zenith  angles  close  to  90"  in  the  lower  layers  of  the  atmosphere,  however, 
refractive  effects  can  cause  a significant  increase  in  the  path  length  (up  to  30  per- 
cent for  . 90"  path  to  space  from  ground  level).  Figure  2 shows  the  effect  of 
atmospheric  refraction  on  defining  the  minimum  height  of  a path  trajectory  f’om 
space.  The  minimum  height  referred  to  here  is  also  known  as  the  tangent  height. 

In  Figure  2,  the  difference  between  the  geometrical  (no  refraction)  and  the  actual 
minimum  height  is  plotted  against  the  actual  minimum  height  for  three  different 
model  atmospheres.  The  sketch  in  the  upper  right-hand  corne~  of  Figure  2 indi- 
cates that  there  is  also  a discrepancy  in  the  earth  center  angle  P subtended  by  the 
trajectory,  when  refraction  is  significant.  The  difference  P - P'  shown  in  f ig- 
ure 2 is  equal  to  the  total  angular  deviation  1 1>  of  the  trajectory  due  to  refraction. 

For  many  applications  it  is  necessary  to  account  not  only  for  the  effect  of  re- 
fraction and  earth  curvature  on  the  transmittance  over  a given  path  tra  jectory,  but 
also  on  the  purely  geometrical  aspects  of  the  trajectory  itself,  l or  example,  the 
total  deviation  angle  of  arrival  6,  or  angle  P subtended  by  the  path  trajectory 
mav  be  required  as  illustrated  in  figure  3.  I.OWTHAN  3 calculates  the  quantities 
4.  P and  slant  range  on  the  basis  of  a layered  atmosphere  in  the  following  para- 
graphs. 
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I- igurc  3.  (ieneral  Schematic  of  a Kefracted  Path  From  Altitudes 
HI  to  112  Showing  the  Angles  Defining  the  Trajectory.  Initial 
zenith  angle  0o  at  III,  angle  of  arrival  6 at  112,  total  angular 
deviation  i,  and  angle  subtended  by  path  at  the  earth's  center.  P 

The  earth's  atmosphere  is  assumed  to  he  divided  into  a series  oi  concentric 
spherical  layers  for  each  of  which  a mean  refractive  index  is  defined.  However, 
the  non-sphericity  of  the  earth  is  taken  into  account  to  some  extent  bv  using  a dif- 
fereci  earth  radius  for  each  latitude  (associated  with  a given  model  atmosphere!. 

Consider  the  trajectory  of  a ray  passing  from  heights  111  to  112  at  an  initial 
zenith  angle  flQ.  I.et  z.  and  z^  define  the  boundary  heights  of  a given  layer, 
and  let  0.  and  tl.  he  the  local  zenith  angles  at  the  respective  boundaries  (see 
F igure  31.  Then  at  a height  of  z.  . , the  angle  of  refraction  is  0.  . The  angle 
of  incidence  oc  at  height  z.+j  can  he  defined  as 

Sin  or  = (Kq  + z()  sin  (?.  /ji:oa  zi+1).  (3> 


17 


Applying  Snell's  law  at  boundary  z.+J  . we  have 

n sina.  = n , sin  t).  , 
i 1 i+I  i+I 

where  ^ and  n)+]  are  the  mean  refractive  indices  of  the  layers  above  z(  and 
z.+l  respectively. 

Substituting  for  sin  3.  in  Lq.  (41.  we  have 

n (It  + z.  i sin  0-  = n.  . jit  + z ,1  sin  . 

1 ' o ii  1 1+1  1 o i+l)  i+l 

It  follows  from  symmetry  that 


11  ( It  + z.  I sin  0 = n.  , (It  + z , | sin  0.  , 
1 ' o 1'  1 1-1  ' o 1-1*  1-1 


= n (It  + Ml ) sin  0 
0 0 1 o 


= const. 


therefore,  the  angle  of  refraction  at  any  level  z can  be  written  in  terms  of  the 
initial  input  conditions  and  the  refractive  index  nQ  of  the  layer  above  HI  as 

Sin  0 = n0  (i?Q  + 111)  sin  ()  /n  ( lto  + zj  . (71 

The  angle  0.  subtended  at  the  center  of  the  earth  by  the  intersection  of  the  ray 
with  the  layer  zf  to  7.  ^ is  given  by 

= 0i  - a|  . (8) 

Thus  the  total  earth  center  angle  subtended  by  the  ray  when  traversing  the  atmo- 
sphere from  111  to  H2  is 
m-  1 

* = ? (fli  ‘“ii  <9’ 

m-  1 . . 

= T |sin'1{A/ni(Ho+  zj}  - sin’1  jA/ni(Ho+  zi+J,U  (10) 


where  in  is  the  number  of  levels  between  111  and  112,  and  A = n (H  HI)  sin  0 , 

00  o 


The  angle  of  arrival  6 of  the  ray  at  i I ^ <s  given  by 

6 = 180°  - sin'1  | A / n , (It  + H2)l  . 

| m-1  ’ o if 


18 


The  total  angular  deviation  of  the  trajectory  ^ is»  given  by 


* = P - * - 0o  * 180  . 


The  effective  path  length  between  levels  7,  and  if  given  by 

PS.  = (It  -17  sin  £ /sin  9 far  0"'  •-  0 < 180' 

l ' o i + l 1 1/  i 


(12' 


(13) 


for  0 - 03  and  180°,  I)S  = a - i . If  we  assume  that  the  equivalent  absorber 

i i+i  i 

amount  per  unit  path  length  <*;<see  section  4.1)  for  a given  gas  varies  exponential!', 
with  altitude,  we  can  write 


\ « d«  = ii,  [wl^i  - «{*,,,)] 


<14' 


where  II,  = (71+]  ' )/loge  (z, ) /w{71+])j  • The  'mount  of  absorber  V\  . along  a 
path  of  length  t)S,  between  altitudes  and  z,^  is  therefore  given  by: 

!)S, 

V\  , s \ OJ  ds 


I)S  1+, 

;—hr  \ a.dz 

! + 1 1 

Z1 

l)S,[^P,|-  “MVi1] 

,0«e[UJlZi)/W(Z.+l  )]' 


(15' 


The  total  equivalent  absorber  amount  V\  for  a given  atmospheric  path  is  given  by 

m - 1 

the  sum  of  the  U values  for  all  layers;  that  is,  tt  = L \5  where  m is  the 

i ! 

number  of  levels  traversed  by  the  path. 
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4.3  Rffitcliir  Index  of  \ir 

The  following  simplified  version  of  r.dlen’s""*  expression  for  the  refractive 
index  of  air  is  used  in  I.iHVTH.XN  3: 


Pjj  ( ) 

/ „ + 5 0.4sn  \ I*  2 

(n,  - ■) 10  ■ ("•">  * -nr)  r ' 


(<,<:  ■ 5-i). 
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24.  Edlen,  H.  < 1 f>66>  Mctrologia  2:12. 


in 


where  P[j9(i  and  I’  refer  respectively  to  the  partial  pressuie  of  water  vapor  and 
atmospheric  pressure  in  millibars.  T is  atmospheric  temperature  in  degrees 
Kelvin,  and  A is  the  w ivelength  in  micrometers  (pm). 

The  above  expression  lias  been  used  over  the  entire  wavelength  range  0.2  to 
28.5  prr.  in  LOW  THAN  3.  Although  fallen's  expression  for  the  refractive  aide., 
of  air  is  widely  used  in  both  the  visible  and  infrared  spectral  regions,  it  is  ques- 
tionable how  far  it  should  be  used  into  the  ultraviolet  and  into  the  far  infrared 
since  the  formula  is  based  primarily  on  measurements  made  in  the  visible  part 
of  the  spectrum  from  L.4.3  to  0.8  pm. 

4.4  (^oniftncil  Path  (iHilipii«liinv 

W hen  using  i.OUTllAN  3,  the  type  of  atmospheric  path  for  which  a calculation 
is  to  be  made  must  be  specified  according  to  one  of  the  three  broad  categories 
listed  below. 

TV  PM  1.  Horizontal  path:  that  is,  a constant  pressure  path  where  tin 
effects  of  earth  curvature  and  refraction  arc  negligible. 

TYHf  2.  Slant  paths  between  two  alt  tudes  from  HI  to  M2. 

TYIJf.  3.  Slant  paths  to  space  from  initial  altitude  iii. 

The  variations  within  the  latter  two  categories  for  both  upward  and  downward 
path  trajectories  can  be  seen  from  figure  4. 

It  will  be  noted  that  two  trajectories  are  possible  for  a given  set  of  input 
parameters.  Ml,  M2,  and  fl  for  a downward  looking  path  (TA  Mil  2',  provided  that 
112  lies  between  ill  and  the  minimum  height.  IIMI\. 

In  most  instances,  the  reader  will  not  he  aware  that  two  paths  arc  possible 
for  a given  set  of  input  conditions.  I or  such  i case,  is'U  Tli\\  3 will  execute 
the  shorter  path  condition  | figure  4id>]  and  print  out  ..  message  to  the  effect  that 
the  case  shown  in  figure  4(e*  does  exist.  Should  the  reader  decide  to  run  the 
latter  cose,  he  need  only  set  the  parameter  I.I1N  equal  to  unity  and  resubmit  the 
cose.  This  will  be  seen  more  clearly  from  the  following  section  (also  Section  6.4'. 

5.  INSTRI  CTIONS  FOR  l SIM;  ..OUTRAN  .3 

The  input  data  for  I.OV  THAN  .3  arc  given  in  Appendix  A.  In  general,  il  is 
only  necessary  to  change  the  last  four  raids  (referred  to  here  as  Nos.  1-4'  + in 
order  to  run  the  program  for  a given  problem.  The  formats  for  the  last  four 
cards  and  their  application  will  next  be  discussed. 


t These  four  cards  were  referred  to  as  Nos, 4, 0-4,3  in  the  l.OUTKAN  2 report. 


I ipure  4.  (iconitlnral  1 ’ it li  t'onfipiirat .01.  for  1.1'  Horizontal  Paths  ITvpe  l1, 
O'  Slant  Paths  llitween  Two  Altitudes  HI  >nd  illi  iTvpe  J1.  ami  let  Slant  Paths 
to  Sp  ire  ITvpe  :i>.  I or  dovuiw  iro  looking  piths  where  II  MIN  •'  111’  ■-  111  . two 
f 1 aiectorics  are  possible  is  indicated  in  t :>  and  U-' 
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5.1  bifHit  Da^  aiui  KmiijI.' 

The  data  necessary  to  specify  a 

given  prouleni 

are  givet.  on  tne  last  four 

1 

? 

cards  as 
(AHI)  1 

follows 
\H  (DM.. 

iiiAzi:.  iTYpp, 

I.FX,  UP,  IM, 

Ml,  M2.  M3,  Ml,  HO 

4 

\ 

i 

( AHI)  2 

111.  M2, 

Asci  i:.  HAW.  12 

i;i:ta.  vis 

{ FORMA!  (1013,  1 ! 0.3') 
{ 1'OHMAT  ((>  P I 0.3)} 

■i 

( AHI)  3 

\ 1,  \ 2, 

DV 

{ FORMAT  (3110.3)} 

i 

(’A  HI)  4 

IS  1 

( FORMAT  ( 13)). 

J 
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Definitions  of  the  above  quantities  will  be  discussed  in  Section  5.2, 

If  the  quantity  MO!)  FI.  given  on  C.AIilJ  1 is  set  equal  to  0 or  7 i which  is  the 
case  if  meteorological  data  are  used  as  input  to  the  program),  then  the  above 
card  sequence  (arid  format  for  CARD  2)  is  changed.  These  oases  mil  be  des- 
cribed in  Section  5.5  and  examples  for  some  typical  problems  will  lie  given  in 
Sections  (i.7  and  (i.H. 

5.2  IUmic  ln*in»t-|M>n« 


The  various  quantities  to  be  specified  or.  each  of  the  four  control  eai  do  (sum- 
marized in  Section  5.1)  will  be  discussed  .n  this  section. 

5.2.1  CARP  1 MODI  I , 'j A /,K.  ITYPF.  I.FX.  .ID,  IM, 

Ml.  M2,  M3.  MI.,  HO 

The  parameter  Mill'll!,  selects  cue  of  the  six  geographic  .1  odei  atmo- 
spheres or  specifies  that  meteorological  data  are  to  be  used  in  pi  ice  of  the  stan- 
dard models.  IIIAZI1  specifies  •..nether  aerosol  attenuation  is  to  be  included  in 
the  calculation  or  not.  l or  nay  problem  the  atmospheric  path  must  be  specified 
as  one  of  three  types  according  to  ITYDI1  and  I.FX,  The  rest  of  the  quantities 
given  on  DA  HI)  1*  (which  can  he  b’ft  blank  if  not  required'  provide  the  user  with 
options  to  suppress  printing  tUI’),  to  intermix  the  six  standard  model  atmospheres 
(Ml,  M2,  MM  a no  to  input  a new  model  atmosphere  (IM,  Ml.'.  The  options  for 
the  above  parameters  and  their  uses  are  stated  and  described  in  detail  below: 

MOI>Fi.  - 0 if  meteorological  data  are  specified  (for  horizontal  paths  onlyH 

- 1 selects  TIMDIDAI.  Mi)Di:i.  ATMOSPIIFRF 

= 2 selects  M II H .ATIII  I >b  si  MAI  I : H 

= 3 selects  M II'l .ATITI  DI 1 U l\T D It 

= 4 selects  St  PARC TK'Sl  MMCI! 

= 5 selects  41  HARCTIC  MIXTFK 

- b selects  1062  I S STANDARD 

= 7 if  a new  model  atmosphere  (or  radiosonde  data)  is  to  be  inserted.  1 

t The  parameters  UP,  IM,  etc.  given  on  CARD  I were  not  included  in  I.OWTRAN  2. 
I In  thise  cases  the  format  for  t’Alih  2 changes  (see  nonstandard  conditions) 

Section  5,3  and  examples  7 and  }.  (Section  6.7  and  6.8). 


■ 

J 
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IHAZE  = 0 means  no  aerosol  attenuation  included  i l the  calculations. 

IHAZE  = 1 or  2 if  aerosol  attenuation  is  required  (see  also  CARD  2). 

If  I HAZE  is  set  equal  to  1 or  2 and  visual  range  (\  IS)  is  not  specified  on  CARD  2, 
then  the  program  will  automatically  select  visual  ranges  of  23  km  or  5 km  re- 
spectively. 

ITYPE  = 1 for  a horizontal  (constant  pressure)  path. 

= 2 for  a vertical  or  slant  path  between  two  altitudes. 

= 3 for  a vertical  or  slant  path  to  space. 

The  TYPE  1 path  should  not  be  confused  with  a long  90°  pat;,  where  the  local 
height  of  the  end  of  the  trajectory  is  at  a significantly  different  height.  In  such  a 
case,  specify  the  pclh  according  to  ITYPE  = 2 (see  Section  5.2.2). 

LEN  = 0 for  normal  operation  of  program. 

LEN  = 1 selects  the  downward  TYPE  2 path  shown  in  Figure  4(e). 

The  parameter  I.EN  can  be  ignored  (that  is,  left  blank)  for  the  majority  of  cases. 

It  need  only  be  used  for  a downward  looking  path  (H2  < HI)  when  f«’o  paths  are 
pojsible  for  the  same  input  parameters  (see  Section  4,4).  In  suen  a case,  a com- 
puter printout  statement  will  be  given  indicating  that  the  user  has  two  choices  for 
t;ie  f lem  and  that  the  shorter  path  |c,ee  Figure  4(d)]  has  been  executed.  Set 
LEN  = 1 for  the  longer  case  (an  example  of  this  is  given  in  Section  6/  . 

J P = 0 for  normal  operation  of  program, 

.iP  = 1 to  suppress  printing  of  transmittance  table  (see  Table  1) 

JM  = 1 when  radiosonde  data  are  to  be  read  in  initially 

JM  = 0 for  normal  operation  of  orogram  or  when  subsequent  calculations 
are  to  be  run  with  MOD  El,  = 7 

MI,  = number  of  levels  to  b'-  read  in  for  MODEL  = 7 


Note  that  I VI  and  MI,  aio  only  used  when  MODEL  - 7 
on  tht  first  calculation  when  the  data  are  read  in. 

ill  — MO  — _ n ~ — 

i»  i J.  — i'i  »•  — **i  o — v i w » livsi  nidi  v/pc  » r\  i >CU  v»l  u 1 vt  i (tin, 


and  then  only 


The  parameters  Ml,  M2,  and  M3  can  each  take  integral  values  between  0 and  6 
and  are  used  to  modify  or  supplement  the  altitude  profiles  of  temperature,  water 
vapor,  and  ozone  respectively,  i'or  any  given  at.no  pheric  model  specified  by 
MODEL,.  For  example: 

Ml  = 1 selects  the  tropical  temperature  altitude  profile 

Ml  = 2 selects  the  MID1.AT1TUDE  SUMMER  temperature  altitude  profile' 

Ml  = 6 elects  the  1962  US  STANDARD  ternpe ’-ature  altitude  profile 


M2  = 1 selects  the  TROPICA i,  water  vapor  altitude  profile 

M2  = 2 selects  the  MIDI.ATITUDE  SUMMER  water  vapor  altitude  profile 

M2  = 6 selects  the  1962  US  STANDARD  wrter  vapor  altitude  profile 
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M3  = 1 selects  the  TROPICAL  ozone  altitude  profile 

M3  = 2 selects  the  MIDLATITUDE  SUMMER  ozone  altitude  profile 

M3  = 6 selects  the  1962  US  STANDARD  ozone  altitude  profile 

RO  = radius  of  the  earth  (km)  at  the  particular  geographical  location 
at  which  the  calculation  is  to  be  performed. 

If  RO  is  left  blank,  the  program  will  use  the  midlatitude  value  of  6371.23  km  if 
MODEL  is  set  equal  to  zero  or  7.  Otherwise  the  earth  radius  for  the  appropri  - 
ate standard  model  atmosphere  (specified  by  MODEL!  will  be  used. 

The  use  of  the  parameters  defined  above  will  become  more  apparent  by  re- 
ferring to  the  examples  given  in  Section  6. 

In  the  case  where  MODEL  = 7.  the  new  atmosphere  (model  or  radiosonde 
data)  is  inserted  between  CARDS  1 and  2 (see  Section  5.3). 

5.2.2  CARD  2 HU.  H2.  ANGLE.  RANGE,  BETA.  VIS 

Card  2 is  used  to  define  the  geometrical  path  parameters  for  a given  problem. 

HI  - initial  altitude  (km' 

H2  = final  altitude  (km) 

ANGLE  = initial  zenith  angle  (degrees'  as  measured  from  111 

RANGE  = path  length  (km' 

BETA  = earth  centre  angle  subtended  by  HI  and  H2  (degrees) 

VIS  = sea  level  visual  range  (km' 

It  is  not  necessarv  to  specify  every  quantity  given  above;  only  those  that 
adequately  describe  the  problem  according  to  the  parameter  ITYPE  (as  described 
below). 

(1)  Horizontal  Paths  (ITYPE  = 1) 

(a)  specify  HI,  RANGE  and  VIS  only 

(b)  if  nonstandard  meteorological  data  are  to  be  used,  that  is,  if 
MODEL  = 0 on  CARD  1),  then  the  following  parame'ers  must  be  specified  on 
CARD  2:  Hi,  p,  T,  DP.  RH,  V\  H , WO,  VIS,  RANGE  according  to  FORMAT 
(3F10.3.  2F5.1,  2E10.3,  2F10.3).  where  P.  T,  DP,  RH,  WH  and  WO  are  the 
pressure  (mb),  temperature  ( °Cr.  dew  point  temperature  ( “O,  relative  humidity 
(%),  HgO  density  (gm  m’"*)  and  ozone  density  (gm  m S)  respectively. 

Note  that  it  is  necessary  to  specify  all  of  the  quantities  underlined  with  a full 
line  and  one  of  the  quantities  underlined  with  a dashed  line.  If  the  ozone  density 
(vVO)  is  not  known,  a value  can  be  chosen  from  one  of  the  standard  atmospheric 
models  by  using  the  parameter  M3  on  CARD  1 (see  Section  5.  2.  1 above). 

Some  examples  of  typical  horizontal  path  calculations  are  given  in  Sections 
6,1  and  6.C , 
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(2)  Slant  Paths  to  Space  (ITYPE  = 3) 

(a!  specify  HI,  ANGLE  and  VIS 

(b)  specify  HI,  HMIN  and  VIS  (for  limb  viewing  problem  where  HMIN 
is  the  required  tangent  height  or  minimum  altitude  of  the  path  trajectory. 

<31  Slant  Paths  Between  Two  Altitudes  (ITYPE  = 21 
(a)  specify  HI.  H2.  ANGLE  and  VIS 
<b)  specify  HI,  ANGLE,  RANGE  and  VIS 
<c)  specify  HI,  H2,  RANGE  and  VIS 

For  cases  <b)  and  'cl,  the  program  will  calculate  H2  and  ANGLE  respec- 
tively. assuming  no  refraction;  then  proceed  as  for  case  (a).  This  method  of  de- 
fining the  problem  should  be  used  when  refraction  effects  are  not  important;  for 
example,  for  ranges  of  a few  tens  of  km  at  zenith  angles  less  than  80°.  It  can 
also  be  used  for  larger  angles  (including  90“)  provided  that  the  path  lies  within 
one  atmospheric  layer. 

<d>  Specify  HI.  H2,  BETA  and  VIS.  Leave  ANGLE  and  RANGE 
blank  in  this  case.  This  method  can  be  used  when  the  geometrical  configuration 
of  the  source  and  receiver  is  known  accurately,  but  the  Initial  zenith  angle  is  not 
known  precisely  uue  io  atmospheric  refraction  effects.  Beta  is  most  frequently 
determined  by  the  user  from  ground  range  information. 

In  the  cases  of  2(b)  and  3(d)  above,  the  subroutine  ANGLE  ia  called  in  the 
program  to  determine  the  appropriate  input  zenith  angle  by  an  iterative  technique 
taking  into  account  atmospheric  refraction  (see  Appendix  O. 

In  the  case  where  MODEL  = 7,  the  new  model  atmosphere  (or  radiosonde  data) 
is  inserted  between  CARDS  1 and  2 (see  Section  5.3). 

5.  2.  3 CARD  3 VI,  V2,  DV 

The  spectral  range  over  which  transmittance  data  ere  required  and  the  spec- 
trai  increments  **t  vhich  the  d&ts  are  to  be  printed  out  is  determine  - by  CARD  3. 

VI  = initial  frequency  in  wavenumbers  (cm 

V2  n final  frequency  in  wavenumbers  (cm~*)  where  V2  > VI 

_ , 

DV  = frequency  increment  (or  step  size'  (cm  ) 

(Note  that  v - 104/A  where  v is  the  frequency  in  cm-1  and  A is  the  wavelength 
in  microns,  and  that  DV  can  only  take  values  which  arc  a multiple  of  5. ) 

5.2.4  CARD  4 IXY 

The  control  parameter  IXY  can  cause  the  program  to  recycle,  so  that  a 
series  of  problems  can  be  run  with  one  submission  of  LOWTRAN  3.  Five  values 
of  IXY  can  be  used  to  provide  the  options  given  below. 

IXY  = 0 or  blank  card  to  end  of  program 

= 1 to  select  a new  CARD  3 and  CARD  4 only  (assuming  other  param- 
eters are  unchanged) 
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= 2 to  select  a new  data  sequence  (CARDS  1.  2,  3,  and  4) 

= 3 to  select  a new  CARD  2 and  CARD  4 only  (assuming  other  parameters 
are  unchanged) 

= 4 to  select  a new  CARD  1 and  CARD  4 only  (assuming  other  parameters 
are  unchanged) 

Thus,  if  for  the  same  model  atmosohere  and  type  of  atmobpher.c  path  the 
reader  wishes  to  make  further  transmittance  calculations  in  different  spectral 
intervals  VI1  to  V2'  etc.  and  lor  a different  step  size  (DV’  etc.),  then  IXY  is 
set  equal  to  1.  In  this  case,  the  card  sequence  is  as  follows  and  can  be  repeated 
as  many  times  as  required: 


CARD  4 

IXY 

= 1 

CARD  5 

VI' 

V2' 

DV 

CARD  6 

IXY 

= 1 

CARD  7 

VI" 

V2" 

DV 

CARD  8 

IXY 

= 0 

The  final  IXY  card  should  always  be  a blank  or  zero  (see  also  Section  6.2).  W hen 
using  the  IXY  = 1 option,  the  wavelength  dependence  of  the  refractive  index  is  not 
changed  (use  IXY  = 2 option  if  this  is  required). 

To  make  successive  transmittanc  ? computations  where  just  the  geographical 
n.cu*l  atmosphere  is  changed  and/or  with  or  without  aerosoi  attenuation,  set 
IXY  = 4 and  construct  a data  card  sequence  along  the  same  lines  as  given  above. 
This  sequence  of  recycling  can  be  repeated  successively,  and  several  examples 
are  given  in  Section  6, 

5.  2.  5 PROBLEM  SEQUENCING 

When  a series  of  problems  is  to  be  executed  (with  one  submission  of  LOW  - 
TRAN  3)  involving  the  standard  atmospheric  models  (MODEL  = I to  6)  as  well  as 
cases  involving  MODEL  = 0 and  MODEL  = 7,  then  the  order  in  which  the  data  are 
set  up  becomes  very  important.  Note  the  following  sequence. 

1.  Run  all  problems  using  MODEL  = 1 through  6 first. 

2.  Secondly,  run  ali  problems  involving  the  use  of  MODEL.  = 0. 

3.  Run  all  problems  involving  the  use  of  MODEL  = 7 last.  The  reason  for 
running  MODEL  = 7 cases  lust  is  that  when  a new  atmospheric  model  is  read  in, 
the  altitudes  may  not  correspond  with  those  given  in  the  standard  models  (see 
Section  2,1)  and  the  program  will  erase  them.  Similarly,  if  a MODEL  - 0 case  is 
run  following  a MODEL  = 7 case,  the  first  level  of  MODEL  7 is  erased, 

5.3  Non-Standard  Gundittom 

Three  options  are  available  if  atmospheric  transmittance  calculations  are  re- 
quired for  non-standard  conditions.  Here  non-standard  refers  ( conditions  other 
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than  those  specified  by  the  six  model  atmospheres  provided  by  LOWTRAN  3 (see 
Section  2.1),  which  ore  selected  by  the  parameter  MODEL  on  CARD  1 (see  Sec- 
tion 5.  2.  1).  The  three  options  enable  the  reader  to  insert: 

(1)  his  own  model  atmopshere(s)  in  place  of  any  (or  all)  of  the  six  standard 
models,  provided  that  the  data  are  in  exactly  the  same  format  and  are  specified 
at  the  same  altitudes  as  the  latter.  In  this  case  the  appropriate  print  statements 
in  LOWTRAN  3 (that  identify  ths  atmospheric  model  used)  must  be  changed  cor- 
respondingly. 

(2)  an  additional  atmospheric  model  (MODEL  7),  which  can  be  in  the  form  of 
radiosonde  data.  The  data  need  not  be  specified  at  the  same  altitudes  as  in  the 
standard  models. 

(3)  meteorological  conditions  for  a given  horizontal  path  calculation  (MODEL  - 
0 case). 

The  first  of  these  options  requires  the  most  effort  and  needs  no  further  dis- 
cussion here,  other  than  a reference  tc  Appendix  A for  a summary  of  the  stan- 
dard model  atmosphere  parameters.  i nits  and  formats  (see  also  Table  A 2). 

5.  3.  1 ADDITIONAL  ATMOSPHERIC  MODEL  (MODEL  = 7) 

As  stated  in  Section  5.  2.  2 a new  model  atmosphere  can  be  inserted  between 
CAHUS  1 and  2 provided  the  parameters  MODEL  and  IM  are  set  equal  to  7 and  1 
respectively  on  CARD  1.  The  number  of  atmospheric  levels  to  be  inserted  (ML) 
must  also  be  specified  on  CARD  1.  The  appropriate  meteorological  parameters 
and  format  for  the  atmospheric  data  are  given  below, 

2,  P.  T.  DP.  RH.  WIL  WO.  AHA2E  (FORMAT  (3F10.3,  2F5.1,  2E10.3, 
2FI0,3>]  where 

Z = altitude  (km) 

P = pressure 

T = ambient  temperature  (°C) 

DF  = dew  point  temperature  <“C) 

RH  = relative  humidity  (%) 

_ 3 

WH  = water  vapor  density  (gm  m > 

WO  = ozone  density  (gm  rn  a) 

AHAZE  = aerosol  number  density  (cm  ) 

Note  that  it  is  only  necessary  to  specify  those  quantities  underlined  with  a full 
line  and  either  of  the  quantities  underlined  with  the  dashed  line. 

If  the  ozone  density  (WO)  ia  not  known  then  a value  can  be  obtained  from  one 
Of  the  standard  atmospheric  models  (for  the  appropriate  latitude  and  season)  by 
using  the  parameter  M3  on  CARD  1 (see  Section  5.  2.  1). 

If  the  aerosol  number  density  was  not  measured  as  a function  of  altitude  and 
the  reader  wishes  to  include  aerosol  attenuation  in  the  calculation,  set  1HAZE  * 1 
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on  CARD  1.  In  this  case  (as  with  the  Ml.  N2,  and  M3  options)  LOWTRAN  3 will 
use  the  aerosol  modelo  already  contained  in  the  program  and  interpolate  to  give 
aerosol  number  density  values  at  the  same  altitudes  as  the  radiosonde  (or  new 
model  atmosphere)  data.f  The  program  wi'i  then  look  for  a sea  level  visual 
range  (VIS)  to  be  specified  on  CARD  2.  If  VIS  is  not  specified,  a 23  km  sea  level 
visual  range  will  be  assumed.  If  aerosol  attenuation  is  not  required,  set  IHAZE  = 0 
on  CARD  1 as  before. 

5.  3.  2 HORIZONTAL  PATHS  (MODEL  = 0) 

If  meteorological  data  are  to  be  used  for  horizontal  path  atmospheric  trans- 
mittance calculations,  then  set  MODEL  = 0 on  CARD  1,  The  following  param  - 
eters can  then  be  specified  on  CARD  2: 

CARD  2 lU,  P,  T.  DP,  RH,  WH,  WO,  VIS.  RANGE  (FORMAT  (3F10.3, 
2F5.1,  2E10.3.  2F10.3))  where  the  above  parameters  refer  to  altitude  (km), 

pressure  (mb)  ambient  temperature  (cC),t  dew  point  temperature  (eC),  relative 

-3  -3 

humidity  (%).  water  vapor  density  (gm  m ),  ozone  density  (gm  m ‘ ),  visual 
range  (km)  and  path  length  (km)  respectively  (.is  previously  defined  in  Sections 
5.3.1  and  5.2.2). 

The  format  for  the  above  card  is  similar  to  that  for  inputting  radiosonde  data 
(MODEL  = 7)  described  in  Section  5,  3.  1 above.  Again,  it  is  only  necessary  to 
specify  the  quantities  underlined  with  the  solid  line  and  one  of  the  quantities  under- 
lined with  the  dashed  line.  The  ozone  density  WO  can  be  specified  using  the  pa- 
rameter M3  on  CARD  1 if  measurements  are  not  available.  In  the  latter  case,  a 
value  will  be  calculated  at  altitude  HI  based  on  the  appropriate  model  atniob^-here 
selected  by  M3. 

6.  EXAMPLES 

The  following  examples  explain  the  problem  parameters  that  need  to  be  speci- 
fied and  card  formats  necessary  for  using  LOWTRAN  3 to  make  a wide  range  of 
atmospheric  transmittance  calculations.  Example  9 describes  how  the  reader 
may  insert  an  additional  model  atmosphere  (for  example,  a set  of  radiosonde  or 
rawinsonde  data)  into  the  program. 

As  mentioned  earlier,  only  the  last  four  data  cards  (see  Table  A2  in  Appen- 
dix A)  are  necessary  to  specify  a given  problem.  It  should  be  noted  that  the 

t A similar  interpolation  scheme  is  applied  when  the  Ml,  M2  and  M3  options  are 
used  to  supplement  MODEL  7 (and  also  for  MODEL  = 0)  using  the  standard 
atmospheric  models. 

t Note  that  temperature  is  given  in  °C  for  MODEL  = 0 option  and  not  ”K  as  pre- 
viously used  in  LOWTRAN  2, 
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terms  ''CARD  1,  CARD  2,  etc."  are  used  for  convenience  to  denote  these  impor- 
tant control  cards  in  their  correct  sequence,  and  do  not  refer  to  t^e  actual  first, 
second,  etc.  data  cards  in  LOWTRAN  3 (described  in  Appendix  Ah 

6.1  Luuuptr  1 - Horizontal  Path  Calculation!  (ITYPE  = 1) 

PROBLEM  AND  CARD  3RMAT 

Calculate  the  transmittance  from  500  cm  * to  1000  cm  * in  steps  of  5 cm  1 
for  a horizontal  path  of  range  10  km  at  an  altitude  of  2.5  krn.  for  a midlatitude 
summer  model  atmosphere  and  the  5 km  (sea  level)  visual  range  haze  model. 
CARD  1 **2*v*i 

CARD  2 ***  ■•••-*2. 500  **''  ****  -****'-■':  : : :=  * : *<-10.000  :***  :■*  ******** 5.000 
CARD  3 --*500.000  -1000.000-  **: *5.000 

CARD  4 *** 

or  alternatively 

CARD  1 * - 2 * - 2 - * 1 

CARD  2 *****2.500  —'  ********  : *****  ***  ****10,0 

CARD  3 ***500.000  1000.000  ***5.0 

CARD  4 Blank 

{'  represents  a space  on  the  card.) 

Note  that  setting  (HAZE  = 2 can  select  a 5 km  visual  range  at  sea  level  without 
setting  VIS  on  CARD  2. 

6.2  Extunpir  2 - Sant  P«th  Calculation  to  Space  (ITYPE  = 3)  for  Two  Spectral  Interval* 

PROBLEM  AND  CARD  I ORMAT 

Calculate  the  transmittance  from  3 to  5 fJ.cn  and  8 to  14  fim  in  steps  of 
~0.05prn  for  a 45°  slant  path  to  space  from  12.2  km,  for  a midlatitude  winter 
model  atmosphere  and  a 23  km  visual  range  haze  model, 

lii  this  case  Vi  = 2000  cm"1,  V2  = 3335  cm  and  DV  varies  from  20  cm  1 
to  55.5  cm  *,  Let  us  choose  the  lower  value,  that  is,  DV  = 20  cm  S also, 

VT  = 714  cm1  and  V21  = 1250  cm  1 with  DV  varying  from  2.55  cm  1 to 
7,81  cm"1.  Thus  DV'  = 5 cm"1,  since  this  is  the  nearest  multiple  of  5 cm 
CARD  1 **3 ' 13 

CARD  2 * * * > 1 2. 2 00 * * * : ************ 4 5 , 0 

CARD  3 **2000.000  *3335.000  ' ****20.0 

CARD  4 <*1 

CARD  5 ***714.000  : 1250.000  *****5.00 

CARD  6 Blank 

Note  setting  IHAZE  - 1 on  CARD  1 will  select  a 23  km  sea  level  visual  range  if 
VIS  is  not  set  on  CARD  2. 
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6.3  Example  3 - Upward  Looking  Sant  Path  Calculation  (ITYPE  - 2) 

PROBLEM  AND  CARD  FORMAT 

Calculate  the  transmittance  from  4 to  5jtm  in  steps  of  *'0,01  pi m for  a slant 
path  from  2.25  km  to  22.8  km  at  zenith  angle  of  65°,  for  z subarctic  winter  model 
atmosphere  and  a 15  km  (sea  level!  visual  range  haze  model. 

In  this  case  VI  = 2000  cm  V2  = 2500  cm  *.  and  DV  varies  from  4 cm  1 

to  6.26  cm  '.  Therefore,  set  DV  = 5 cm  ' since  this  is  the  nearest  multiple  of 

, *1 
5 cm 

CARD  1 ' ' 5 : 1**2 

CARD  2 ******2.25  ••*****8.50 ' :****65.0‘ ***********************15.000 

CARD  3 **2000.000 ••'*2500.000: ******5.0 

CARD  4 Blank 


6.4  Example  4 - Slant  Path  Hrtween  Two  Altitude*  tanking  Uownwaid  (ITYPK  = 2) 


PROBLEM  AND  CARD  FORMAT 


Suppose  the  transmittance  from  10  to  20/im  is  required  for  a downward  path 


i .5  .-H.muuiu  inouei 


Atmosphere  and  no  haze. 

In  this  case  set  VI  = 500  cm  1 <20Mm',  V2  = 1000  cm  5 (lOpm)  and  DV  = 
5 cm  *. 

CARD  1 *6  *0  2 

CARD  2 ***•■  10.000  ***8.000- ***92.000 

CARD  3 * * *500.000  * 1000.000  "■■■"•  *5.000 

CARD  4 Blank 


When  this  case  is  executed,  however,  there  will  be  a message  printed  out  tc 
the  effect  that  a longer  path  is  also  possible  for  the  same  input,  conditions  [for 
example,  see  Figures  3(d>  and  3(e)].  If  you  wish  to  rerun  the  conditions  for  the 
longer  path,  the  card  sequence  is  as  follows; 

CARD  1 **6'r*0"*2'*l 

CARD  2 **  + <10.000- ****8.000  ***92.000 

CARD  3 ***500.000 : * 1000.000  : * + **5.000 

CARD  4 Blank 


6.5  Example  5 - Repeated  Problem*  in  Sequence 

PROBLEM  AND  CARD  FORMAT 

Calculate  the  transmittance  for  the  4 to  Bjim  region  for  <1>  a 15°  zenith 
angle  slant  path  from  0.23  km  to  8.55  km  for  tropical,  midlatitude  summer  and 
subarctic  summer  conditions,  and  (2)  for  horizontal  paths  of  1.  5.  and  10  km  at 
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5 km  altitude  for  midlatitude  summer  conditions.  In  both  cas-'.  , assume  the  sea 
level  visual  range  to  be  15  km. 

The  set  up  cards  for  the  above  problems  are  as  follows: 


CARD  1 
CARD  2 
CARD  3 
CARD  4 
CARD  5 
CARD  6 
CARD  7 
CARD  8 
CARD  9 
CARD  10 
CARD  11 
CARD  12 
CARD  13 
CARD  14 
CARD  15 
CARD  lfi 


* * 1 ^ * 1 ' « 2 

*#*  **0.230**’**  *8.550**  * *7  5, 000’  ************  *****  ******  j5<  Q00 
**1250.000" *2500.000*****5.000 

* *4 

**2**1" *2 
*'*4 

* *4 **1**2 
<■*2*  *1**1 

* $ * *5,000**  * ************** ********1, 000^ 
**1250.000**2500.000*****5.000 

**3 

* * # > * *5.000** *********************  **5,000^ 

* *3 

^*<.**5000******  ********  **********'10.  OOOt 

*** 


6.6  Example  6 - To  Calculate  Viewsig  Angle,  Taking  Into  Account  the  Effect  of  Refraction 

PROBLEM  AND  CARD  FORMAT 

Suppose  the  exact  position  of  a source  and  receiver  is  known;  that  is.  HI, 
H2.  and  geometrical  range.  From  this,  one  can  calculate  the  apparent  zenith 
angle  (assuming  no  refraction)  and  the  total  angle  aubtendej  by  the  path  at  the 
center  of  the  earth,  P (see  Figure  3).,  It  the  apparent  zenith  angle  is  such  that 
atmospheric  refraction  could  cause  appreciable  bending  of  the  path  trajectory, 
then  the  receiver  could  be  considerably  off  target.  To  take  into  account  the  effe." 
of  atmospheric  refraction,  one  can  use  the  subroutine  ANGLE  (with  HI,  H2,  and 
BETA  as  inputs)  which  is  called  by  the  program  to  determine  the  correct  zenith 
angle. 

Suppoee  Example  4 was  the  above  case  where  = 5.3715°  is  known  specifi- 
cally for  the  shorter  path.  The  card  sequence  would  then  be: 

CARD  1 **6**0**2 

CARD  2 ** **10.000*****8.000* *************** ********5. 37  15 

CARD  3 ***500. 000**1000. 000**+**5. 000 

CARD  4 Blank 

t Since  the  aerosol  attenuation  models  are  independent  of  sea  level  visual  range 
above  5 km  altitude,  it  is  not  necessary  to  specify  VIS  in  this  case. 
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6.7  Exmplr  7 - Earth  Limb  Viewing  Cut 

PROBLEM  AND  CARD  FORMAT 

Consider  an  earth  limb  viewing  problem  where  it  is  required  to  calculate  the 
atmospheric  transmittance  from,  say,  an  altitude  of  100  km  to  space  passing 
through  a tangent  height  (HMIN)  of  12  km  (see  Figure  4'  , ..  the  2 to  5 pm  region 
for  a subarctic  winter  model  atmosphere.  In  this  case  let  us  assume  that  aerosol 
attenuation  is  required.  As  described  in  Section  5 , 2.  2 , only  the  quantities  III 
and  HMIN  are  specified  on  CARD  2.  The  appropriate  zenith  angle  (taking  into 
account  atmospheric  refraction)  is  calculated  by  the  program  using  the  subroutine 
ANCL. 

CARD  1 *5  I ' <3 

CARD  2 ***100.000  **  *1  2,000 

CARD  3 **2000.000  *2500.000  ”**5.000 

CARD  4 *** 

6.8  Example  8 - Horizontal  Path  with  hon-Standanl  Condition* 

If  the  appropriate  conditions  are  known  for  the  horizontal  path,  then  a differ- 
ent procedure  to  that  given  in  Examples  i to  7 can  be  used,  t or  this  case, 

CARD  2 (see  Sections  5.1  and  5.21  contains  the  following  information:^ 
altitude  (km),  pressure  (mb),  ambient  temperature  (°C),  dew  point  temperature 
CC).  relative  humidity  i^),  water  vapor  density  (gm.m  ),  ozone  density  (gm.m  ), 
sea  level  visual  range  (km),  and  path  length  (km)  according  to  format  '31  10.3, 
2F5.2,  2E10.3.  2F10.31. 

PROBLEM  AND  CARD  FORMAT 

Calculate  the  transmittance  from  3 to  5pm  for  a 10  km  path  at  sea  level 
(midlatitude  winter  environment)  for  the  following  conditions:  Pressure  = 1000  mb, 
ambient  temperature  - I0"C,  relative  humidity  = 40To,  and  sea  level  visual  range  = 
50  kin. 

CARD  1 ■ >0  ' *1 1 ' *0  0 *0  *0- *0  - 3 

CARD  2 *****0.000  *1000.000  *00.000  **40,00  ' ** '^  ***■  *****';-****<•  : *^ 

50.000  *‘-*10.000 

CARD  3 **2000.000  *3330.000  ****5.000 

CARD  4 *** 

The  index  3 on  CARD  1 corresponds  to  the  parameter  M3  (see  Section  5.21 
and  selects  the  ozone  density  corresponding  to  the  midlatitude  winter  model  atmo- 
sphere, since  the  latter  quantity  is  not  specified  in  the  above  problem.  If  M3 

t Note  that  this  problem  is  specified  differently  in  I.OWTRA  2. 
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were  left  blank,  then  ozone  attenuation  would  not  be  included  in  the  transmittance 
calculation. 

If  a number  of  different  problems  are  being  run  with  one  submission  of  the 
program  involving  the  standard  model  atmospheres  as  well  as  problems  similar 
to  Example  6.8,  then  the  order  in  which  the  problems  are  run  becomes  important 
(see  Section  5.3). 

6.9  Lump!?  9 - To  faiwrl  or  Change  Model  Aimoephere  Data 

If  the  reader  wishes  to  include  his  own  model  u*mosphere  data  in  l.OWTRAN  3 
two  options  are  available.  Either  replace  one  or  more  of  the  standard  models  by 
the  new  data  or  use  the  MODEL  = 7 option  which  is  more  convenient.  If  the  for- 
mer approach  is  taken,  it  is  recommended  that  the  new  atmospheric  data  be  digi- 
tized at  the  same  altitudes  and  according  to  the  same  format  as  the  remaining 
standard  models  (see  Appendix  A and  Table  A2),  Note  that  each  of  the  standard 

model  atmosphere  data  cards  contain  the  following-  altitude  'km',  pressure  (mb), 

* 3 - 3 

temperature  ("K),  water  vapor  density  (gm. in  ),  and  ozone  density  (gm.m  ) 

with  two  models  specified  or.  one  card  (see  Table  A?!  according  to  format  (K6. 1, 

2E10.3.  F6.1.  2E10.7). 

In  the  above  case,  the  reader  should  modify  the  appropriate  print  statements 
in  the  program  accordingly  in  order  to  identify  the  new  model  atmospheres,  other- 
wise the  printout  from  L.OWTRAN  3 will  identity  them  according  to  Me  standard 
atmospheric  models  which  they  replaced. 

If  the  reader  wishes  to  make  a limited  number  of  atmospheric  transmittance 
calculations  for  one  or  more  sets  of  radiosonde  (or  rawinsonde)  data,  it  is  more 
convenient  to  use  the  MODEL  = 7 option;  an  example  of  which  is  given  below. 

PROBLEM  AND  CARD  FORMAT 

Calculate  for  a given  set  of  radiosonde  data  the  transmittance  in  the  4 to 
8pm  region  for  (1>  a slant  path  from  0,21  km  to  8.55  km  at  a zenith  angle  of 
35, S'  and  (2'  horizontal  paths  of  i,  5,  and  10  km  at  0.21  km  altitude.  Assume 
the  sea  level  visual  range  to  be  15  km  and  the  ozone  distribution  to  be  represen- 
tative of  midlatitude  summer  conditions.  In  this  example,  the  radiosonde  data 
consists  of  21  levels  and  only  the  following  parameters  are  given:  altitude  (km), 
pressure  (mb',  ambient  temperature  (°C),  and  dew  point  temperature  <°C). 

The  set  up  cards  for  the  above  problem  are  as  follows: 

• ft ##*0.000* 4 10 15.000* ***24.4  '‘** *21.4 

0.136  1000.000  22.0  10.4 

0.560  950.000  17.8  16.1 
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1.080 

892.000 

14,8 

1 1.9 

1.5  20 

850.000 

12.8 

5.8 

1.650 

832.000 

12.8 

-6.2 

2.270 

77  5.000 

11.8 

-18.2 

3. 140 

700.000 

7.2 

-20.8 

5.820 

500.000 

-10.1 

-28, 1 

5.900 

488.000 

- 1 1.5 

-27.5 

7.510 

400.000 

-19.5 

-31.5 

8.720 

338.000 

-28.5 

-4  1.5 

9. 180 

318.000 

-32.7 

-39.7 

9.590 

300.000 

-35.3 

-43.3 

9.720 

294,000 

-34.7 

-42.7 

10.020 

281.000 

-38.7 

-45.7 

10.330 

250.000 

-44.7 

-50.0 

12.290 

200.000 

-57.1 

-50.0 

13.600 

161.000 

-69.5 

-50.0 

14.050 

150.000 

-71.1 

-50.0 

16.450 

100.000 

-70.9 

-50.0 

**»*»y. 210«***»u. 550 ‘*+*3 5. 500 4 ****** ******  -t ******** *15^ 000 
' 4 1 250.000  2500.000  <5.000 

* 2 

‘<7  110  0 O 0 0 0 21 


*‘*•■•0.210  ' •■  '•‘1.000  - : ‘ ~'J-  '■  15,000 

■■'1250.000  2500.000  '5.000 

' *3 

***  -0.210  '5.000  ‘ ‘“15.000 

* < 3 

-tf.  ',0.2  10  1 »-  1 0.000  -'■*  ‘ ••'•'•****•*15.000 


Note  that  problems  utilizing  the  standard  modt1  atmospheres  must  precede 
problems  similar  to  Example  0 (sec  Section  5,3>. 

6.10  Tjpifil  Output  frum  LO^TRAN  3 

The  output  for  a problem  similar  to  Example  3 (.Section  0.3)  is  given  in 
Table  1.  The  parameters  defining  the  atmospheric  path,  mode!  atmospheres 
and  frequency  range  will  first  tie  printed  out.  Following  the  heading  HORI/.ONTA1 
PROFILES  there  are  12  columns.  The  first  column  gives  a running  integer  asso- 
ciated with  each  level  (level  indicator).  The  second  column  gives  the  level  altitude 
in  km.  The  next  8 columns  give  the  equivalent  absorber  amounts  per  km  for  the 
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following  absorbing  species:  water  vapor,  uniformly  mixed  gases,  o/one,  nitro- 
gen continuum,  water  vapor  continuum,  molecular  scattering,  aerosol  extinction, 
and  L'V  ozone,  respectively.  The  last  two  columns  give  the  mean  refractive  index 
modulus  from  that  level  to  the  level  above  and  the  refractive  index  modulus  (mul- 
tiplied bv  lb®)  at  that  level. 

A heading  VERTICAL  PROFILES  is  then  printed  followed  by  ,5  co  umns. 

The  first  and  second  columns  give  the  integer  associated  with  the  leve'a  traversed 
by  the  path  and  the  height  of  the  level.  Then  follow  8 columns  which  give  the  in- 
tegrated equivalent  absorber  amounts  from  the  initial  altitude  to  the  level  above 
(in  the  same  order  as  indicated  above!.  The  next  4 columns  are  labelled  I’SI, 

PHI,  BETA,  and  THETA,  and  correspond  to  the  angles  similar  to  4,  d,  0,  and  9 
given  in  Figure  3 (Section  4.2).  Columns  PSI  and  BETA  give  the  accumulated 
values  of  V and  0 to  the  level  above.  Columns  THETA  and  PHI  give  the  local 
zenith  angle  (V  corresponding  to  that  level  and  the  angle  of  arrival  at  the  level 
above,  respectively.  The  accumulated  slant  range  is  printed  out  in  the  last  col- 
umn under  RANGE. 

The  total  equivalent  absorber  amounts  for  each  absorber  species  are  then 
summarized  below  in  their  appropriate  units, 

A transmittance  table,  containing  12  columns,  now  follows.  The  first  3 col- 
umns give  the  frequency  (cm  1 ) , wavelength  (microns),  ana  total  transmittance. 
The  next  7 columns  show  trie  individual  transmittance  due  to  water  vspor,  um- 
i-  rmly  mixed  gases,  ozone,  nitrogen  (4  pm!  continuum,  water  vapor  (10pm)  con- 
tinuum. molecular  scattering,  and  aerosol  extinction.  The  last  2 columns  give 
the  absorption  due  to  aerosols  and  the  cumulative  total  integrated  absorption. 

The  latter  quantity  can  be  used  to  determine  the  average  transmittance  over  any- 
given  spectral  interval  within  the  spectral  range  covered  by  the  calculation. 
Finally,  the  total  integrated  absorption  from  Vl  to  V2  is  printed  out  (units  are 
cm  S together  with  the  average  transmittance  "ver  the  band. 

7.  LIMITATION  AMD  COMMENTS 

It  should  be  remembered  that  the  transmittance  values  obtained  from 
LOWTRAN  3 are  at  a spectral  resolution  of  20  cm”  , although  the  output  can  be 
obtained  at  5 cm  * intervals. 

The  program  will  round  off  input  frequencies  to  the  nearest  frequency  at  which 
spectral  data  are  given. 

The  overall  accuracy  in  transmittance,  which  this  technique  provides,  is 
better  than  10  percent.  The  largest  errors  may  occur  in  the  distant  wings  of 
strongly  absorbing  bands  in  regions  where  such  bands  overlap  appreciably. 
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Table  1.  Typical  Output  From  I.OWTRAN  3 
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The  reason  for  this  error  is  twofold.  First,  a unique  spectral  .urve  in  Figures 

2 

13  to  18,  McClatchey  et  al,  is  based  on  a single  absorber  parameter  and  cannot 
be  defined  for  a wide  range  of  atmospheric  paths  without  some  loss  in  accuracy. 

Secondly,  xhe  transmittance  in  me  window  regions  between  strong  bands  gen- 
erally lies  in  the  weak  ,;ne  approximation  region,  where  the  transmittance  is  a 
function  of  the  quantity  of  absorber  present  and  not  of  the  product  of  absorber 
amount  and  pressure.  The  one -dimensional  prediction  scheme  presented  in  this 
report  is  not  accurate  for  such  conditions,  which  in  general  give  transmittance 
values  greater  than  0.  99.  The  digitized  data  were  obtained  for  conditions  repre- 
sentative of  moderate  atmospheric  paihs  and  will  tend  to  overestimate  the  trans- 
mittance for  very  long  paths  and  underestimate  the  transmittance  for  very  short 
paths,  in  the  spectral  regions  described  above. 

As  the  transmittance  approaches  1.  0,  the  percentage  error  in  transmittance 
decreases  toward  zero  but  the  uncertainty  in  the  absorptance  (or  emittancto  in- 
creases greatly. 

Because  of  the  nature  of  the  program  - which  uses  a layered  atmosphere  ■■ 
errors  can  be  introduced  into  the  refraction  calculation,  since  we  assume  each 
layer  to  have  a mean  refractive  index  associated  with  it.  This  is  particularly  true 
for  a long  path  in  one  layer  near  ground  level  where  one  would  expect  refraction  to 
be  a maximum.  But  ti  fact,  for  such  a condition  the  program  may  indicate  no 
refraction  at  all.  If  problems  like  these  are  encountered,  the  number'  of  levels 
must  be  increased  in  the  altitude  region  of  interest. 

The  effect  of  wavelength  on  refraction  has  been  included  (see  Section  4.3). 

When  a particular  problem  is  set  up  for  a frequency  range  VI  to  V2,  however,  the 
average  frequency  is  used  in  the  refractive  index  calculation.  If  the  reader  is 
using  LOWTRAN  3 for  trajectory  analysis  where  wavelength  dependence  of  re- 
fractive index  may  be  important,  it  is  recommended  that  dummy  values  for  VI  and 
V2  be  used  such  that  the  average  frequency  corresponds  to  l he  frequency  of  interest. 

If  LOWTRAN  3 is  to  be  used  for  trajectory  analysis,  it  is  recommended  that 
the  number  of  levels  in  the  model  atmospheres  be  increased  substantially  between 
the  altitudes  of  interest  The  computer  program  (with  34  levels)  will  tend  to  under- 
estimate the  refraction  of  the  atmosphere  by  as  much  as  20  percent  for  the  worst 
c ise  (that  is,  a 90c  path  to  space  from  sea  level). 

If  LOWTRAN  3 is  not  heitig  used  for  problems  similar  to  Examples  6 and  7 it 
is  recommended  that  the  subroutine  ANGL  be  removed  from  the  program,  there- 
by reducing  the  size  of  the  card  deck  by  approximately  one -third. 

Some  examples  of  transmittance  spectra  obtained  from  LOWTRAN  3 are  pre- 
sented in  Figures  5 through  11.  From  Figure  5 one  can  see  the  separate  contri- 
butions to  the  total  transmittance  due  to  the  various  absorbing  gases,  molecular 
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scattering  and  aeroaoi  attenuation.  The  other  figures  show  the  variation  of  at- 
mospheric transmittance  with  (1)  model  atmosphere  for  fixed  paths,  (2)  altitude 
and  path  for  a fixed  model  atmosphere,  and  (3)  range  and  2enith  angle  for  a fixed 
altitude.  In  Figures  6 through  11  the  aerosol  attenuation  is  calculated  on  the  basis 
of  a 23  km  sea  level  visual  range. 

In  general,  fairly  good  agreement  has  been  found  between  experimental  field 
measurements  and  LOWTRAN  3 predictions  although  a lack  of  accurate  meteoro- 
logical data  (including  visual  range  and  actual  aerosol  characteristics)  gives  rise 
to  some  uncertainty. 


WAVELENGTH  ini  CRONS) 


Figure  5.  Vertical  Path  to  Space  From  Sea  Level  for  the  Midlatitude  Winter 
Model  Atmosphere  Showing  the  Separate  Contributions  to  the  Total  Trans- 
mittance 
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Atmospheric  Transmittance  for  a 1 km  Path  at  Sea  Level  for  Six  Model  Atmospheres 
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•e  9.  Variation  of  Transmittance  With  Altitude  for  a Vertical  Path  to  Space  for  the  1982 
Standard  Atmosphere 
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& COMPARISONS  OF  LOWTRAN  3 PREDICTIONS  WITH  MEASUREMENTS 

Figures  12  - 14  show  some  comparison  of  LOWTRAN  3 predictions  with  labors- 

25  26 

tory  measurements  of  Howard,  Burch  and  Williams  and  Burch  ° for  some  im- 
portant water  vapor  and  carbon  dioxide  bands.  It  will  be  seen  that  the  LOWTRAN 
3 calcvla'ions  agree  closely  both  spectrally  and  in  integrated  absorption. 

In  Table  2 a comparison  is  made  between  five  integrated  absorption  measure- 
ments of  HBW^  in  the  2.7  pm  H^O  band  and  those  predicted  using  LOWTRAN  3, 
The  measurements  were  chosen  to  be  representative  of  conditions  occurring  in  the 
first  50  km  of  the  atmosphere. 


4 

if 
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Table  2.  Comparison  of  LOWTRAN  3 Predicted  and  Measured  Integrated 
Absorptions  in  the  2.7  pm  Water  Vapor  Band 


A 

= Integrated  Absorption  (3450  - 3850  cm  *) 

— 

Pressure 

(mb) 

h2o 

(pr  cm) 

ALOWTRAN  3 
(cm-1) 

AMeas. 

Icm'*) 

Case  No,  1 

1073 

1.68 

399 

399 

Case  No.  2 

997 

0.21 

382 

374 

Case  No.  3 

4 13 

0.034 

27  3 

258 

Case  No.  4 

413 

0.017 

133 

131 

Case  No.  5 

90.6 

0.017 

226 

206 

27 

Figure  15  shows  a transmittance  spectrum  measured  by  Gebbie  et  al  for  a 
1 nautical  mile  path  over  water,  covering  the  0,  5 to  15  pm  region.  In  the  LOW- 
TRAN 3 simulation  of  the  spectrum  the  visual  range  was  assumed  to  be  13  km 
which  gives  a transmittance  of  60  percent  at  0.  6 pm.  The  effect  of  ozone  absorp- 
tion in  the  10  pm  region  can  be  seen  in  the  LOWTRAN  3 comparison,  which  shows 
a midlatitude  winter  concentration  in  one  case  and  no  ozone  in  the  other.  The  no- 
ozone case  appears  to  agree  better  with  Gebbie's  measurements. 


25,  Howard,  J.  N. . Burch,  D.  L. , and  Williams,  D.  (1955)  Near-Infrared  Trans- 

mission Through  Synthetic  Atmospheres,  AF'CRL-TR-S3-213,  Geophysical 
Research  Pipers  No.  40. 

26,  Burch.  D.  E. , Gryvnak,  D. , Singleton,  E.  B. , France,  W.  L. , and 

Williams.  D.  (1962)  Jnfrared  Absorption  by  Carbon  Dioxide,  Water  Vapor. 
and  Minor  Atmospheric  Constituents,  AFcRC-62-59f. 

27,  Gebbie.  H.A.,  Harding.  W.  Rv . Hilsum,  C, , Piyce,  A.W.,  Roberts,  V,  ( 195 \) 

Proc.  Roy.  Soc.  206 A: 87 . 
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Figure  15.  Atmospheric  Transmittance  for  a 1 Nautical 
Mile  Path  (Water  Content  1.7  pr  cm) 

2g 

Figures  16  -22  show  a series  of  measurements  by  Yates  and  Taylor  cover- 
ing a wide  spectral  interval.  The  path  lengths  in  these  measurements  are  0.  3,  5.5, 
16.25  and  27.7  km  and  were  made  at  sea  level  (over  water)  with  the  exception  of 
the  27.7  km  path,  which  w«J  at  an  altitude  of  10,^00  feet.  No  aerosol  attenuation 
was  used  in  simulating  trie  0.3  urn  path  using  LOWTHAN  3.  The  apparent  wave- 
length shift  at  short  wavelengths  in  Figure  16  is  due  to  a calibration  error  in  Yates 
and  Taylor's  measurements.  The  obvious  discrepancy  between  prediction  and 
measurement  in  the  10  un  is  due  to  the  fact  that  5' ales  and  Taylor  artificially  set 
the  transmittance  level  to  be  100  percent  in  this  window  region  (since  they  were 
unable  to  estimate  the  water  vapor  continuum  contribution).  The  water  vapor  con- 
tinuum given  in  LOWTRAN  3 extends  from  7 to  15  u m.  There  is  evidence  that  the 

12 

effect  of  continuum  absorption  bv  H,0  extends  to  longer  wavelengths  (Bignell, 

13  L 

Burch  et  al  ‘ ),  and  this  is  currently  being  investigated. 

28,  Yates,  H.  W. , Tavlor.  H.  (1060)  Infrared  Transmission  of  the  Atmosphere, 
NR  Report  5453,  U.  S.  Naval  Research  Laboratory,  Washington.  IJ.  C. 
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Figure  18.  Comparison  of  I/WTRAN  3 Predictions  With  Atmospheric  Transmittance  Measurements  Over: 
(at  5.5  km  ; (b)  16.25  km  in  the  Chesapeake  Bay  Area  on  19  June  1956 


Figure  21.  Comparison  of  IJOWTRAN  3 Predictions  With  Atmospheric  Transmittance  Measurements  Over 
27.7  km  in  the  Hawaiian  Islands  on  1 September  1957 
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’n  Figures  17  - 22  the  transmittance  values  at  0.55  pm  given  by  Yates  and 
Taylor  were  used  to  estimate  the  visual  range  used  in  the  LOWTRAN  3 simulations. 
In  some  cases  (see  Figure  22  for  example),  choosing  a lower  visual  range  gives 
far  better  agreement  between  the  measurements  and  LOWTRAN  3 predictions. 
Generally,  the  predictions  and  measurements  disagree  moat  in  the  10  pm  region, 
where  it  is  extremely  difficult  to  obtain  absolute  transmittance  values  from  long 
path  measurements  in  the  atmosphere.  Because  of  the  importance  of  the  10  pm 
window  region  in  many  applications,  it  is  strongly  suggested  that  future  field  mea- 
surements of  this  kind  be  carried  out  with  extreme  care,  in  order  to  obtain  reliable 
atmospheric  transmittance  data  against  which  a model  such  as  LOWTRAN  3 can  be 
tested. 

It  is  also  apparent  that  there  are  some  inconsistencies  in  the  measurements 
shown  in  Figures  16  - 22  which  are  the  cause  of  some  apparent  deviations  between 
the  LOWTRAN  3 predictions,  and  the  measurements. 

Figures  23  - 26  show  some  more  recent  sea  level  measurements  made  by 
2d 

Ashley  et  al  (General  Dynamics),  These  measurements  have  been  made  at  some- 

- 1 

what  higher  resolution  (—4  cm  1 and  cover  the  spectral  regions  1.8  -5.4  um  and 
4.8-14  pm  using  two  interferometers  with  different  detectors.  In  Figures  23  and 
24  the  transmittance  measurements  were  made  over  short  and  medium  path  lengths 
(that  is,  "'U.045  km  and  3.25  k*n)  and  covered  the  Spectral  region  frGiii  1.8  iti 
5.4  pm.  A 1.3  km  sea  level  path  transmittance  spectrum  covering  the  1,8  to  14  pm 
region  is  shown  in  Figure  25.  In  the  latter  spectrum  the  transmittance  was  nor- 
malized to  unity  in  the  10  pm  region,  The  anomalous  spike  shown  at  — 13  pm  was 
due  to  electrical  noise.  Figure  26  shows  a long  slant  path  measurement  (~  12  km) 
from  800  ft  to  3187  ft  altitude.  Again  the  agreement  between  these  measurements 
ana  LOWTRAN  3 predictions  is  good.  The  improvement  of  LOWTRAN  3 over 
LOWTRAN  2 in  the  2 - 3.5  pm  region  can  be  seen  from  Figure  25. 

Figure  27  shows  three  aircraft  measurements  made  by  Camming  et  aF 
(CARDE)  in  the  2.7  pm  region  at  an  altitude  of  13.7  km.  The  measurements  were 
made  over  Florida  using  the  sun  as  a radiation  source. 

In  the  LOWTRAN  3 predictions  shown  in  Figure  27  the  sensitivity  of  the  trans- 
mittance to  atmospheric  conditions  is  indicated.  For  example  the  upper  figure 
shows  the  transmittance  obtained  from  LOWTRAN  3 u8ing  the  tropical  and  1962 
U. S,  Standard  Atmospheres  compared  with  the  CARDE  measurement.  The  second 
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29.  Ashley.  G.W.,  Gastineau,  L„  , and  Blay,  D.  (1973>  Private  Communication. 

30.  Cumming,  C. , Hawkins,  G.  R. , McKinnon,  D.  .1,  G. , Rollins,  J. , and 

Stephenson,  W.  R.  (1965)  Quantitative  Atlas  of  Infrared  Stratospheric 

Transmission  in  the  2.7  Micron  Region,  Canadian  Armament  Research 

an<3  Development  'Establishment,  cAHDE  T.  R.  546/65,  Project  D46-38-01-19. 
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figure  shows  a LOWTRAN  3 comparison  using  the  Tropical  Model  Atmosphere 
with  and  without  aerosol  attenuation.  The  latter  curve  appears  to  agree  better 
with  the  measurement. 


GEN.  DYNAMICS 
LOWTRAN  3 


Figure  23.  Comparison  Petween  l.OWTRAN  3 and  General  Dynamics  Measure 
ments  150  ft  (45  m)  Path  at  Sea  I.evel  (H./)  - 0.35  pr  cm/km' 
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Figure  24.  Comparison  Betv  .n  LOWTRAN  3 and  General  Dynamics  Measurements ; 
10,679  ft  (3.25  kn  1 Path  at  S - Level  (H^J  = 0.35  pr  cm/km) 
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Figure  28  shows  two  atmospheric  transmittance  measurements  made  by 
3 1 

Farmer  et  aP  (EMI)  for  horizontal  paths  at  5.2  km  altitude  in  the  Bolivian 
Andes.  The  measurements  were  made  over  path  lengths  of  1.9  and  3.4  km  and 
covered  the  4 um  spectral  region. 

In  general,  LOWTRAN  3 predictions  compare  fairly  well  with  available  at- 
mospheric. transmittance  measurements.  Further  laboratory  studies  will  be 
conducted  to  determine  more  accurately  the  effect  of  continuum  absorption  by- 
water  vapor  primarily  in  the  10  pm  and  18  - 30  pm  regions.  It  is  hoped  that 
further  refined  atmospheric  transmittance  measurements  over  long  atmospheric 
paths  will  be  made  in  the  near  future,  with  much  attention  paved  to  obtaining  ac- 
curate absolute  transmittance  measurements  in  the  window  regions. 

As  more  measurements  become  available,  it  is  planned  to  update  LOWTRAN. 
A version  of  the  program  which  is  also  capable  of  predicting  atmospheric  "clear 
sky”  and  earth  backgrounds  will  be  published  in  the  near  future.  Two  examples  of 
some  typical  radiance  calculations  using  a modified  LOWTRAN  program  are  shown 
in  Figure  29.  Figure  29a  gives  the  clear  sky  radiance  for  a vertical  path  from  sea 
level  to  space  for  six  model  atmospheres  3nd  Figure  29b  shows  both  the  upward 
spectral  radiance  (as  viewed  from  space)  for  a midlatitude  winter  model  atmos- 
phere and  ihe  downward  radiance  at  sea  level.  The  back  body  spectral  radiance 
curve,  superimposed  on  Figure  29b  corresponds  to  a ground  temperature  of 
272”K  (assuming  a ground  emissivjty  of  unity). 


31.  Farmer.  C.  B. , Berry,  P.  .1. , and  Lloyd,  D.  B.  (1963)  Atmospheric  Trans- 
miesion  Measurements  in  the  3.5  to  5.5  Micron  Band  at  5200  \Tetors 
Altitude,  FTR7  fcleclronics  Limited,  Hayes,  Middlesex,  England,  Ttepo rt 
No'.D'MP  1578. 
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Comparison  Between  IOWTRAN  3 and  E 


mrl  *»'l  3JN»I0*» 


TO  r 


y 


spectral  radiance 

TOR  A 272**  Black  BOOT 


-SPECTRAL  RADIANCE  LOOKING  DP 
E RON  SEA-LEVEL 
MiOLATlTUDE  WINTER  MODEL 


* 

\ 


SPtCTfiAL  RApiANCC  LOOKING 
DOWN  PrtOy  SPACE 
MIDLATiTUOC  W‘NT£R  MOOCL 


II)  1?  IS 

WAvELESGTh  |mtt  ) 


2! 


25 


25 


! 

i 


\ 

I 


Figure  29.  Two  Kxampltrs  of  Some  T ypical  Radiance  (.'alculations  l mil 
Modified  I.OWTRAX  I’rOL’iani 
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Appendix  A 

Listing  of  Program  and  Data 

A listing  of  the  Fortran  pi'-ugi-am  LOWTRAN  3 is  given  in  Table  Ai  together 
■with  the  two  subroutines  POINT  and  ANGL  (see  Appendix  C for  an  explanation  of 
subroutine  ANGL).  The  input  data  for  the  program  is  given  in  Table  A2.  A gen 
eral  flow  chart  for  the  main  program  is  presenter  in  Appendix  B.  and  definitions 
of  the  symbols  used  in  the  computer  codes  are  summarized  in  Appendix  D. 

The  subroutine  POINT  has  a twofold  purpose.  When  the  subroutine  is  called 
for  a giver  altitude  X,  it  is  used  to  determine  the  mean  refractive  index  (1)  in  the 
layer  between  X and  the  level  above,  TX{9',  and  (2)  in  the  layer  between  X and 
the  level  below,  YN,  In  addition,  an  interpolation  scheme  is  used  to  determine  the 
effective  absr  rbtr  amounts  per  km  at  altitude  X for  each  absorber.  When  the 
parameter  IP  is  set  equal  to  z'  > only  the  mean  refractive  index  above  and  below 
altitude  X m determined  from  T. 

The  subroutine  ANGL  is  usea  ■'olely  for  the  purpose  of  calculating  the  initial 
zenith  angle  (0  or  ANGL)  by  an  iteraiive  scheme  taking  into  account  refraction, 
given  (1'  the  initial  and  final  altitudes  of  the  path  (HI  and  H2  respectively)  and  the 
angle  subtended  at  the  earth's  center  (f>  or  BETA!  Dy  the  trajectory,  or  (2)  the 
initial  altitude  and  tangent  height  (HI  and  HMIN  respectively'.  Examples  of  two 
typical  problems  involving  the  use  of  the  subroutine  ANGL  are  given  in  Sections 
C.6  and  6.7  . An  explanation  of  the  iteration  scheme  is  given  In  Appendix  C. 

The  changes  necessary  to  update  LOWTRAN  2 to  LOWTRAN  3 are  indicated 
by  the  symbols  *,  + ,A,B.  C etc.  against  the  card  sequence  numbers  in  Table  AI, 
The  - symbol  indicates  that  the  following  card  (in  LOWTRAN  2)  has  been  removed. 
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Table  Al.  Listing  of  Fortran  Code  LOWTKAN  3 


PROGRAM  LONTRANd  NPUT  , 0UTPUTvTA3E5:INPUT,TAPE6=0UTPUT)  A 

COMMON  Zl  341  (7,34)  . T(7,  34  T ,FH<10, 3 4 J ,MH(7, 34  I ,M,Nl.  , RE, CM, CO, PI  A 

01 ME NS I ON  M0(7,34> ,MZ1 (34) ,HZ2 (5) ,AHAZE( 34) , AHZ2I 20>  A 

DIMENSION  TPI5  7) ,FH<47) ,F3(67I ,HZ(2) . TXI 1CI , VHI10 I ,W(1 j| , E 119)  A 

DIMENSION  Cl(?5fl8),  C2I1S75),  C3(540l,  C4(133),  C5(15),  C8(132»  A 

DIMENSION  VXC45T  ,C7(45I,CTA(45)  A 

F(A)  - FX«M1  8.9266-14. 9595*A- 2 . 4381 2* A »A > » A A 

DATA  H2(1)/5H?3  KM/,KZ(?)75H  5 <M/  A 

!•••••>.  >1  M««  M > .11  •>«>«•  .Ml  A 
PROGRAM  L0NTRAN3  CALCULATES  THE  TRANSMITTANCE  OF  THc.  ATMOSPHERE  A 

FROM  3 SO  OM~l  TO  40001  CM-1  <0.25  TO  28. ST  HI'RONS)  AT  20  CH-1  A 

SPECTRAL  RESOLUTION  ON  A LINEAR  WAVE  NUMBER  SCALE.  A 

REFRACTION  AND  EARTH  CURVATURE  EFFECTS  ARE  INCLUDED.  ATMOSPHERE  A 
IS  LAYERED  IN  ONE  KM.  INTERVALS  9ETWEEN  C AN0  25  KM.,  5 KM.  INTER-  A 
VALS  TO  Cl  KM.,  A TWENTY  KM.  INTERVAL  TO  7C  KM.,  AN0  A THIRTY  KM.  A 
INTERVAL  TO  100  KM.  A 

PROGRAM  ACTIVATED  9Y  SUBMISSION  OE  -OUR  CARD  SEQUENCE  AS  FOLLOWS  A 

A 

CARD  1 MOIKL. IHAZE. ITYPF.LEN, J P, TM, Mi , N2 ,M 3 , ML , RO  FORMA T C 1 C 13 . F 1 0 . 3)  A 

CAPO  2 Hi, M2, ANGLE, RANGE, 9ETA, VIS  FORMA  T(ZFJ0.3)  A 

C A°C  3 VI,  V?,  OV  FORMAT  C ZF  V 0 . 3 i A 

CARO  4 rxY  FORMAT (13)  A 

A 

HC0EL=1,2,3,4»5  OR  6 SELECTS  ONE  OF  THE  FOLLOWING  MODEL  ATMOSPHERE  A 
TROPICAL, MIDLATITUOE  SUHME « . HI 0L A T I T UDE  MI NT E? , SU BARCI I C SUMMER,  A 
SUBARCTIC  WINTER, OR  THE  1962  U.5.  STANDARD  R-SPECTIVELY  A 

MOOFL-P  FOP  H0RI7.  payh  MHEN  METEOROL.  DATA  USED! INSTE  AD  OF  CARD  2 A 
RFAO  «l,PIMB»,TfOfr.  C!tOEW  PT.TCMP'OEC  C?,*tr\CL  riU*TIDITY,H<0  ug.  N j t 4 
IGM.N-7),03  O^NSITYJGh.m-JI , V 1 S ( KN) , RANGf I K M)  Ml TH  FORMAT  429,  A 

MODE  I.  = 7 WHEN  HEM  MOOFL  ATMOSPHERE  IE.G.  RADIOSONDE  DA  T A I USED.  A 

OATA  CARDS  A Rc  READ  In  BETWEEN  CAROS  1 AND  2,  AND  SHOULD  CONTAIN!  A 

ALTITUOEtKM. ), PRESSURE, TFMP, OEW  PT. TEMP, REl.  HUMIDITY, M?0  DENSITY,  A 
OZ  DENSITY, AEROSOL  NO,  DENS  I TY  I CM -31  ACCORDING  TC  FORMAT  429,  A 

NOTE  THAT  EITHER  OEM  PT.  TEMP. OR  REL.  HUMIDITY  CAN  9E  USED.  A 

A 

Ml, M2, M3,  APE  USED  TO  CHANGE  TEW»,M20,  AND  03  ALTITUPE  PROFILES.  A 

A 

If  IhAZE=C-  no  ATROSOL  SCATTERING  IS  COMPUTED  A 

IHA7F  =1  IF  «Tt(05Ci  ATTENUATION  REOJIRED  (THIS  IS  USED  IN  A 

CONJUNCTION  MI TH  VISUAL  RANGE  I SEE  CARO  2»)  A 

1 H A *E  - 1 OR  2 ALSO  GTVF  AEROSOL  A TT  ENUA  T ION  c OR  23KH  AND  5KM  VIS.  A 
HAZE  MODELS  R:SPLCriVrLY  IF  VIS  -0  ON  CARD  2 A 

A 

ITYPf  =1,2  OP  3 INDICATES  THF  TYPE  OF  ATMOSPHERIC  PATH  A 

TTYDr--3, VERTICAL  OP  SLANT  PATH  TO  SPACE  A 

ITVPE-2,VFRTICAL  OR  SLANT  PATH  9ETWEEN  TWO  ALTITUDES  A 

ITY»E  = 1,  C0PRrSP0N0S  TO  A H ORI  70N  T A 1.  (CONSTANT  PRESSURE)  PATH  A 

A 

*■ 1-QBEFRVFR  ALTTTUOE  ( KM i A 

H2- SOURCE  AlTITUUE  {KMT  A 

ANGLF=  7FN I TH  ANGLc  AT  HI  (DEGREES:  A 

RANGE=PATH  LENGTH  (KMT  A 

RET A = E APTH  CENTPE  ANGLE  A 

VIS  = VISUAL  RANGF  AT  SEA  LEVEL  (KM)  A 

(IF  I T YPF  = l «rAO  HI  ANO  RANGE!Ir  ITrPE=3  RFAO  HI  AND  ANGLE.  A 

IF  I TYPE- 2 RF6D  HI  ANU  TWO  OTHER  PARAMETERS  E.G.  H2  AND  ANCLE)  A 

A 

VI - INI T I AL  FPtQUENCY  ( WAVrNUM8ER  CM-1  ) INTEGER  VALUE  A 

V2=FINAL  FREQUENCY ( M3 VENUMBEP  CM-1  ) INTEGER  VALUF  A 


1 

2* 

?• 

4 

5A 

53 

5C 

6 

7 

8 
9 

:e 

n 

12 

13 

14 

15 
lb 
1? 
18» 
19 


21 

22 

23 

24 

25 
26* 
2/» 
26* 
29A 
293 
?9C 
29D 
29E 
29F 

29G 

30 

31 
3? 

33 

34 

35 

36 

37 
3 = 
39 
'.C 

41 

42 
-3 

44 

45 

46 

47 
-8 
49 
5 C 
51 


52 
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Table  Al.  Listing  of  Fortran  Code  LOWTRAN  3 (Cont) 


0V=  FREQUENCY  INTERVALS  AT  HH1C9  TRANSMITTANCE  IS  PRINTED 
NOTE l DV  MUST  8E  * NULHPLE  OF  5 CN-1 

JXY  = 0 TO  END  DATA  , = 1 FOR  NEW  Y1,¥2,0V  ONLY  , s 2 TO  CONTINUE  DATA 
IXT-r-l  FOR  NFM  CtiRO  2 ONLY,  = 4 FOR  NEK  CARO  1 ONLY. 


I ATM,  ML 

tHZiCII .1=1, NL> 
(H72tXI ,1=1,5) 


Z(I)*(P<K«XI»r<K»l»,ICH(K,I),UO<X, I)  ,K=K1,K2I 
(VXII».C7lXl.C7A(X)  .1*1,44) 

(TRIII  ,FH(n,EotI),l  = l,67> 
fCtfII»I=l»258QI 
(C2( II ,1=1,15751 
(C3(I>»X=1.540> 

(C4(I> .1=1.133) 

<c5(it«x=i«i5> 

( C8( I) ,1=1, 192) 


IXV -0 

PESO  15,4081 
PE«0  (5,401) 

READ  (5,491) 

00  1 J = 1,  3 
K2=2*J 
Kl*«2-i 
00  1 1=1, NL 
READ  (5,452) 

»E«0  (5,431) 

REA0  (5,403) 

READ  (5,4041 
READ  (5,404) 

READ  (5,494) 

READ  (5,405) 

READ  (5,404) 

RE *0  (5,495) 

PI=?.0*ASIN(1. 0) 

CA=PI/180. 

IP=9 

• CONTINUE 

RE =6371. 23 
IFIND=# 

; JR  WE  0 SUFRESS  PftiN) 

READ  480, MODEL, IHAZE, I TYPE , LFN.JP,tK»Hl,N2»N3»NL,PO 
PP I NT 40  0, MODEL ,TNAZE,ITYPE,L£N,JPrXM,Hi,H2»K3»NL,  RO 
: PRINT  424,  M03EL,IHAZF,ITYPE,LCM 

260  M=fiODEL 

If  (M.EQ.l)  RE  =637 6*35 
IF  (H.EQ.4)  RE=6356.91 
IF  (M.FQ.S)  RE*6J56.91 
IE(I)(Y,GT.3)  GO  TO  6 
IF(RO,NE.0.*>  RE=»0 
IF(t4.E0.7.  AND.  IM.«C.0)GO  TO  4 
IF  (NOOH..ea.O)  GO  TO  4 
380  READ  406,  Ml , H2» ANGLE, RANGE, BE TA , VI S 

PRINT  425,  HI, H2, ANGLE , RANGE  * BE FI • ¥1 S 
»Mi 

IF  (ITYPE.EQ.3)  GO  TO  560 
IE  (lTYPE.E0.il  GO  TO  8 
X2=RF*M2 

IF  (RANGE. ea.0.)  GO  TO  5 

PRINT  428,  Ml, H2, ANGLE, RANGE, BETA, YIS 

IE  (M2. EQ.O. AND. ANGLE. NE.C)  GO  TO  S 

ANGIE=AC0S(Q.3*((M2-Hl)*(l.+X2/Xl>/IANGE-R*NGE/Xi>l/C» 

GO  TO  7 

3 X2  = SQRT ( ( Xi/RANGE»RANGE/Xt  *2.G*I0S( ANGLE*CA) ) • XI 'RANGE ) 

M2=X2- RE 

GO  TO  7 

4 CONTINUE 
IE(NL.LE.0INL=1 
DO  540  K = l,  M*. 

AMAZE (K)=C.O 

IF(M.EO.9>REA0  429,  HI, PI  7,1 1 ,TN|S,0P,RH,WM(7,KI  ,NO  ( 7,  K)  , Yl  S .RANGE 
IE (M.EQ.O) PRINT  438,Hl,Pir,l»,TNP,DJ,«H,NH(7,K»,M0(7,(O,YIS,RANGE 


A 

53 

A 

54 

A 

55 

A 

56 

A 

57A 

A 

578 

A 

5 7C 

A 

58 

A 

59 

A 

60 

A 

61 

4 

62 

A 

63 

A 

64 

A 

65 

A 

66* 

A 

67 

A 

6fi 

A 

69 

A 

?n 

A 

71 

A 

72 

A 

73 

d 

74* 

A 

75 

A 

76 

A 

77 

A 

78 

A 

79 

A 

79* 

A 

80 

A 

cl 

A 

81 

A 

9 2 

A 

83 

A 

64 

A 

65A 

A 

858 

A 

85C 

A 

650 

A 

86 

A 

87* 

A 

88 

A 

c 9 

A 

9C* 

A 

91 

A 

92 

A 

93 

A 

9- 

A 

5* 

A 

36 

A 

97 

A 

96 

A 

99 

A 

100 

A 

101* 

A 

102* 

A 

f d JA 

A 

1038 

A 

103C 

A 

1030 
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$00 


5 no 


C 

C 

560 


580 


IF(h.gt„o » Rfarj  429, z<  k)  ,Pt7  , 0 ,th?,op,rh,hh  <r,  ki  . wo  1 7 ,ki  , ahaz?  < ki  4 U3£ 

JsIFITIZtiO  M.0E-6)*1.  4 103F 

IMI.ECLOJZIKMHl  4 i:jg 

IFf?(KI  .GF.25.0t  J=IZ<tCI-25. 01/5. 0*26.  * 1QJH 

IMZ  OO  .GE.50. 0)  J~«Z(KI-50.C>  /2O.0+31.  4 1 : 31 

rrtroo  .gf.to. 0 j=(Z(<i~?D,ct/io.a»32.  * 103J 

IfU.GT.33l  J*33  4 1 3K 

F4C  = 7(tCI-FL0AT  « J- 1 » A 133L 

IFU.LT.26t  GO  TO  500  » 103M 

FAC=(7(K)-5.0*FLOAT(J-26)-25.l/5.  » 133N 

IFU.GE.J1I  FAC=<Z(KI -SCI, 31/20.  4 1 = 30 

IFU.GE.3~5  F4C=  IZtKl -TO  /30  . A 103P 

IF(FAC.GT.l.O)  FAC=1.0  4 113Q 

l=j»l  4 X 0 35t 

T(7,K)=TMP*?73.15  * 103S 

IF<Hl,GT.O»T<T,IO  = TtMi,J»»irCHi,l.»/T«Hl»JH»*F4C  A 103T 

TT=2?3.15/T 17, K>  A 103U 

IF<«tH.LE.  0.01  TT  = ??'  3.  15/I?73.15*0P)  * 103V 

TFINHf/.Kt  ,L£. 0.05  WH ( 7.  K)  = F ( T T ) 4 1C3H 

IF<H2.GT.0l  NH(  7.M  --HMM2,  J)  * < WHt N 2 , L > / WH ( H2 , J)  ) »»F»C  4 10JX 

It- (RH.GT.0,0)  HHI7,K)=0.Ol*RH*tH(7,<)  * 1 3 3T 

IF(M?.GT.0)WOt7,K»=HOIH3cJ|»l«O«13,L»/MO(H3,J)ls»FAC  4 1j3Z 

IF(Z(K» .GE.5.0 IGO  TO  520  4 1044 

IF(AMA7E(K) .£9.0,0) AH72<K>=HZ2 (J)  *IHZ2(L)/H72< Jt t **FAC  » 10AB 

IFtAHAZEtK) .EO.O.O) AHAZf CK)  = H?1(  J)  •(  HZltLI/HZl (J»> **F4C  4 1J4C 

IFtHOOEL.Ea.05  GO  TO  8 4 1040 

IFtK.FQ.DPRlNT  *.41  4 i04E 

PRINT  429. Z(K)  ,P(7,K»  , TMP, OP ,R H, WH J 7 , * I , "0 ! 7 . K } , AM A ZE J KS  « 1U4F 

CONTINUE  4 1S4G 

TN=S  4 104H 

NL=*tl  4 1041 

K1=0  A 1 04 J 

H2  = 3 4 1S4K 

M’sO  4 134L 

NOTE  THAT  Ztll  HAT  NOT  CORRESPOND  TO  THE  VALUES  GIVEN  FOR  STANDARD  A 104H 

MODEL  ATHOSPH-RFS  A 104N 

GO  TO  300  4 1140 

IF  (R4NGE.GT .0.0)  GO  TO  5=0  A 134P 

IF  ( H£ . GT  . 0 » 0 o AND  .H2.LT.Hil  Ifi*(3  = i t 1043 

GO  TO  8 A 1 C4R 

ITYPFi?  4 104S 

eFTA=ACOStC.5,,tRANGE*RANGE/t)Cl**2)-<2/Xl-Xl/K2  5 5 /CA  £ 104T 

IF  <OFTfi.E5.c.5  GO  TO  6 4 IS  5 

IFINO-1  4 106 

BETsC/i’BFTA  4 117 

X2=RE»H2  4 10? 

ANGLF=r4T4NtX2*SIN«BET)/tX2»C0S(3Er)-Xll)/C4  4 109 

RANGr=X?*SINtBEr?/SINt4NGLE*C4)  4 110 

BF’irBFTfc  A 111 

GO  TO  8 4 112 

PANGE=(X2/Xlt**2- (SINSANGLE*CA)I **2  A 113 

IF  f RANGE  . GE. 3 . 0 5 RANGE-X1 * C SORT t RANGE  I - 4B5 fCOSt ANGL E*C4) > 5 4 114 

IF  TANGLE. NE.3..0R„ANC'LE,NE.18C«t  9E I = 4S IN (RANGE* S IN (4 NGL E *CA5 /X? 5 A 116 

IF  (4NGLE.LT .3 .5  ANGLE^AMGLETPI  4 116 

IF  (RANGE. I.T. 0.0)  R4NGE  = -»4NGF  4 117 

BET  =BET/CA  4 iU* 

PRINT  428,  HI, H2, ANGLE, RANGE, BET,  VI3  4 119 

CONTINUE  A 120A 

SUH4-0.  A 12CB 

IFtIXY.LE.2)  READ  406,V1,V2,OV  4 121* 

IF(IXY.LE.2)PRIMT  406.V1.V2.0V  A 122* 
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IF  (TTVPE.EQ.l)  PRINT  40  7.  Hit  RANGE  A 123 

IF  «ITVP€.EQ.2>  PRINT  408.  Ml, M2,  ANGLE  * 12** 

IF  ( ITVPE.E0.3  > PRINT  409,  Hi, ANGLE  * 125 

IF  (MODEL. £0.9?  N=7  A 126A 

IF  tM.EQ.7)  PRINT  417,  VIS  A 1263 

IFIVIS.LT,  2.  0.  ANO.VIS.GT.O.OI  PRINT  442  A 126  C 

IF  (M.EQ.ll  PRINT  410 , N A 127 

IF  (H.EQ.2?  PRINT  Nil , N A 128 

IF  (H.EQ.  3)  PRINT  <*12,  H a 129 

IF  fN.Ea.9l  PRINT  Nil,  N * 130 

IF  (M.EQ.5)  PRINT  <*15,  M A 131 

IF  (H.EQ.6)  PRINT  414,  N A 1 32 

IF  (IHAZE.EQ.fi  , ) PRINT  <*26  A 133 

IF  (N.NL.7.AND.IHAZE.GT.0)  PRINT  416,  INAZE.  I- Z ( IHAZE)  A 134* 

SVN=10000./V1  A 135 

ALAM*ififiO0./V?  A 136 

PRINT  416,  VI, V2»OV,ALAN,AVW  A 137 

AVM*fl.5E~4*(Vl*V2»  A 138 

AVM*AVW«AVM  A 139 

C0*77.46*.459*AVM  A 140 

CK*43. 4*7-0. 34?3*AYM  A 141 

IF  i IrINO.EQ.l > GO  TO  15  A 142 

9 IF  (IFINO.EQ.il  CALL  ANGL  (HI, M2, ANGLE. BETA, LF N, HU  A 143* 

IFIND=fi  A 144 

IF  tJP.EQ.fi>  ’PINT  427  A 146* 

IF  (ITVPE.EO.il  GO  TO  15  A 147 

00  11  K*l,ifi  A 148 

VHIKIsO.fi  A 149 

11  CONTINUE  A JSC 

OETAsC.fi  A 151* 

SRsfi.S  A 153 

IP*fi  A 154- 

0****  NOH  OFFINE  CONSTANT  PRESSURE  PATH  QUANTITES  £H(l-a>  A 156 

V*CA*4NGLE  A 157 

SPHIsSIN(r)  A 158 

Ri*(RE*Hll*SPHI  A 159 

IF  (Hl.GT.ZtNLI I GO  TO  13  4 160 

GO  TO  15  A 161 

13  X*(RE«-Z(NL»»/(RE*H1>  A 162 

IF  (SPHI.GT.XI  GO  TO  14  A 163 

M1*Z(NL>  A 164 

J1»NL  A 165 

5PHI*SPHI/X  A 166 

AWGLE*18#.0-ASIN(SPHII/C:a  A 167 

R|* <RF*H1!*SPHI  A 166 

CO  TO  15  A 170 

14  HMINsRl-RE  A 171 

PRINT  433,  KMIN  A 172 

GO  TC  95  A 173 

15  00  17  1*1, NL  A 174 

PS*P4N» 117 1C13.0  A 175 

TS*273.15/T(«,I>  A 1 76A 

IFtNi.GT.fi. AMD. N.LT.7»T5  = Z7  3.157T(!U,I>  A 1768 

X*PS*TS  A 177 

Pf*PS*SCRT(TS>  A 178 

0*fi.l*«MtN,It  A 179 

IF(H2.GT.C.AN0.H.LT.7>  0*0 . 1*NH(N2« I > A 160* 

EH(l,n*0*PT*»fi.9  A 181* 

EHC2,I)*X*PT**8.7S  A 182* 

EM(4,n««.8*PT*X  A 183 

PPM*4,S6C-5*0*273, 157T5  4 164* 
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EMC5. Tt=tPPN+0 .UOSMPS-PPWM  *0 
EM(6,T»=X 

haze=kzi(ii 

iFtH.ea.?!  haze=ahaze»ij 
jrizm.GE.s.o*  co  to  iso 

IFITHAZE.EQ.2J  HAZE=HZ2fI* 

IFIIHAZE.EQ.2.  AND.H.EQ.7>NAZE  = AHZ2m 
IFIVIS.LE.C.OI  CO  TO  ISO 

HAZE*  6.389*(CHZ2(I»-HZ1 (III /VIS*HZ1  1 1 1 1 5 . B-HZ 2 ( 1 1 /23 . 0 J 
IF  (H.NE. 7)  GO  TO  150 

MAZE=6.389*4  «AHZ2(II-AHAZr t I J I /VI S + * HA  ZE  C I > / 5.  Q -AHZ21 1 1 /2  3.  3 1 
158  IF(NA7E.LT.0.0 I HAZEsC.C 
EH(7f  11=3. 5336 £-6*HAZE 
FH(3,TI  *^6.666 7*U0(N, I) 

IF(N3.GT.C.AN0.«.Lr.7l  EH  (8,  1 1 =1*6 . 65  /*MO I N3 , 1 J 
EH(7,T}rfHt8,r}*PT**0  .<* 

EHtq.ri =1.0 

EH(lO„II=l.OE-6*(CO*X*10H.  0/273.  15-PPN*CW) 

IF  (I.EQ.NLI  CO  TO  16 
IF  (MODEL.Ea.O.ANO.I.GE.ll  GO  TO  26 
T2=T(H, I*i» 

W2=NH(M,I*1I 

IF(Hl.GT.0»T2=TtHl,IMI 

IF0H2.Gr.P9W2=HHtK2.I*lJ 

PPH=«*.56E-6*M2*T2 

rH(9.n*d.5*tFH<l0,I»  ♦1.0F-6*tCO»PM.I*-l)/T2-PPM»CJt)  J 

16  IF  (I.EQ.NLI  EH(9,II=C. 

IF  (Hl.GE.zmi  Ji=T 

IF  (IFINO.EQ.O -OP.  JP.EQ.O!  PRINT  •tiki  1 * 2(1  J , i £Kt  K.  I i *K=i  r iui 

17  CONTINUE 


A 1 05 
A 186 
A 1&7 
A 188* 
A lo9* 
A 19:  A 
A 19G6 
A 190C 
A 1500 

a i<5: e 
A 1 OOF 
A 190  G 
A 191 
A 192 
A 19  34 
A 1936 
A 193C 
A 1930 
A 196 
A 1 95A 
A 1958 
A 195G 
A 1960 
A 1 9SE 
A 196* 
A 197* 
A 19‘* 
A 199* 
A 200* 
A 2:1 
A 202 


IP=1 
TK  = 0 
Xl=Hl 

CALI.  POINT  (Hl.VN.N.NPl.Tx,  IP! 

J1  = N 

TX1=TX(9I 
00  10  K=1,S 
18  EtKJsTXKJ 

IF  (ITYPF.EO.l)  GO  TO  26 
IF  !ITVPE.EQ.*I  H2=Z(NLi 
IF  (ANGLE. GT. 90. 01  GO  TO  28 
15  IF  ( ANGLE • CT .90. s.ANO . NP1 . G V . 0 t J1=j1*1 
J2  = Nl 

IF  (ITVPE.Ea.3J  GO  VO  20 
CALL  POINT  (H2,YN,NtNP,TX,IPJ 
J2  = N 

IF  (NP.CT.9)  J2= J2-1 

20  DO  21  K=1 . 8 
FH ( X. J1 1 - E ( XI 

IF  (ITVPE.EQ.3I  GO  TO  21 
EH(X,  J2U»=TX(K) 

21  CONTINUE 

IF  (J1.EQ.92I  TXl  = rXKYN-EH(9,  Jll 
C***»  NOW  DEFINE  VERTICAL  PATH  QUANTITIES  VM(1-8I 
IF  (JP.EQ.OI  PRINT  <*2  0 
00  ?5  I=J1,J2 
X1=Z(II 
X2*Z(I*1I 

IF  (I.EQ.J1J  Xi=Hl 


A 2 06 
A 2C5 
A 206 
A 207 
A 208 
A 2:9 
A 210 
A 211 
* 212 
A 213 
i 21  6 
A 215 
A 216 
A 217 
A 218 
A 219 
A 220 
A 221 
A 222 
A 223 
A 226 
A 225 
A 226 
A 22 7 
A 228* 
A 229 
A 230 
A 231 
A 232 
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I r (I.EQ.J2I  X2»M2  A 233 

0Z*X2-X1  A 236 

TF  (I.EQ.Nl)  B2*zm-z«-l>  A 235 

OS*OZ  A 236 

C****»  UPWARD  TRAJECT ORV  A 23/ 

RX*fffE*Xl)/<RE*X2l  A 23# 

THETAxASINCSPMIl/CA  A 239 

PHI»ASIN(SPHI*RX)/GA  A 260 

«ET*THFTA-1HI  A 261 

SALP*RX*SPHI  A 262 

IF  CSPHI.GT.l.E-lOl  0S*IRE*X2> *SINCSET*CA>/SPHI  A 263 

BETA*RETA*8ET  A 266 

PSI»BETA*PHI-*NGLE  A 265 

PMIslflO.-PHI  A 266 

SP*SR*OS  A 26/ 

00  26  Ml,#  A 268 

EV*0S*FM(K,f»  A 269 

IF  (I.EQ.NL)  GO  TO  22  A 250 

IF  IEHttC,I>.EO.O.O.OR.EHIK,I*l>.£Q.a.O>  GO  TO  23  A 251 

IF  <EHCK,I>.£Q.£H<K,I*in  GO  TO  26  A 252 

[ EV*0SMEHCK,II  -EK(K(r+llt/At.eG<£H(X'I>/EHtK'X»il>  A 253 

GO  TO  26  A 256 

22  IF  <£HCK,I).EQ.0.0>  GO  TO  23  A 255 

IF  CEHCK.I-U. £0.0.0)  GO  to  23  A 256 

IF  <EH<K,I).E3.EH(K,I-in  GO  TO  26  A 257 

EV«6V/AL0GSEMtK,I-i»/EHIK,in  A 258 

« GO  TO  26  A 259 

? 23  e**a.  a 26o 

26  VHIK)xVH(K>6EV  A 261 

IF  CJP.EQ.OI  PRINT  635,  ITXi|(V1(L)tL=l,6),PSI ,PHI , BET A, THETA, SR  A 262* 

IF  CI.GE.NL)  GO  TO  25  A 2G3 

IF  CI*i.EQ.J2>  EHC9,I*1>=TN  A 266 

IF  CI.FCt.Jll  EM (9, 1 > * TX1  A 265 

1 RN*EHI9,I*1»/EHC9,II  A 266 

SPHI*SPHI*RX/«N  A 267 

i IF  CSALP.GE.KNf  5PKI=SALP  A 268 

< 25  CONTINUE  A 269 

i GO  TO  67  A 270 

I HORIZONTAL  PATH  A 271 

> 26  00  27  K«l,8  A 272* 

miO*R*NGE*EHIK,l>  A 273* 

IF  ! MODEL  . GT  .fl  ) HCK3 -RANGE*TXCK>  A 276* 

27  CONTINUE  A 275 

j GO  TO  69  A 276 

26  CONTINUE  A 277 

C'"5  OOwNWSRu  f K* JfcGf ORT  A 276 

KZ~9  A 279 

IE  INPl.EO.il  J1*J1-1  A 280 

J2«J1*1  A 261 

VN1«YN  A 282 

J*J1*1  A 233 

IF  CH2.GT.2C Jl*l) .0R.M1.E3.H2>  GO  TO  30  A 266 

IF  CNP1,EG>.1.AND.H2.GE.Z<J1*1>)  GO  TO  30  A 285 

CALL  POINT  (H2,VN,N,NP2,TX, I!»>  A 266 

00  29  K*1 » 6 A 287 

29  NIK)  *TX  CK>  A 268 

TX2«TXC9>  A 269 

VN2*TN  A 290 

IF  (H2.LT.H1)  H*H2  A 291 

J2=N  A 292 

IF  (J1.EQ.J2)  TK2»TX1*TN2-EHC9,N!  A 293 
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30 


31 

32 


33 


34 


O*** 

35 


IF  (H2.GT.H1I  TX1=TX2 

IF  (Jl.EQ. J2.»N0.H2.Lr, Hl»  YNl*TX2 

AO  - <9EfH1 1 *SPMl»rNl 

IF  (H2.GE.H1)  YN2=YN1 

00  31  1=1, J1 

MMIN-sA0/EH(9fl}-(?E 

IF  ( I .EQ.  J1  ) H*1lN  = 40/VNl-IF 

JMIN=T 

IF  (HNIN.lE.ZI  mil  GO  TO  5? 

CONTINUE 

x=mmin 

IF  (HNIN.LE.O)  GO  TO  34 
CALI  POINT  IX, YN.N.NP.TX.IPI 
JHIN=N 
TX3=TX (91 


IF  I J2 • EQ. N. 0? . Ji . EQ. Nf  TX 3 = YN2» TX ( ) I -£M  ( 9 
IF  CJ1.EQ.N .AN0.M2.GE .HI*  GO  TO  33 
HHlN=40/TX3-8E 


,N> 


IF  IANS  CX-HHIN ).GT. 0.00011  GO  T3  3? 
IF  (J1.EQ.N.AN0.H2.GF.H1*  YN1=TXJ 
IF  < J2.  EQ.N.  AH0.  Jl.  NE  . J2I  YN2=TX3 
tF  (H2.GE.H1)  TX2=TXJ 
IF  (H2.GE.H1I  J2=n 

IF  (H2.GE.Ht.08.H2.LT .mhIN)  H=H9IN 
PRINT  436,  MNTN 

IF  (H2.LT.HMINI  PRINT  440,  HHIN 

GO  TO  35 

PPINT  436,  HMIN 

IF  IH2.LT, HI)  GO  TO  35 


TF  *ITVPF.EQ.3.0R.N2.GE.HH 

ITYPEtJ 


PSINT  43? 


TXJsfH ( 4, l* 
JBI N=0 


J?  = l 


H2=0.a 


H = 0.0 


NOW  DEFINE  VERTICAL  PATH  QUANTITIES  VH(l-P) 

IF  (JP.EQ.0)  ORINT  420 
00  40  I-l.NL 
J = J-1 

RFF  = FH  O,  J) 

IF  (I.EQ.l)  PEFiYNI 

IF  It. FQ.l. AND. H2.eQ.il  R=F=VN2 

tC  t I C f*  I *%  A H A «#•»  r A A.  — 

<v«v  •n.iw«Ac«r.jt  vi  iir  si  ^ 

if  (I.NE.n  Xl=2iJ»l* 

X2=7 ( J ) 

IF  (J.EQ.  J2.ANO.K2.EQ.O)  X2  = H 
IF  (J. F0.JNlH.AN0.K2.EQ.il  X2*MNIN 
HH=r|REFXl)*SP“fI-PE 

IF  (HH.GT.2IJ) .AN0.HN.GT.X2I  X2  = HN 
«5X=  IRE»X1)  / (R=*X2) 

05=X1-X2 

ALP=90.0 

TMETsASINISPHII /CA 
S«LP=RX*SPHI 

IF  (A<tS(X?-HM)  .GT.1.0E-5)  ALPiAjlNISALPIXC# 
BET-ALP-THET 

IF  (5PHI.GT.1. OE-10)  OS=(RE»X2l»SIN(BET*CA*/SpHl 

TMETA=180.0-THET 

bfta=retafbet 

P5I=BET*-ALP-ANGLEf1S? .0 


A 294 
A 295 
A 2 j6 
A 297 
A 29# 
A 299 
A 300 
A JJ1 
A 302 
A 323 
A 3. 4 
A 305 
4 3.6 
A 3:7 
A 3 J 8 
A 3S9 
4 31C 
A 3 1 1 
A 312 
A 313 
A 314 
A 315 
A 316 
A 317 

* 31# 
4 319 
A 32. 

A 321 

A X81 

-»  SA  «-  <U 

4 323 
A 324 
A 325 
4 326 
A 327 
A 32# 

A 329 
4 33f 
4 331* 
A 332 
4 333 
4 334 
A 335 
4 336 
A 337 

* 33# 

A 339 
A 340 
A 341 
A 342 
A 343 
A 344 
A 345 
A 346 
A 347 
A 348 
A 349 
A 350 
A 351 
A 352 
A 353 
A 354 
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Table  A 1.  Uatin*  of  Fortran  Code  LOW  TRAN  3 (Cont) 


SR=SR*DS 

A 

355 

00  39  X»i,« 

A 

356 

AJ=EH4K,J1 

A 

357 

8JsEH4K,J*l> 

A 

356 

IF  IJ.EQ.J1)  8 J=E  4 Kt 

A 

359 

IF  IJ.  Ea.J2.ANO.H2.LT. HI. ANO.H2.GT.O. 01  AJ*KIKJ 

A 

36r 

IF  4J.EQ. JNIH.6N0.H2.GE.H1)  AJ 

= TX  <<t 

A 

361 

IF  4J.EQ.JHIN.AHD.*BS<H2-HHt .LT.1.0E-51  AJ-TX(K> 

362 

IF  4K2.EQ.0t  00  TO  36 

A 

363 

IF  IJ.FQ.J2>  9J=KIKI 

A 

366 

IF  4J.EQ.JMIN)  AJ  = TX  C K) 

A 

365 

36 

IF  lAj.EQ.O.O.OR.BJ.EQ.C  .■■)  GO 

TO  39 

A 

366 

IF  4AJ.EQ.BJ)  GO  TO  37 

A 

367 

FV*0S®  4 AJ"BJI/ALOG ( A J/BJ > 

A 

368 

GO  TO  39 

A 

363 

37 

F**  OS»AJ 

A 

3 7 0 

GO  TO  39 

A 

371 

36 

EtfcO.Q 

A 

372 

39 

VM4K)=WH(K»6EV 

A 

373 

IF  4 JP. EQu  01  PRINT  639,  J, X 1 , « VN4LI , L = 1, 8» ,PSI , ALP.BFT A, T NETA , SR 

376* 

IF  4J.EQ. J2.AN0.H2.GE.M1)  GO  TO  69 

A 

375 

IF  4J.FQ.JNIN.ANO.K2.FQ.lt  GO 

TQ  63 

A 

376 

IF  (J.NE.l!  RN=REF/EH  49. J*1 1 

A 

377 

IF  (J.EQ.J2*1I  RN-RFF/T12 

A 

378 

IF  f J.EO. J2. AN0.K2.EQ.0)  «N=REF/tN2 

A 

379 

IF  <J.FQ.4JNIS*lt. AnO.K2.EQ.lt 

*N=IEF/TX3 

A 

383 

IF  I5ALP.GE.RNI  RNsl.C 

A 

331 

5PNlsSALR*RW 

A 

382 

IF  4J.EQ.J2.ANO.K2.EQ.Ol  GO  TO 

61 

A 

303 

60 

CONTINUE 

A 

386 

61 

IF  IHNIN.LE.et  GO  TO  67 

A 

385 

IF  ILEN.CQ.8l  “RINT  67# 

A 

386 

IF  4LFN.EQ. 0 1 GO  TO  67 

A 

307 

IF  4LEN.E0.lt  PRINT  679 

A 

388 

K2=l 

A 

389 

XlsX2 

A 

3-70 

IE  4 AGS  4X1— HN I Nt.LE.O.tClt  GO 

T3  67 

A 

391 

HiHNIN 

A 

392 

JXJ2*1 

A 

393 

IF  INP2.EQ.lt  J=J-1 

A 

396 

b*?ifta 

A 

395 

RN*i80.0-*iSIN4  SPHII/C* 

A 

396 

TS-SR 

A 

397 

PS*PSI 

A 

3?  S 

nn  *a •*  ■ 

A 

5*9 

62 

E«KI*VH4KI 

A 

600 

GO  TO  35 

A 

631 

63 

BET  A*  2. *9F  TA-9 

A 

632 

P5I*2. *PST-PS 

A 

603 

$R=2.*SR-TS 

A 

606 

C 

LONG  RATH  TAKEN 

A 

605 

PHIsPM 

A 

606 

DO  66  K»l,8 

A 

607 

66 

VHIKt  «2.*VN4K) -t.KI 

A 

608 

GO  TO  67 

A 

609 

69 

00  66  K*l,ie 

A 

610 

66 

VHIKt  *2.0®VMiO 

A 

611 

BETA*2 . 8*BCT  A 

A 

612 

SR *2. 8* SR 

A 

613 

IF  IH2.EQ.Mlt  GO  TO  67 

A 

616 

RN*TXl/rNl 

A 

615 
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Table  A 1.  Listing  of  Fortran  Code  LOWTRAN  3 (Cont) 


SPMI«SIN( ANGLC*CA) 

416 

IF  (SPHI.LT.RX)  SPNIrSPMI/RN 

*1? 

GO  TO  19 

418 

4/ 

continue 

4 

419 

PRINT  496,  HN 

4 l 9« 

00  4*  K=l,1.0 

422 

N(KliVHIK) 

421 

ad 

CONTINUE 

422 

49 

MPITE  <6,4191 

423 

WRITE  (6,421)  (4(1), 1=1, 9) 

<*  2 4 

1 = 1 

425 

L = 1 

426 

TVl=Vl/5. 0 

427 

142=42/5. *.99 

428 

IV1  = 5*IVJ. 

4 2'* 

IV2=5*TV2 

43C, 

IF  (IVl.LT.35fl)  141=350 

431 

IF  (IV2.GT.50flaO)  142=5001 0 

432 

IF  (0V.LT.5.)  nv- 5. 

433 

IOV=DV 

4 34 

IV=I41-I0V 

4 35 

ICOUNT=0 

436 

c**»* 

BEGINING  OF  T WANSNT  T T ANCr  CALC'JLATI  ONS 

4 37 

SI 

IV*t4-H0V 

4 38A 

IF(JP.NE.O)  GO  TO  52 

4 38B 

IF  (ICOUNT.FO. 0)  GO  TO  51 

s,3Q 

IF  t I COUNT • FQ • 50 ) GO  TO  5i 

44  C 

GO  TO  52 

441 

SI 

ICOUNT=0 

442 

PRINT  4 22 

I'yJ 

5? 

00  63  K-1.10 

444* 

TX«K)  = (J.0 

445 

IF  (K.LT.41  TS  <K) =1.0 

446 

S3 

CONTINUE 

447 

ICOUMY  = ICOUNl 4 1 

443 

5UN=(S.,0 

449 

V=I»r 

452 

1=  < l'V-  350  1/5  41 

451 

IF  (IV. IT. 1403)  GO  TO  61 

<•52 

IF  iIU.LT.2749)  GO  TO  68 

45  3 

MOUTCIRAR  SCATTERING 

•4  5 

C6=9.«i«7£-20*I  4**4,01171 

455 

TX(6)*C67M(6> 

456 

SUM=  SUM*T  4(6) 

u 4 7 

Ir  ( IW.  t ' , si 2 80  ) GO  TO  72 

453 

IF  (IV. LT.  13000)  GO  VO  6*3 

4^ 

£•*••«»*•  Uif  0 7 ONC 

4o? 

IF  ( IV. LE » 2 S4P 0)  GO  TO  54 

4bl 

IF  (14. GF. 27590)  GO  TO  55 

462 

Gu  TO  87 

4o3 

54 

XK=  260 , Q 

464 

Xl*« V-lSfcOO.O)  /XX*1.0 

465 

L 1 * 1 

<♦66 

L2  = 53 

467 

GO  TO  56 

A 

4t  8 

55 

XX=5li. 0 

469 

XI=(4-2 7500.0) /XX757.0 

470 

Ll*57 

4/1 

L2* 102 

472 

55 

DO  57  N=L1,L2 

•♦73 

XOsSI-FLOATIN) 

474 

74 


Table  A 1.  Lilting  of  Fortran  Code  LOWTRAN  3 (Cont) 


IF  (VO)  59,56,57 

a 

475 

57 

CONTINUE 

A 

476 

56 

*X(6)*M(9)*C6(N) 

a 

hr  7 

GO  TO  60 

A 

476 

55 

TX(8)*C6(N)*X0*(C8(M)-C8(N-l)t 

a 

479 

TX(6)*MI8)*TX(8) 

A 

460 

68 

SUM*SUtt»TX(8) 

A 

461 

IT(IV.GT.1459»1G0  TO  97 

A 

462* 

GO  TO  69 

A 

463 

MATFR  VAPOUR  CONTINUUM 

464 

61 

IF  (IV.C6.670)  GO  TO  72 

A 

485 

IF  (IV. LT .760  GO  TO  66 

A 

466 

Xl«(V-7ft9,)/50.»i. 

A 

487 

00  63  NN* 1,15 

A 

466 

XM*XI-FLO£T(NN* 

A 

469 

IF  (XM)  66,66,63 

A 

4-?0 

67 

CONTINUE 

A 

491 

66 

TX 1 5 ) «C5( NH1 

4 

492 

GO  TO  67 

A 

493 

65 

TX(S)=CS(NH|*XH*(C5<NM)-C5(NH-1I) 

A 

494 

GO  TO  67 

A 

495 

66 

TX (5 ) = ( V-670 . ) *0  .69 

A 

4 96 

67 

TX(S)=M(5)*TX(5) 

A 

497 

SUN*SUN*TX(5) 

A 

4'J» 

GO  TO  72 

A 

499 

<;•»* 

NITP )GEN  CONTINUUM 

A 

500 

66 

IF  (IV. LT. 2995)  GO  TO  72 

A 

501 

*4=1-346 

A 

502 

n 

SUN«SUM*TX(4> 

A 

504 

GO  TO  72 

A 

505 

c*** 

•**»•  HATER  VA°OUR 

536 

69 

IF  (IV. LT. 12600. AND. IV. Gt. 9875)  GC  TO  70 

A 

507 

IF  (IV. LE. 14520. AMU. IV. GE. 13400)  GO  TO  71 

A 

50  6 

GO  TO  76 

A 

539 

7® 

1*1-135 

A 

510 

GO  TO  72 

A 

511 

71 

1*1-255 

A 

512 

77 

*1*1 

A 

513 

IF  (Mill. LT.l. PE-211  CO  TO  76 

514 

MS1*AL0G1C  (M(l  )MC1 (I) 

515 

IF  (MSI. LT, -2.  3*>66)  GO  TO  76 

A 

516 

IF  ( MS l.GT.7.0)  *1*60 

t 

5 1 6 

IF  (MSI. GT. 3. 5662)  GO  TO  75 

A 

517 

DO  73  ***1,67 

A 

519 

IF  (NSl.LE.FUdtn  GO  TO  74 

520 

73 

CONTINUE 

A 

521 

76 

TX(1)*TR(*)  ♦tTR(*-l)  - T«(*|  l*IFM(*)-MSl)7(FM(K)-FM(*-m 

A 

527 

GO  TO  76 

A 

52J 

75 

TX ( >*8„e 

A 

524 

76 

CONINUE 

A 

5 25 

A 

526 

IF  (IV.LT.666a.AN0.IV.GF.50Q)  GO  TO  77 

A 

527 

IF  (IV. LT. 131)9. AND. TV. GT. 12970)  GO  TO  76 

A 

5 28 

GO  TO  83 

A 

529 

77 

J«I-39 

A 

530 

GO  TO  79 

A 

531 

76 

J* (TV-12956) /5»1516 

A 

532 

79 

IFIK(2).CT.l.SE-28)  GO  TO  63 

A 

5 33* 

*1*1 

A 

534 

MS?*  AL  O'*  16 (M(2))«C2(J) 

A 

535 

75 


Table  Ai.  Listing  of  l-ortran  Code  LOWTRAN  3 (Cont) 


IF  SWS2.1T.-2, 34681  GO  TO  83  A 536 

IF  SWS2.GT .3.5682)  GO  TC  ?2  A 537 

IF  TITS?. G*. 2.8)  K 1 - 4 0 A 538 

PO  8T  K=Kl,67  A 539 

IF  <WS?.|,F  ,FW<«I  | GO  to  J1  A 543 

88  COMTTNUE  A 541 

81  TX  (21  sTRCXl  * r»m<-  II -TRtKII  • »F  M«<)  -WS?I/<FWS<1  -FWSK-lt  I A 54? 

GO  TO  6 3 A 543 

82  TXS?1=C.C  A 544 

63  CONTINUE  A 545 

079NC  A 546 

IF  IfV.i  T. 575. OR. IW.GT. 32701  GO  TO  87  A 5 5 7 

L-I~45  A 547 

61=1  A 549 

IF  <W« Sl.LT. 1, QF-201  G3  TO  87  A 5*9 

WS3  = ALOGIO CV S3  1 1 *C3 H 1 A 553 

IF  <W53.LT. -1.6779)  GO  TO  67  A 551 

IF  SWSJ.GT. 3.9345}  GO  TO  56  A 55? 

IT  SWSJ.GT. 1.51  Ki--36  A 553 

no  84  K = K 1 » 6 7 A 554 

TF  SW53.LE.FO! <(1  GO  TO  63  A 555 

84  CONTINUE  A 556 

85  TX!  31  =TRI*<»-«rilC«l-TH<<-ll  I * <F0«<» -hSJI/ «F0  IK)  -FO  IK-1 ) ) A 558 

GO  TO  67  A 555 

86  TXSJJ=C.O  A 560 

87  CONTINUE  A 561 

c«. AEROSOL  EXTINCTION  A 562 

ALAN=l.9C»%/W  A 563A 

X X = 0 . 8 A 5639 

VY=0.!l  i 563C 

IF  (IHA7E.EO.O.I  GO  TO  93  A 564 

00  88  Nn,44  A 565* 

XO=ALAM-VX(N)  A 566* 

IFSXD)89. 86,88  A 567* 

«8  CONTINUE  A 568A 

69  XX=  (C7SNI -07  IN-11)  «n/!VX«NI -VXSN-ll  » *C7  fN)  A 5b89 

TY=(C7A<N»-C7A<N-:n»X0/!WXINI-ifX(v-l»»*C7A!N)  A 5b8C 

90  TXS10»=YY*W(7)  A 5680 

TXI7I=XX*W(7)  A 569* 

SUN=SUN»TV!7t  A 573 

TXS91=SL‘N  A 5 71 

nc  94  K =4  * 1 0 4 5 72* 

IF  <TX <«> .EO.O.OI  GO  TO  9?  A 5’? 

IF  ITX<K» .LF.O  ,i»  GO  TO  91  A 574 

IF  S * X ? X!  • G ' . r q . ! GO  i u 9‘  A575 

TXCM  rfXP  1-T /S  K)  1 A 576 

CO  TO  94  A 577 

91  Txm  =1  .0-TXfX)*C.5*TX|iO»TX  IK1  A 578 

GO  TO  94  A 579 

92  TXIK1  = 1.0  A 6 3? 

GO  TO  94  A 5S; 

9?  T X 1 K ) =0.  A 552 

94  CONTINUE  A 583 

T»  <1CI  sl.O-r  J!l  101  A 5 o ?♦ 

TXI9) = TX J 1» *TX S?) *TX S T» *T* t 91  A 564 

IF  (IV.GEw>1303Q)  TX(JI=TX!8)  A 555 

IF!  JP.eo.  *»  TX!9>  =TXS7J  A 5S6A 

A9=i.-TX|9l  A 5863 

IF!  IY.n.Ifl.  OR.  I*.  E0.IY21  AB=0,5*A3  A 586C 

5 UNA  r SUM4 ♦ AS  *3  Y A 5860 

IFSJP.E9. 0»  W»ITFI6.4?JI  IW, ALA8, TX ( 91  . 1 T«  (<),<- 1 . 7( , T X! 1 0» , SJNA  A 5 8 7* 
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Table  Al.  LlBtinv  of  Fortran  Code  LOWTRAN  3 (Cent) 


95 


9* 


9T 

96 


160 
Ml 
961 
91? 
96  3 
919 

965 

966 
987 

916 

969 

918 
91 1 
91? 
913 

919 

915 

916 

917 
916 

919 


9?f 

9?  1 
9?? 


9?  3 
9?9 

935 
9 36 
937 

936 

939 
9 33 


IF  (IV.GE.IV2)  GO  TO  95 
GO  TO  58 
*E«0  981,  IKY 
•a-l.*SUR«/V3~  61 

PRINT  9?9,  I Vl.IV3.SJrtA.A9 

PRINT  918, IKY 

IF  I IKY*  EQ.  6 I GO  TO  HC 

GO  TO  (96, 2, 97 ,98. 100), IXY 

RE 60  986,  V1,V?,0V 

AVW*16606. / VI 

ALAH*lf 060./V? 

PRINT  916,  W1,V?,0Y,4L6N,6VM 

50*6*1.8 

GO  TO  99 

IF  ( NOOEL.EQ. 0 ) GO  TO  286 
GO  TO  368 

RE 60  960,HODEl ,IH6?E, I TYPf.LEN.JP. X*,*l, *3 v *3.NL,R0 
PRINT  980, NOOEL. IM4ZE, IT YPE, LE N, JP, IN, Nl, N3.N3.NL, RO 
GO  TO  ?60 

STOP 

FORN6TIiai3.Fia.3l 
FORMAT  (8C10.3I 

FOR *6T  (F6.1.3 (CIO. 3, F6. 1, 3E 18.1 ) ) 

F0RH6T  (9IF6.T.3FT.91 I 
F0RN6T  (15F5.3I 
F0RN6T  (6E9.3I 
F0RN6T  I7F16.3I 

F0RN6T  (//10X, 38K  HORIZONTAL  PATH,  6LTITU0E  =,F7.3,11K  UN, RANGE  *, 
1F7.3.3N  KM) 

F0RN6T  I//18X,  50M  SL6NT  PATH  BETWEEN  6LTITU0ES  Ml  6N0  *3  WHEP.E  Ml 
l = ,c7.3.e«  ttN  H J S.F7.3.16H  K*,?EMIT9  ANGLE  *,F7.3;6M  DEGREES) 
FORNAT  (// 10 X, 39H  SLANT  PATH  TO  SPACE  FROM  ALTITUDE  Ml  *,F7.J.19M 
1KN,  7ENITH  ANGLE  =,F7.3,8H  OCGREES) 

FORNAT  I/30X.16H  MODEL  ATMOSPHERE  .I1.11M 

format  (/Z8X.16H  MODEL  A f mOSPhERE  ,Il,3lH 

FORNAT  (/36X.16H  MODEL  ATMOSPHERE  .I1.31H 

FORMAT  «/>0K.l6H  MODEL  A TN05PHETE  ,I1,31M 

FORNAT  I/38X.16M  MODEL  ATMOSPHERE  .I1.31H 

FORNAT  (/36X.18H  MODEL  ATMOSPHERE  .I1.31M 

FORMAT  (/76X.16M  HAZE  MODEL  .I1.3M  = .A5.13M  VISUAL  RANGE  ) 

FORNAT  (/?5X*MAZE  MODEL  -*,F5.l.*  <N  VISUAL  RANGE  AT  SEA  LEVEL*) 
FORMAT  ( / 1 • X . ? I.H  FREQUENCY  RANGE  VI*  ,F7. 1,139  C*-l  TO  V3*  , FT.  1,1 
19M  CH-1  FOR  DV  *,F6.1,9X  GM-i  (,*6.3,*  - *,E5.2,*  NICRONS  )*) 
FORMAT  I / 1 0 X , 3 SH  FQUILAVENT  SEA  LFVEL  ABSORBER  AN0UMTS//2 IXlljMWAT 
1FR  TiruiR 


TROPICAL) 

NI0LATITU3E  SUMMER) 
Ml  OLA  T ITiJDE  VINTER) 
S JB-ARC TIC  SUMMER  ) 
1963  US  STANDARD  ) 
SUB-ARCTIC  WINTER  > 


LUC  Z 1 U . 


^A  INUWLN  « b WVI  • 


t f*  TUiT  A 


? MOL  SCAT  AEROSOL  0 Z ON-  l’J-W)/29X.  TMGN  CN-3.1CX,  2MKM,  1 

30X.6MATM  CN.13  X,?HKM,9X,7HGN  CM- 3 , 1 0 K , 3HKM,  1 3«  , 3HiCM,  1 1 X , 6H AT N C*l 
FORMAT  I1N1.///16X,*  VERTICAL  PROFILES  * ,69X  ,*PS  J • , bX,  * PHI*  , 6X , * 
*,META».9X,*TMeTA  RANCF*) 

FORMAT  (/18X.9M  M(l-i) *6 (-19.31/1 

FORMAT  ( 1 Ml , /I  OX , 32m  FREQ  WAVELENGTH  TOTAL  M30. 5X9MC03* , 5X , 6 

19H070MF  N»  CONT  M?0  CONT  MOL  SCAT  AEROSOL  AEROSOL  INTEGRATED 
3 /1M.19H  CM-1  MICRONS, 8<9X5MTRANS) .9X.Z6H  APS  ABSORPTION  I 

FORMAT  <16X,I6,16F9.9.F12.?) 

F.1R14T  (*  INT-GRATET  ASOP°riON  FROM*  , I i , * 13*. I>,*  CH-1  **,r18,3, 
1*. AVERAGE  TRANSMTTTA  ICC  =•/!..■.) 

TOP  MAT  (lJX.TFia.I) 

FORMAT  ( /zcx  ,*A£PO jOl  SCATTERING  NOT  COMMUTE 0, 1 9 A ZE  = 6*1 
F 0*  M ! T (in  , ///1CX,  *0M  MOP!  I3NTAL  PROcIlE$/> 

FDTMAT  Hi*,*  Ml** 7. *«W, H3s»,F7. 3 ,*SM, A*GtE* * ,F8.o,  »iEOM.  RANG 
It  **, FT. >,*<*,,?£  T*-T,F8,S,»,  VI  S**,F>. It 
f 0* i»T ( 3rlC . 3,3FS. 1 , 3r  1 J . 3, 'Flu.  3) 

FORMAT ( 13*,* INPUT  HETCDROlOGICAL  0 A T A I • 7 11 X , • Z * • , F 7 . 3 . * <H,  P**,F7 


5tM 

569 

5 SC 

SRl» 

5 91 B 

591C 

591U 

591E 

59? 

593 

599 

5 *5 
596* 
597 
596A 
59PB 
598C 
5?t0 
596E 
599* 
63C* 
b:i 
6:2 

6 3 3 

6 C 9 
635 
606 
637 
639 
639 
610 
oil 
612 
613 

619 

615 
olb 
617 

616 
619 

630 

631 
633 
633 

OL 

635 

636 

637 
63S 
639 
630* 
631* 
633* 
633* 
639  A 
63  5 9 

635 

636 
6 I 7 
63  6 
639 
69  0 * 
6*1* 


77 


9 


? 


Table  Al.  lasting  of  fortran  Code  l.OVVTHAN  3 (Cent) 


i.?,»  N8,T  = *.F=i  .1,*  C,  OFW  PT  ,T£NP»,- 5 . 1,*  C,  REL  HUMIDITY^, F5.1. 
Z*  X,  h?o  Of  NS  I T»:,i  IP'S.  ?•  * GM  N-3»/10X,*  OZONE  DENSITY:*, E9. 2,*  G 
3M-3,  VISUAL  RANGE:*, 0PF6. t, * K N.  R A <3  E = » , r 1 0 . J , • KM  * ) 

431  FCRMATT4TP6.Z,  ZF7.511 

4jz  forhat  c*  starting  parameters  hi  and  anill  have  been  red? mneojhi- 

1 •,f10.3.*ANGl  E :*,F1'.6I 

433  F0?1«T  (•  TRAJECTORT  '•tSS'-'S  EARTHS  ATNOSPMERP.  CLOSEST  DISTANCE  OF 
1 APPROACH  IS»,F10.?,l*./.l*. *END  0r  CALCULATION*) 

434  FORMAT  Y10X, 14 ,*fc. 1, 11 TF1C. 31  I 

4J*  FORMAT  U5,f  7. 1.8E10  . T.4P3.4.F7.  1» 

43*  FORMAT  (*  MMIV  r * , F 1 " . 3 I 

437  FORMAT  <•  PATH  INTERSECTS  EARTH  - PATH  CHANCF D TO  TYPE  2 WITH  H2  = 
1 C.  0 KM*1 

4 3 3 FORMAT  !•  CHOICE  0‘  T HO  PATHS  FDR  THIS  CASE  -SHORTEST  PATH  TAKEN. 

1 TOR  10NGCR  PATH  SET  LEN=1.*1 

439  FORMAT  <•  CHOICE  of  tko  paths  PJR  TJiS  case  -loncfst  path  TAKEN. 

1 FOR  SHORT  PATH  SET  L f N : fl  *1 

440  format  f*  HZ  HAS  SET  L r S S THAN  H N I N ANO  HAS  BEEN  RESET  EQUAL  TO 
1 shin  r.E.  H?  : *,F|B,3» 

441  FORMAT  I * N0DrL  ATM0S°HfP-  NO.  7*,F  4 X,*Z  TKN  *,3X,*P  IMB>»,4X, 

1 *T  TCI  Of*  PT  7.RH  H'0(Gh.M-3)  03TCM.N-3)  NO  DEN.*> 

4<*2  FORMAT  ('  rOG  ICUUOITIJTS  MIT  EXIST  AT  S'l  LEVEL  P0k  THIS  VISUAL  PA 
14G£»./.*  Ip  Su  ThE  , ASS'JH  Th;.  T RA  Nfrll  T T A CE  CJE  TO  EOS  IS  GIVEN 
1 By  T hp  TRANSMITTANCE  AT  C.55  MICRONS*) 

ENO 


A 64?* 
A o4.3* 
A 644* 
A 045* 
A o 46 
A E 4 7 
6 *.9 
6m9 

650 

651 
65? 

653 

654 

655 

656 

657 
65S 
659 
66C 

a 66i* 

A 66?* 

A 66  3* 
A 644* 

A 5o5* 
A 6 66* 


SUBROUTINE  POINT  I X , Y N . N ,NP , Tx , I a ) 

COMMON  7M41, "17, 34), T(7, ’41. E HT  10, 34), WH|7,  34  >,H,NL,  RE  *CW,  CO,  Pi 
DIMENSION  TXT10) 

SUBROUTINE  POINT  COMPUTES  THE  MEAN  RSFRATIVE  INDEX  ABOVE  AND  3 EL  0 W 
A 5 T VrN  AUTIT'JOP  and  INTERPOLATES  E X PONENT I ALL  T TO  DETERMINE  THE 
EQUIVALENT  ABSORBER  AMOUNTS  AT  THAT  ALTITUDE. 


X IS  THf  HEIGHT  IN  QUESTION 

TXI9)  4 NO  Y.N  ARE  ThF  MEAN  REFRACTIVE  INDICES  ABOVE  AND  8EL0H  X 
N IS  THF  LFVEL  INTESFR  CORRESPONDING  TO  X OR  THF  LEVEL  BE LON  X 
N"  - 1 *P  X CO  INCI9E5  Mi»H  MODEL  ATMOSPHERE  L-VEL  , IE  NOT  N P = I 
TXT  1-8 ) ARE  ABSORBER  AMOUNTS  "E  R <«  AT  HEIGHT  X 

N-NL 

NP-0 

IF  TX.  IT.  0 . 0 I * = 0. 

IF  TX.GT.7(NL» I GO  TO  4 

00  1 Isl.NL 

N-T 

IE  TX-7TD)  7,4,1 

Continue 

J?:N 
N = N-  1 

FAC- TX-’TNl l/T  Z I J2 1 - Z TN)  I 

P*l  = PTM,N)*fPTN,J?)/PTN,N)l  **EAC 

Txi=T(M,NI»TTIM,j?(/rTM,N«)»*PA: 

MXlsNHTB,  N)»TMHTM,JZ)/HH|NtN)l»»PAC 


9 

3 

9 

3 

3 

3 

B 

3 

B 

3 

3 

3 

a 

a 

9 

B 

B 

B 

8 

3 

B 

B 

8 

B 

9 

9 

B 

3 

B 

a 


i 

3* 

3 

4 

5 

6 
7 
« 

9 

x c 

ii 

i? 

I 3 
IU 

1 5 
16* 
17* 
IS 
19A 
19B 
20 
21 
22 

2 3- 

25 

26 

27 

28 

29 

30 


7 a 


^ i'V^e Ir***-1***-  '. 
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Table  Al.  Listing  of  Fortran  Code  IjOWTRAN  3 (Cont) 


* 


TX(  3»=C0*PXl/TXl-4.56F-6<1HXl»rXltCW 

TX(21=C0*P(N,J2)  /TIM,  J2I  -W^eE-a^KM  M,  J2)  ®T(M,  J2)  *CW 

TX(1) =CO*P(H,N)/T(M,N)  -4 . 56E -6*WH IN,  N>  *T (H,NI*CW 

TXI9I  = 0.5F-6MTX(2f«-TX(3l) 

rN=c,FE-6*irxii5  4-rx(3n 

IF  IT' .EQ.O > 50  TO  9 

00  3 K = l,8 

TXIKI =0.0 

IF  IFH(K,N) .EQ.O.C)  GO  TO  3 
IF  <EH(K,N)  .GT.100Q  . 0)  CO  *0  3 
TXIKI =EH(K,N>* (£H,K, J2I /EH(K,N)> **F4C 
x CONTINUE 

GO  TO  9 
4 NP=i 

IF  IIP.FQ.O)  50  TO  6 
00  5 K = t,  9 
6 TXIKI =EHIK,N) 

6 TX(9t=FH(9,N)-l. 

YN  = 0 • 0 

CAPOS  B 24  A NO  50  THROUGH  59  ARE  NO  LONGER  REQUIRED 
If-  IN.GT.ll  YN  =FH  !9 . N-ll  -1 . 0 
9 CONTINUE 

IF  (IP.EQ.il  PRINT  402,  X, N, NP, T X HI , VN, IP, I TX I Kl , K = 1 , »> 
TXI9I =TX( 91 f 1. 

YN  = YN«-1, 

RFTUPN 


c 

400 


FORMAT  (/,»  FROM  POINT!  HE I 5 HT =* , F 1 3 . 4, * KM, N= * . I 3 , * , NP=* , 12 , * , RE F 
1.  TNOFX  ABOVE  0 HFLOW  X**,  2E 11  .4,  *, I P = * , T7 , / , 1 2K , *EQl'I sf. 
zamjunts  pfr  km  at  x=*,8Ei:.3> 

ENO 


B 

B 

8 

8 

B 

8 

S 

a 

B 

0 

B 

a 

B 

B 

B 

B 

B 

3 

8 

a 

B 

B 

8 

5 
B 
8 

6 
8 


31 

32 

33 

34 

35 

36 
374 
37B 
3* 
39 
4 C 

41 

42 

43 

44 

45 

46 

47 
49 
40f 
49 
6D 
61 
62 

63 

64 

65 
65 


B 68 
3 69 


C* 

C 

C 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  AHGl.  I HI  , H 3 , A NGLE  ,B1  »LEN,  ML  ) 

COMMON  2(341  ,P 17,34) ,T (7,341 ,EHI 10, 341 ,MH 17,34} ,H,NL, RE,QM,CO, PI 
DIMENSION  TX(lS) 

*»¥»»¥»» ¥¥¥¥¥¥¥*¥' 

THIS  SUBROUTINE  CALCULATES  THE  INITIAL  ZENITH  ANGLE  (ANGLE) 
TAKING  INTO  *CCOUNT  REFRACTION  IFFE3  TS  GIVEN  H1.H2,  AND  BETA 
(WHERE  BETA  IS  THE  EARTH  CENTRE  ANGLE  SUBTENDED  BY  HI  AND  H2  ), 
ASSUMING  THE  REFRACTIVE  INDEX  T3  0F  CONSTANT  IN  A GIVEN  LAYER. 
h OR  GREATER  ACCURACY  INCREASE  THE  NJMBER  OF  LEVELS  IN  THE  MOOFI. 
ATMOSPHERE, 

THIS  SUBROUTINE  CAN  BE  RE'OVED  FROM  THE  PROGRAM  IF  NOT  REQUIRED. 

I ,* ,7 , « 

IF‘=99 

CAsPI/180 . 

X1=RE*H1 
X2=RE»H2 
LFN=C. 

IT-0 

BI=B1»CA 

IF  (Bl.LQ.C.O)  BIrACDSIKZ.  XI.* 

Y AKG*X2*SIN( Bt )/(X2®C0S(B) » -XII 
THET=ATAN(TANG) 

I F (THFT.LT, o. 01  THET=THCTfPI 
SP;i  I = S I N(  THEY) 


c 

c 

C 

c 

c 

c 

c 

c 

r** 

c 

c 

c 

c 

c 

c 

c 

c 


1* 

2* 

3 

4 

5 

6 

7 

8 
3 

1C 

11 

3.2 

13 

14 

15 

16 
17 

A.  1 ’ 

19 

20 
21 

2 IB 
22 
23 
r.  4 
25 
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Table  Al.  Listing  of  Fortran  Code  LOWTHAN  3 (Com) 


ANG=THET/CA 

PP 1 NT  404,  8i  , ANG  , T 4 N4 
TN-THF T 

tm=tn-p .s»ca 

anglf=thft 

FPT=C. 

«FTA=0  . 

6FT1-0 

BfT2=m 

FpTl-0 

fpi?=o 

FpT  3 = 5.0 

IF< Ri.Le.I.O)  GO  TO  z 
PPlNT  400,  IT 
r=?.*THFT 

IF  (Y-PI.GT.l.  OF-fl)  GO  TO  9 

IF  (IP.E0.100)  GO  TO  h 

XHTNrYJ^COStRl l-PE 

IF  CXMIN-H1)  0,4,4 

HMIN=H2 

H2  = H1 

HI -MIN 

BNGLF-C .S*PI 

THFT-ANGLr 

SPHI=i.O 

ANG=BNGLF/CA 

PPINT  404,  B1,ANG,3PHI 

IP=100 

CALL  POINT  (H1,TN,N,NP,TX,IP» 

J1  = N 

TX1=TX (9) 

CALL  point  (H?,YN,N,Np,TX*I0) 

IF  (NP.FQ.l)  Nr.N-1 
J2*N 

IF  (Jl.EQ.J2t  TXJsTXl»TN-FH(9, Jl ) 

oo  y j=ji,jz 
XI =RF+7(J> 

X2  = 9F ♦ 7 ( Jr  1 1 
IF  (J.FQ.J1)  Xl-REt'31 
IF  (J.FQ.J2)  X2-PF*N? 

SALd"X1*SPHI/X2 
AL  r l SGIp1  i 5 4LPi 
PN=fH(9,Jei)/rH(9,J) 

IF  UJUl.EO.J?)  9N-TN/FH(9,  J) 

IF  fJ.FQ.Jl)  PN=FHt9, J*1»/TX1 
IF  < (J*l)  ,tO.  J2.  AND.  J.EO.  Jl  ) RN-YN/IX1 
BF  T = THE  T- A L P 
FP  = -TAN  (ALP) 

IF  (J.NE.J1)  F9-FBHANCTH- T1 
FBV=F<JT*F8 

bfta-pfta+bft 

TH1-THFT/CS 
et  = BF.T/CA 
C* AL  P/C  A 

PRINT  402,  J,7«J»,TMPT,ALP,BFT.lFTA.reT,FH,THl,3f,C 
IF  «X?,FQ.RF*H2)  C-PI-ALP 
IF  (GALP.GE.PN)  PN-1. 

GPH I = S A LP/P  N 

vhft  = asinc gpht i 


0 Zb 
C 27 
r ?4 
C 29 
Z 3 0 
C Jl 
C 37 
C JT 
Z 34 
4 35 

C 14 
C V 
C 3 7* 
C 33 
Z 79 
: 4: 

C 41 
C 47 
C 4 3 
C 44A 
C 4 49 
C 4.0 

C 49 
3 9n 
r;  4’ 
c 4 3 
c 4 9 

C 9 3 
: si 

c r ,? 

V 3 ! 

3 Zu 
C S3 
C 59 


j > ■ 

C 5 9 

C 5 9 


C L 3 
<Z  64 
0 65 

c be 
c er 
c be 

C 4 5 

c ’0 
0 71 

C 7)’ 
3 7 T 

C 74 
C 75 
C 7 f 
3 7 7 

C 7ft 
3 7 4 

3 A 


Table  Al.  Listing  of  Fortran  Code  LOW  TRAN  3 (Cont) 

7 CONTINUE 

TF<B1<LE<0<91  GO  TO  29 
GO  TO  26 
6 CONTINUE 

TANGs-TANG 
ANGLE=FI» ANGLE 
TN=ANGLE 
ANG=ANGLE/CA 

C PRINT  40*,  01, ANG.TANG 

IF  IMl.LE.OcO*  GO  TO  3 

9 CONTINUE 
IP*t01 

CALL  POINT  <Hl,YN,N,NPl,Tx< IP1 
TXi*TX (9) 

VN1*YN 

IF  (NP1.EQ.1I  N*N-1 
J2*NL 

IF  (N.EQ.n  J2*NL 

Ji*N 

J*Jltl 

IF  (H2. GE. HI I GO  TO  13 
CALL  POINT  <H2,YN,N,NP,TX,IP) 

TX2*TX ( 91 

VN2=YN 

J2*N 

IF  <J1<EQ<J25  TX2=VNltTX<9» -EH<9, Jl) 

10  J«J'l 
X1=9E+2(J*1» 

X?*RE+7<J! 

1~  iJ.eQ.Jl>  Xl=RFFMl 
IF  <J<E<J<J21  X2sRE»H2 
SALP=X1*SPHI/X2 
H*IN*Xl*S?Hl-RE 

C PRINT  *02,  J,X1,2<JI,SPHI,SALP,HHIN,RE 

IF  <SALP.LF.l. 01  GO  TO  11 
SAt P.SPHI 

IF  <mhin.GT.H21  GO  TO  18 

11  ALP*AblIN(SALP» 

TNET=ASIN<SPHl) 

0EI»ALP-THET 
0ETli0ETi»8ET 
FBaTANf ALP> 

IF  < J.NE. J1 ) r8^F8-rAN!THETI 

FBTliENTI  *FO 

TH1=TMFT/CA 

0E=BET/CA 

AL*ALP/CA 

C PRtNT  *02,  J,X2,THET,ALP,BFTl, BET, BMINjHNIN, Foil, THl,8E,«l 

IF  tX2.EQ.RE<-H21  C*PI-ALP 
REF*EH  <9, J) 

IF  <J.EQ.JU  R£F*YN1 
IF  IJ.CQ.J2)  REF*TX2 
IF  < J.EQ< 1 } GO  TO  12 
PN«EH(9, JJ /EHI 9, J-ll 
IF  <J,EQ.J1>  RNsYN1/EH<9,,1-1) 

IF! J.E0.J2frll  RN«REF/TX2 
IFCJ.EQ.J21  RN.RCFFYN2 
IF  <SALP.GE.RN1  RN«1. 

SPH I*SALP*RN 
IF  (IIJI.LE.H?}  GO  TO  12 


C U 

C Bit 
C iZ 
C 8 3 

u 8* 

C 35 

V c 6 
C 87 
C 98 
C 39 
C 90 
C 91 
C 92 
3 93 

C 94 
C 95 
C 96A 
C 960 
C 97 
C 98 
C 99 

c ica 
c 121 

C 102 
C 123 
C 134 
C 105 
C 106 
C 12  7 
c ioe 
C 139 
C 110 
C 111 
C 112 
a ii3 
C 114 
C U5 
C 116 
C 117 
C 118 
C 119 
C 120 
C 121 
C 122 
C 123 
C 124 
C 125 
C 126 
C 127 
C 128 
C 129 
C 130 
C 131 
C 132 
C 1J3A 
C 1338 
C 133C 
C 134 
C 135 
C 136 
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Table  A 1.  Listing  of  Fortran  Code  LOWTRAN  3 (Cent) 


GO  TO  19  0 13? 

12  XI -*2  C 138 

IF  (A8St?(JI— t2t.LT.i.OE-iO. ANO. J.NZ.l)  GO  TO  13  C 139 

GO  TO  14  C 140 

13  J-J-l  C 141 

X1=PETZIJ«TJ  C 142 

IF  tJ.ea.Jl)  XliRF+Hl  C 143 

IF  (J.EO. J2. ANO.J.NE. JU  X1=RE»H2  C 144 

14  X2=RFFZ(J1  C 145 

HMIN=X1*SPHI-<?E  C 146 

IF  (F'MIN.Lf  .0.  0)  GO  TO  25  C 14? 

IF  I7IJT.LT. HMIN)  GO  TO  13  C 148 

RFF=FH<9,J1  C 149 

IF  (J.F0.J2)  PFF--YK  C 15C 

SALP=X1»SphT/X2  C 151 

ALP=AGIN<SALP)  C 152 

THET=*5IN(SPHT > C 1 r 3 

9ET=ALP-TMET  C 154 

FB=TAN(ALPI -TANtTHEn  C 155 

FBT2=FBT2fFM  C 156 

BET2=9FT2+BFT  0 157 

BMIN=BF  fl  I-8ET2  C 15e 

AL=AtP/CA  C 159 

TH1.  = TMFT/CA  C 1 6 c 

C PRINT  432,  J,X2,THET,4LP,9FT2, 827, SHIN, HHIN,reT2,THl,9E,/L  C 161 

PN-9Er /FH  C9 , J- 1>  C 162 

TF  (SALP.GE.R4)  RN=1,0  C 163 

SPt;I  = 5fiLP'RN  C 164 

GO  TO  13  C 165 

17  TX3*YN1*TX(9) -EH(9, J1 J C 16b 

YN1 5 TX  3 C 167 

IF  tA8S<H2-7fJ*l)>.LF.1.0r-5)  Y9i=TX<9)  C 16* 

TF  fBBSfHi-ZIJMU  .LF.i.OE-5)  Y91^TXC9)  C 169 

BN- 1 . 0 C 17C 

GO  TO  19  C 171 

18  CALL  POINT  (HNIN,YN,N,NP,TX, IP)  C 172 

IP  = 1 0 2 0 173 

TXJzTXI9)  C 174 

IF  tJ.EQ. J1.AN0.M2.GE.H1 ) GO  TO  17  C 175 

IF  ( J.fQ. J1.0R . J.EQ. J2)  TX 3 YNZ* T X< 9 ) -EH  (9 , J ! 0 176 

IF  tMMIN.GT.H2)  TX3=TX»9)  C 1?7 

IF  CJ.t'j. Jl.ANO.HHIN.GT.H2)  GO  TO  17  C 17s 

RN*REF/T)T3  C 179 

IF  C S ALP. GF • RN ) KN  = 1.  C 18C 

SPHI=SALP*RN  Z lsl 

X=X1*SPMI-PE  G 182 

0TF=A8SfHMIN-X)  C 183 

MMI N~  X C 184 

IF  (OrF-l.OE-5)  19,19,18  C 155 

19  X?s»E  *HMIN  C 186 

C PRINT  403,  HHI N, DIF , RN  0 187 

TMET=ASINISPHT > C 188 

IF(RN.FO.l.O)  F3T3  = .-TAN(THFT)  C lct*B 

IF  (PN.E9.1)  GO  TO  20  C 189 

DNX=ITX3-1.C)*AL0GI  (5  KT-l.O ) /(REF-1.  C ) )/ TX2-XI  ) C lxT 

F8T3  = ~TANtrHFT)M1.0-t.0/<l.QTTX37<X2»DNX)>>  0 191 

20  8FTs0.5*PT-THFT  Z 1 iZ 

BET2*BE 1 2 »8FT  C 193 

BHIN^BETI *BET?  c 194 

IF  (M2.GE.H1)  GO  TO  23  C 195 

8ET=BET1*2.*BET2  C 196 

OB1-B1 ~6E  T l 

H2 


Table  Al.  Listing  of  Fortran  Code  LOWTHAN  3 (Cont) 


C 197 

DB2*8FT-B1  C 198 

21  0B3  = A8StBHlN'*Bl>  C 199* 

IF1UP3.GT.0B1.AND.DB2.GT.QB1>  GD  TO  25  C 1*38 

TFt0B2.GT.0B3>  GO  TO  22  C 1990 

1FIOR2.GT.OB1I  GO  TO  25  C 200 

BETAxRfT  C 231 

FBT=F8Tie2.0*(FBT2*F8T3>  C 202 

LEN*1.  C 2C3 

GO  TO  26  C 204 

22  eer*=RFTl»6ET2  C 235 

FBTsFBTleF8T2*FBT3  C 206 

C PRINT  401.  J.9£T*,PBT,FaTl,FBT2,cBT3»TXl,TNl  C 20 7 

GO  TO  26  C 20o 

23  BFT»=2.0*<BET1+BET2>  C 239 

LEN  = 1 . C 210 

FBT=2.fl*tFBTi*F0T2*-FBT3)  C 211 

PRINT  401  * J.BET*,FBT,F8Tl,FBT2,FBrS*TXl,TNl  C 212 

IF  tH2.Eq.Hl>  GO  TO  26  C 213 

IP=10!  C 214 

IF  t MPJ  .EQ. 1 ) J1*J1M  C 215 

SPHl=SlNt  ANGLE  > C 216 

IF  tZ(JlMI.LE.H2)  GO  TO  24  C 217 

RNxTXl/VNl  C 218 

IF  (SPHI.CE.KN)  PN=1.  C 219 

$RHT=S2HI/RN  C 220 

ThETxASTN  tSPHI ) C 221 

GO  To  5 C 222 

24  CALL  POINT  IHZ.VN.N.NPfTX.IPI  C 223 

TXlsTXlMN-EMt9,  Jl>  C 224 

RNsTXl/'IN.l  C 225 

J2--JJ1  C 226 

IF  t3PHI.GE.RNt  RN=l.  C 227 

SPMTxSPHI/RN  C 228 

THET*ASIMISPHI)  C 229 

GO  TO  5 C 230 

25  8FTMBFT1  C 231 

LEN*0.  C 232 

FBT*F9T1  C 233 

26  THE V= ANGLE* { 81 “SET Ap /( i • *F5T/F *NG»  C 234 

08ET*  = BET A/CA  C 235 

BxBETl/0*  C 236 

TH1*THET/C»  C 237 

PRINT  404,  BET*,08ET*,FBT,TH1,T*NG  C 238 

IF  tTHET.GT.TN.OR.THET.LT.TH)  THE T= f TN»TM> /2 . C 239 

TH1«THET/C*  C 241 

PRINT  404,  8ET1.B.F8T ,TK1  C 240 

TN1*TN/C*  C 242 

THlxTHZC*  C 243 

PRINT  405,  TN, TN, TNI, TNI  C 244 

SPHI*SIN5THETJ  C 24  r 

TANG*TAN(TM£T>  C 246 

ITxJTM  C 247 

0BE=*BStBl-BET»)  C 248 

OTM«*BSt*NGLE»THET)  C 249 

IF  (IT.EQ.10)  THETx0.5*tANGLE*TMET>  C 250* 

IF  tIT.EQ.lO>  GO  TO  29  C 251 

IF  tO9E.GT.1.0E-7.*NO.OTH.GT.l.»E-7)  GO  TO  1 C 252 

28  ANGLExTHETZC*  C 253 

PRINT  406.  4NGLE.IT  C 254 
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RETURN  C 255* 

29  HJ-H2  C 2559 

ANGLF=C/CA  3 2 55C 

PRINT  *06,  ANGLEtU  C 2550 

RETURN  C 25DL 

C C 256 

*80  FORMAT  I//*  If  TER#  'ION  NUMBER  C 25? 

*C1  FORMAT  416, £16  .7,8*13.8)  C 258 

*62  FORMAT  4'*,F10  .*,6E13.*,4rl!)  .*/)  C 259 

*03  FORMAT  I*  HMIN=*,F1%.6,*  0XF=*Ei*.6,*  PR^»,El6.8>  C 26C 

* 0*  FORMAT  |*  TOTAL  SF  T A = * , r 1* .6 , f l 5 . & , * , F BT  = *«E1*.6,»  THET  =*,c10  C 261 

1.6.*TANG=*,F10 .6,*)  C 262 

*05  FORMAT  45*12.6)  C 26? 

*06  FORMAT  <8*, /JH*, ’ZENITH  ANGLE  =*,*2.3.*  DEGREES  « RECOMPUTED  5 26* 

1 FROM  SU0ROUTINE  ANGL  4 r TT£R AT  ION*. I 3 , •» •>  C 265 

ENO  C 266 


The  input  data  given  in  Table  A2  can  be  summarized  as  follows: 

(11  The  first  card  gives  the  number  of  model  atmorpheres  IATM  to 
be  read  in  and  the  number  of  levels  NL. 

(2)  The  next  six  cards  contain  trie  haze  number  densities  (HAZE  1 
and  HAZE  2>.  Units:  cm  3. 

(31  The  next  102  cards  contain  the  model  atmosphere  data  for  the  six  geo- 
graphical models  (with  two  models  for  each  altitude  on  one  card).  Units: 

• _a 

altitude  (km',  pressure  (mb),  temperature  ( K),  H.,0  den  ity  (gm  m 1, 

- 3 * 

density  (grn  m 1. 

(4)  The  next  11  cards  contain  the  aerosol  extinction  (C7>  and  absorption 
(C7A)  coefficients  as  a function  of  wavelength  VX.  The  o’-der  and  units 
are  as  follows:  VX  <um>,  07  (km'  'l  and  C7A  (km  '). 

(5)  The  next  17  cards  contain  the  transmittance  scale  (TR)  and  logarithmic 

scaling  factors  for  water  vapor  and  the  uniformly  mixed  gases  (F\V>,  and 

_ a 

for  ozone  (FO).  Units:  H9<  > (log  jq  gin  cm  ),  uniformly  mixed  gases 
Hog  kml,  <)„  (log  jq  atm  cm)  respectively. 

(6)  The  next  344  cards  contain  the  spectral  data  for  the  ’'arious  molecules 

in  the  following  order:  O 1 (l^O),  02  (uniformly  mixed  gases),  C3  (ozone), 
C’4  (nitrogan  continuum),  05  (HjO  continuum),  08  (UV  ai  d visible  ozone). 
The  parameters  Ol,  02  and  03  are  in  the  form  of  logariti  mic  absorption 


coefficients,  with  units  log  jq  (gm 


•1 


1. 


- 1 2 - 1 

cm  ).  log  jq  (km  ) and  log 


(atm  ' cm  ‘I,  respectively.  The  units  of  04  to  OH  are  km 

(7)  The  last  four  cards  contain  the  operational  instructions  for  executing  the 
program  and  are  discussed  in  detail  in  Section  S. 

A wavenumber  (cm  ')  ident.fication  is  given  in  the  last  5 columns  of  the 
spectral  data  cards  described  above.  Labels  identifying  the  various  card  groups 
are  also  given  in  Table  A2. 
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Table  A2.  Listing  of  Data  for  IXOWTRAN  3 


6 3*. 

2.8306*03  1.245E*8J  5.3746*82  ?. 2576*32  1. 
6.0736*0 1 f;.  8226*01  5.6796*01  S.320£*3i  5. 
4.514E*81  4 < 460E*01  4.317E*01  ’.6366*31  2. 
8 • S3 16*00  ’.4346*80  2.2396*00  5.893E-01  t. 
1.970E-08-0. 

1.3796*04  5.0346+83  1.0456*03  6.7356*02  2. 


8.0 

1.0136*93 

380  .0 

1.9E  01 

5. 56-05 

1.0 

9.0406+02 

294.0 

1.36*01 

5.36-05 

2.0 

8 . 050F ♦0? 

288.0 

9. 36*05 

5.4E-05 

3.0 

7. 150E*02 

284.0 

4,7E*0: 

5.1E-35 

4.0 

6. 3306*02 

277.0 

2 . 2 E ♦ 0 0 

4.7E-05 

5.0 

5 .590E  *02 

’70.0 

1.5E*00 

4.5E-05 

6.0 

4,9206*02 

264.0 

8.5E-C1 

4.36-05 

7.0 

4. 323E  *02 

257.0 

4.76-01 

4.16-05 

8.0 

?.78?E*0? 

250.0 

2.5E-C1 

3.9E-35 

9.0 

3. 2°0E*0? 

244.0 

l.’E-Gl 

3.36-05 

10.0 

2.8606+0? 

237.0 

5.0E-02 

3.9E-05 

11.0 

2.4706+02 

230.0 

1.72-0? 

4.1E-05 

12.0 

2. 1 JOE  402 

224.0 

6.0E-C3 

4.3E-J5 

13.0 

1.8206*02 

217.0 

1.8E-03 

4.5E-35 

14.0 

1.5606  *0? 

210.0 

1, 06-03 

4 .5E-05 

15.0 

1 . 320  E 402 

209.0 

7.6E-04 

4. re-05 

16.0 

1.110E402 

197.0 

6.4E-G4 

4.7E-05 

17.0 

9.  3?G £+0 1 

195.0 

5.66-04 

6. IE-05 

10.0 

7.B90E401 

193.0 

5.06-04 

9.0E-05 

19.0 

6. 660  E 40 1 

203.0 

4.9E-04 

1.4E-04 

20.0 

5.650F401 

207.0 

4.5E-04 

1. 96-04 

21.  t 

4.800E401 

211.0 

5.1E-04 

2.46-04 

22.0 

4. 090  6 40 1 

215.0 

5. IE -04 

2. 86-04 

23.0 

3.500E401 

217.0 

5.46-04 

3.26-04 

24.0 

3.000E401 

219.0 

6,06-04 

3.46-04 

25.0 

2.570E401 

221.0 

6. 76-04 

3.4E-Q4 

30.0 

1.2206+01 

232.0 

3.66-04 

2 .4E-04 

’5.0 

6. OOOE490 

243.0 

1. IE-04 

9.2E-05 

40,0 

3. 05  0 E 400 

254.0 

4.3E-05 

4.1E-05 

45.0 

1. 590E  *30 

265.0 

1 • 9C " 0 5 

1.36-05 

?Q.O 

B.540E-01 

270  .0 

6. 76-C6 

4.3E-35 

”0 . 0 

5 . 7 9C6-62 

2i9  .0 

i.4t-l)7 

8 .56-39 

100.0 

3.0006-04 

210.0 

l.CE-09 

4.36-11 

99999. 

210. 0 

0.0 

1,0186+03 

272.2 

3.56+00 

6.0E-05 

1.0 

8.973E*02 

268.7 

2.56*0: 

5,46-35 

2.0 

7,997C4fl? 

265.2 

1. 86*00 

4.9E-05 

3.0 

6.9386*02 

261.7 

1.26*00 

4.9E-35 

4,0 

6, 08  IE  *0  2 

255.7 

6.6E-C1 

4.9E-05 

5.0 

5. 3136*02 

249.7 

3.8E-01 

5.8E-3S 

6.0 

4.6276*02 

243.7 

2.1E-C1 

6.46-05 

7.0 

4.016F40? 

237.7 

8.5E-02 

7 ,76-05 

8.0 

3. 473E  *02 

231.7 

3.5E-02 

9.0E-05 

9.0 

Z.992F402 

225.7 

1.66-02 

1 . 76-34 

10.0 

2.5686*0? 

219.7 

7.56-0 J 

1 .56-34 

11.8 

2.1996*02 

219.2 

6.9E-03 

2 .IE " 04 

12,8 

1.8826*02 

219.7 

6.0E-J3 

2.6E-04 

13.0 

1.610E*02 

210.2 

1.86-0* 

3. 06-14 

14.0 

1-3786402 

217.7 

l.GE-03 

3. 26-04 

15.0 

1.  1286*02 

217.2 

7,66-04 

3.46-84 

16.0 

1.0076*02 

216.7 

6,46-04 

3.5E-04 

17.0 

8.6106*01 

216.2 

5.6E-C4 

3.96-04 

18,  6 

7.3506*51 

215.7 

5.36-04 

4.1E-54 

19.  • 

6. 2806*81 

215.2 

4.96-04 

4.3E-04 

1936*02  8. 

9926*01 

6.3416*01 

5.8936*01 

5 69E*C 1 5. 

150fc*Ci 

5, C 52E*  C 1 

4.7476*01 

6636*01  1. 

935E*C 1 

1.456E*01 

1 • l 1 4E*0 1 

551E-01  4. 

064E-02 

1,0786-02 

5.5  5 3E-05 

4546*02 

1.0136*32 

294.0 

1.46  .1 

6.3E-05 

9.0206*02 

29Q.3 

9.36*30 

6 . 3 £ -3  5 

-.020E+02 

285.0 

5. 9E*0  0 

6 • 0 E -0  5 

7 ,100E*02 

279.0 

3. 3E*0  0 

6.  ? E-u  5 

6.28C6+02 

273.0 

1.9E*0 , 

6 » 4 E -0  5 

5.5406*02 

26’. 0 

i.oe*o : 

6.5E-05 

4 .6  7CE*  02 

261.0 

6.1E-01 

6.  iE-05 

4.2606*0? 

255.  G 

3.76-31 

7.3E-35 

3.72DE*ii2 

248.0 

2.16-31 

7 , IE -05 

3.2406*0? 

242.0 

1.2E-J1 

9.5E-05 

2 .0 1 CE*  02 

235.0 

6.46-02 

9. 3E-05 

2 .430E*  0? 

229.0 

2.2E-02 

1.1E-04 

2 .0906*02 

222.0 

6,06-33 

1 . 2 E -0  4 

1.79EE+02 

216.0 

1.6E-03 

1.  5£ -04 

1 .5306*02 

216.0 

1. 86-0  3 

1.8E-04 

1.3006*02 

216.3 

7.6E-04 

l.JE-04 

1 . li0E*0? 

216.0 

6.46-34 

2. IE-04 

9.5CCE+01 

3i  f r 
fc.  4.  Ki  « u 

5 . 6E-3  ** 

2, 4£  -3  4 

6 .1206*01 

216. C 

5.UC-04 

2.SE-04 

6.95C  E*  01 

217.0 

4.9E-04 

3.2E-04 

5 .950E*  01 

218.0 

4.5E-0-. 

3.4E-J4 

5.1006*01 

219.0 

5.16-04 

3. it  - J 4 

4,3706*81 

220,0 

5.1E-04 

3 • 5 E - 0 4 

3.7636*01 

222.0 

5.46-04 

3. « 6 -04 

3.2206*01 

223.0 

6.06-34 

3.2E-04 

2.770E*0i 

224.0 

€ ,76-Ow 

2.0E-04 

1 .3236*  ul 

2 34.0 

3.  66-3<* 

2.3E-0  + 

6.520E*  00 

245.0 

1.  16-34 

9.2E-05 

3.330E+GC 

258.0 

4.  36-35 

4. IE-05 

1.760E*00 

270.0 

1. 96-05 

l.JE-05 

9.51CE-01 

276.0 

6.3E-36 

4, 3E-.6 

6.71PE-0Z 

213.3 

1 . 46  ~37 

o . 3E -0  , 

3 ,0  COE- 04 

”C.C 

1.0E-39 

4.3E-U 

210.0 

* 

1 .010 E *03 

287.0 

9.16*0 

4.9E-05 

6.960E+02 

202.3 

6.06*03 

S.4E-05 

7 .929E*  02 

276.3 

4. 26*00 

5.5E-05 

7.0006*32 

271.0 

2.7e*’J 

5.5E-05 

6 . 16CE»02 

266.3 

1 . 76*  0 J 

6.0E-05 

5.41CE+02 

260., 

1.0£*J. 

6.46-05 

4 .7  30E*  02 

253.0 

5.4E-I!  1 

7.1E-05 

4 . 1 30E*02 

246.0 

2.9E-01 

7.56-05 

3 o 59BE*  02 

239.0 

1.3E-31 

7.96-05 

3.1 C7E*0? 

23?. 0 

4.26-32 

1.16-04 

2.6?7E*0? 

225.0 

1. 5* -32 

1.36-04 

2.3006*02 

225.0 

9.4E-33 

1,86-04 

1 . 9776*02 

225.0 

6.0E-8J 

?. 16-84 

1 .7006*82 

225.0 

1.6E-B3 

2. 3 £-04 

1 .460E*0? 

225.  J 

1.06-3  3 

2.86-04 

1.250E*02 

225.0 

7.6E-34 

3.26-04 

1 .0  886*02 

225.0 

6.4E-34 

3.46-04 

9.2806*01 

225.0 

5.66-34 

3. >i -34 

7 .98CE  » 01 

225.0 

5.06-34 

4.1E‘'2  * 

6 .860E*  91 

226.0 

4.96-04 

4. IE-04 

8 

a 

« 

r* 

tn 
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MODEL  ATMOSPHERES  i & 2 MODEL  ATMOSPHERES  3 


Table  A2,  Listing  of  Data  for  IjOWTRAN  3 (Conti 


20.0 

5. 370E*0 1 

215.2 

4.5E-04 

4.56-04 

5.8956*01 

225.3 

4. 56-04 

3.92-04 

21.0 

4.5606*01 

215.2 

5.16-04 

4.36-34 

5.02r£*il 

22  8.0 

5.16-04 

?.Sr-04 

22.0 

3.9106*01 

215.2 

5. 16-04 

4.3E-04 

4 .3  60  E *01 

22  6, 

5. 1E-Q4 

3.22-3- 

23.0 

3.3406*01 

215.2 

5.46-04 

3.96-04 

J.25CE*01 

225.  u 

5 , 46 -0  * 

.3. 02-34 

24.0 

2.0606*01 

215.2 

6.36-04 

3.56-04 

3.227E+01 

22 6. -1 

6.06-34 

2.3E-J4 

25.0 

2.4306*01 

215.2 

6.26-04 

3.46-04 

Z.T8C£*01 

228,0 

6. 76-3* 

2.52-34 

30.0 

1.1106*01 

21?.4 

3.6E-04 

1 .96-04 

1 .3426*01 

235,  J 

3.66-34 

1.4E-J4 

35.0 

5.1006*00 

222.0 

1. IE-04 

9 .26-05 

6.6106*0: 

247. c 

1. 16-24 

9,26-05 

40.0 

?.530E*00 

243.2 

4. 36-05 

4 . 1 fc  -3  5 

3.4006*80 

262.3 

4. 36-05 

4.12-35 

45.0 

1.2906*30 

258.5 

1.96-05 

1.3E-05 

l.filC6*8C 

274.0 

1. 96-05 

1. 32-05 

50.0 

6.0206-01 

265.2 

6.36-06 

4.SE-36 

9.8206-31 

277.3 

6. 36-06 

4.32-36 

rs.o 

4.6206-02 

230.2 

1.4E-C  2 

8 .56-08 

7.020E-D? 

216.3 

1.46-37 

a. 52-0 " 

100.0 

3. 0006-04 

219.2 

1.06-0 < 

4.3E-11 

3.0C66-34 

213.0 

1.06-09 

4,  3 2-11 

99595. 

210.0 

21C  .3 

* 

0.  0 

1.0136*03 

252.1 

1.26*00 

4 .16-35 

1.0136*03 

288.1 

S. 9E*  3 ! 

5.46-05 

1.0 

0.8206*0? 

259.1 

i.?f*o: 

4.1E-05 

8.986E*S? 

281.6 

4. 2E*0  1 

5.4E-05 

2.0 

2.2256*0? 

255.9 

9.4E-C1 

4. 16-05 

7.95t.6*C2 

275.  1 

2. 96*0  ; 

5.46-05 

3.0 

6.2906*0? 

252.2 

6. 8E-C1 

4.36-35 

7 . 0 1 26  ♦ 0 2 

26*. 7 

1. 6£*0  3 

5.02-05 

4.0 

5.9326*02 

242.2 

4.16-01 

4.56-05 

6 . 166E  * 02 

262.2 

1. 16*03 

4.56-05 

5.0 

5.1506*02 

240.9 

2. 1C- 01 

4.26-05 

5.4056*0? 

255.7 

6. 46-0 1 

4. 52 -y  5 

6.0 

4.4626*02 

234.1 

9.86-02 

4.36-35 

4.2226*0? 

24  j.  2 

?. 86-31 

4.56-05 

r.o 

3.0536*0? 

222.3 

5.46-02 

2.16  35 

4.1116*0? 

242.7 

2.  16-01 

4.96-35 

0.0 

3.5006*02 

220.6 

1.16-02 

9.  JE-05 

3.5  65  F *1)’ 

236.2 

1.2E-01 

5 . 2 C- 3 5 

5.  S 

2.029* *0? 

212.2 

0.46-03 

1.56-04 

3.05CE*62 

229.7 

4. 66-32 

7.16-05 

10.0 

2.4156*0? 

212.2 

5.56-03 

2.4E-04 

2.6506*0? 

223.2 

1.8E-3? 

9.02-J6 

11.0 

2. 0626*02 

21 2.2 

3.86-0? 

3.26-04 

2 .22C6*0? 

216.8 

8. 26-3  i 

1.36-34 

12.0 

1.2666*02 

212.2 

2.66-0? 

4.36-04 

1.940E*0? 

216.6 

3. 7E-J  3 

1.56-04 

13.0 

1.5106*02 

212.2 

1.86-0? 

4.26*04 

1.6566*0? 

216 .6 

1. 86-J  J 

1.76-0* 

14.0 

1.2916*82 

212.2 

1. 06-0  3 

4.3E-U4 

1 .4176*0? 

216.6 

8.46-04 

1.92-34 

15.0 

1.1036*02 

212.2 

2.66-04 

5.56-04 

1.2116*0? 

216.6 

7. 26-04 

2.16-04 

16.0 

9 4316*01 

216.6 

6.4E-C4 

6.26-04 

1.0S5F*0? 

216.6 

6. IE-04 

2.46-04 

ir.o 

6.0506*01 

216.0 

5.6E-0L 

6.26-04 

B.65CE*01 

216.6 

5. 2E-C4 

2.86-04 

10.0 

6.0026*01 

215.4 

5.06-04 

6.26-34 

2.5656*01 

21  6.6 

*. 46-0* 

3.2E-:4 

15.0 

5.8256*01 

214.8 

4.9E-04 

6,36-34 

6«467E*01 

216.6 

4.46-34 

3.56-34 

20.0 

5.0146*01 

214.1 

4 . 5E-04 

5.56-0. 

5.5296*01 

216.6 

l»,  46-04 

3.52-04 

21.0 

4. 2226*C1 

213.6 

5. 1E-C4 

5.16-04 

4.2296*01 

217.6 

4. 8E-04 

3.56-04 

22.0 

3.6426*01 

213.0 

5. 16-04 

4,26-34 

4,042E*01 

218.6 

5. 2E-J4 

3.96-34 

23.  n 

i.iAar.m 

212.6 

5 . 46-04 

4,36-34 

0,41>7t*Cl 

219.6 

5. 76-04 

3 . 5 E - i 4 

24.0 

2.6496*01 

211.8 

6.0E-04 

3.56-04 

2.9226*01 

220. b 

6. 16-34 

3.56-0'. 

23.0 

2.2566*01 

211.2 

6.26-04 

3.26-04 

2.5496*01 

221.6 

6.66-04 

3.42-04 

30.  0 

1.0206*01 

216.0 

3 .66-04 

1.56-04 

1 ,197E*01 

226,5 

3. 06-04 

■Z.  36-04 

35.6 

4.2016*00 

22c. 2 

1.  IE-04 

9.26-35 

5. 7466*0; 

236.5 

1. E_-Qw 

1.16-34 

40.0 

2.2436*00 

?34  .2 

4.36-35 

4,16-35 

2.8716*03 

25  3,4 

6. 76-05 

4.9E-C5 

45.0 

1.1136*00 

242.0 

1.9E-05 

1.36-05 

1.4916*03 

264.2 

3. 26-05 

1 . 7 F - 'j  5 

50.0 

5.2196-01 

253.3 

6.3E-06 

4.36-06 

7 <.97bf-0  1 

27  0.6 

1. 2E-05 

4.56-06 

20.  f 

4. 0166-02 

245.2 

1 • 46- 0 2 

8.56-38 

5.5206  *02 

219.7 

1.5E-C7 

6.5E-08 

100.0 

3.0CQ6-B4 

210.0 

1.06-09 

4.36-11 

3.CO0E-C4 

2 1 C . 3 

1.06-39 

4.32-11 

99999. 

210  .0 

21C.0 

* 

.20 

.28600  .09530 

.25  .28000 

,05660 

.31  ,26200  .02063 

.34  . 

245C0  .01450* 

.49 

.10500  .01050 

.51  .12600 

.0 1800 

,63  .14600  .09914 

.69  . 

13433  .30914* 

.06 

.10800  .01020  1 

.96  .08910 

.01380 

1.54  .05790  « CO  924 

2.00  . 

C 351 3 • 0 J 348* 

2.50 
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-0.63-0.40-8. 36-C. 28-6. 16-0. 06  0.00  0.20  0 -20  0.41  ;.5t  9 . 6 i 3,  8 0 0.92  1.34  6575 
1.19  1.19  1.0  0.49  1.0?  1.19  J . ? 9 1,30  1.29  1.30  1.19  1.33  1,4?  1.4J  1.70  6653 

1.62  1.54  1.41  1.53  1.86  1.96  1.97  2.02  2. Cl  1.94  1.94  1.03  2.03  2.21  2.42  6 ?25  J* 

2.30  2.16  2.02  2.8?  2.0?  2.13  1.90  1.71  2.01  1.56  1.50  1.51  1.30  1.63  1.54  6803  - 

1.67  1.78  2.22  2.39  2.38  2.30  1.93  239  2.49  2.52  2.57  2.21  2.10  2.40  2.41  6i75 

2.45  2.51  2.23  2.41  2.30  2,61  2.7?  2.  52  2,63  2.56  ?.S1  2.73  2.62  2.62  2.0ti  6950 

2.74  2.79  2.74  2.73  2.68  2.01  2,7?  Z.re  2,84  2.32  2.38  2.85  2„#S  3.02  3.08  7025 

3,26  3.03  3.1*  3.21  3.93  3.11  3.15  3. 1C  3.31  3.22  3.30  3.06  3.34  3*43  3.37  7101 

3.32  3.08  3.09  3.09  3.01  3.07  3.07  3. < 1 3.21  3. *1  3.67  3.5*  3.79  3.7u  3.49  7175 

3.39  3. 11  3.1?  3.01  3.10  3.01  3.18  3.42  3.43  3 5 3.40  3.39  3.39  3.51  3.54  7250 

3.42  3.50  3.67  3,59  3.63  3.66  3.40  3.39  3.29  >.31  3.41  3.2 3 3.32  J.12  2.9 1 7325 

2.91  2.75  2.70  2.7?  2.62  2.58  2-32  2.22  2.00  1.97  1.60  l.b?  1.64  i.53  1.56  7403 

1.51  1.5?  1.40  1.4?  1.42  1.40  1.41  1.43  1^56  1.52  1,51  1.52  1.39  1.39  1.33  7475 

1.89  1.16  1.21  1.20  1.22  1.20  1.10  1.20  1.19  1,17  1.10  1.10  1.09  1.19  l.U  7550 

1.04  0.98  0,99  i.§5  0.93  8.90  Q.90  0.86  3.71  0.79  0.70  0.7l  0.67  0.62  0,53  7625 

0.42  0.31  0.20  8.01 -0.00-0 .17-0. 26-0. 35-0 .44-0.53-0 .63-0.7 J-0. 03-0.93-1.04  7700 

-1.16-1. 24-1. 34-1.44-1.54-1.64-1. 74-1. 84-1 . 94-2. 0 4-2 . 14 -2 . 24 -2 , 34-2 . 44-2. 5 4 7775 
-2. 64 -2.74-2. 94-2. 94 -3 .04-3,14-3. 24-3. 34-3, 44-3. 54-j. 64-3, 04-3. 94  -4.14  7150 

-5, 80-5. 00-5. 80 -5. 80 -5. 09 -5. 00-5. 00-5. 00-5. 90-5. 00-5. OC-5.0j-5, 00-5.01-5. 00  7925 
-5.U0-5. 00-5. 00-5. 00 -5. 89-5. 00-5. 00-5. 00-5. 00-5. 00 -5. 00 -5. 60 -5. 00-5. 8 0-5. 00  8009 
-5.80-5. 00-5, 00 -5.0 0-5.  C 1-5. 88 -4.  00- 5. 00-5. 00-5. 08-5. 00-5.  If  3-5. 0t-5. 00 -5. 09  8,75 
-5.80 "5. 80-5. 90 -5.00 -9.00 -5 .00-5,  00-5,  JO-5, 00-5. 60 -5. 00 -5. 03  -5, 00 -5. 03- 5, JO  0150 
-5. 00 -5. 00-5. 00-5. 03 -5.00-5. 00 -5. 00- 5. 98 -5. 00-5. 90 -5. §0-5. 00*5. 00-5. 00-5. B*  0225 
-4. 15-4, 06-3. 97-3. 04-!- 71-3, 70 -3.6i-3. 52-3. 4 3-3. 34-3 .25- 3 , 16- J. 0 7-2 . 98-2 . 19  830. 

08 


SPECTRAL  DATA: 


Table  A2<  Listing  of  D«t*  for  LOWTRAN  3 (Cont) 


-*.00-2. 7 1-2. 6*- *.53-2. 66-* 
-i.66-1. 37-1.  28-1. 11-1.  10-1 

-i.u-mvm;  i. as  o.io  a 
a.M  i.5s  a.ri  a.n  a. as  a 

1.63  1.9;  1.9*  1.99  1.9*  1 
l.sa  l.ra  i.*3  1.90  i,99  i 
1.**  1.61  1.19  1.3?  1.9?  1 
f.il  *,*0  ?.?1  2.13  ?. 09  1 
1.99  1.27  1.20  1.18  1.11  0 

-2.9l-2.99-3.U-3.3r-3.96~3 
-9.ca-9.aa-9.ao-3.o0~9.ao~9 
-2.09-2.29-2. 29-2.6 3-2.62-? 
-*.*2-1.96-1. 22-1.20-1.63-1 
-0.93-0. *2-0. 91-0.03 -0. 21-0 

-O.Ol-0.0O-a.26-O.lfc-O.fl?  0 
0.90  0.83  0 . 66  a. 39  0.30  0 
0.32  0.32  0.30  0.66  0.62  0 
1.60  1.32  1.3?  1.37  1.6?  1 
1.6*  1.60  1.3*  1.36  1.63  1 
1.6*  i.6J  1.2|  1.7?  1.70  1 
1.26  1.20  1.13  1.16  1.19  1 
2.0*  2.19  2. *6  2.*?  2.02  1 
8.90  0.65  0.63  0.29  0.16  * 
-O.CS-O. 16-0. 21-0. 00-0. 06  0 
0.61  0.33  0.3?  0.91  0.90  0 
0.39  0.6?  0.39  0.3?  0.30  0 
-0.03  0.21  0.36  0.39  0.67  0 
0.99  0.93  0.90  0.32  9.2*  0 
-1.63-1. 69-1. 70- 1.91 -*.01-1 
-2. 37-2, 69-2. 66-2. !6-*.32-2 
-2. 60-2. 36-2. 60-2. 6»-2,89-? 
-3, 00-5. 00-5. 00-9. 00 -9. 03 -9 
-9. 00-9. 00-9.09-9. 03-9. 00-9 
-9. 18 -9,95-9.05=?. SC -5 .6 5 -5 
-9.00-3. 00-9. 00 -5. *0-3. 6«-5 
-3.21-3.96-3.66-3.21  -3.06-2 
-1.63-1. Of- 1.71-1. 59-1. 69-1 
-6. 91-1. #1-0. 99-6.07 -6.92-0 
-I. 30 -0.2f -0.30-1. *8-0. 01-1 
-6. 10-6. 06-0.05-0.06 -6.16-1 
-6.96-0.26  0.11  0.51  0. 01  0 

a. 09  a. 06  e.ao  i.n  0.21  a 

0.6?  0.3]  0.63  0.36  0.2?  I 
-0.00-1. 01-1. 10-1. 13-1. 29-1 
-1, 62-1. 60-1. *6-1. 38 -1.31-1 
-1.69-1. 60-1.36-1. 6? -*.66-1 
-2. 00-2. 00-2. ?2-2. 38-2. 35-1 
-1. 96-1.62-1. *9-1. 39-1. 20-1 
-?.0a-r.l*-l. 99-2, 21-2. 06-2 
-*.6*-*. 26-1  99-2,61-2. 16-2 
-1.71 -1.70-1. 72-1. *8 -1.70-1 
-0.90-1,19-1.08-0.79-0.60-0 
-9. *1-0.61-0. 29-0. 29 -0.61-0 
-0.91  0.10  0.20  0.11  a.  -0 
0.22  0.11  0.11  0.21  0.89  0 
-0.00-1. 01-0.06-0. 56-0. 19  C 
a. 61  0.69  0.60  0.66  6.6*  0 
-2.06-2.06. 2.16- i .96-2.00-2 
-*.36-2.36-1.91-1.02-1.69-1 
-1.69-1. 56-1, *6-1. 21 -1.36-1 


.55-2.26-2.10-2.09-2.00-1 
.81-0.92-0.03-11 ,26-0.65-0 
.17  0.??  0.39  8.35  8.61  C 
.11  1.06  1.87  1.82  1.81  1 
.51  1.61  1.50  1.7*  1.*?  1 
.55  1.60  1.3?  1.39  1.53  1 
.70  1.90  2.81  1.92  1.91  2 
.91  1.92  1.97  l.tO  1.91  1 
.99  3.06  0.71  8.68  8.66  8 
.11-1.26-1.61-1.57-1 .73-1 
,75-3. 96-6. 13-6.31-6.69-6 
.80-5. 00-5. <0-6. 68-6. 26-3 
.29-2. 20-2. 1 7-?.?3-2. ?7-2 
.60-1,59-1.63-1.21-1.15-1 
.62-C. 59-a.58-0.6S-0.50-: 
.18  0.3?  0.62  0.37  0.23  0 
.35  :.23  8.26  0.19  fl.CO  0 
.90  1.06  1.02  1.13  1.08  1 
.50  1.62  1.30  1.36  1.30  1 
.66  1.5*  1.69  1,69  1.52  1 
.78  1.67  1.6?  1.6*  1.70  1 
.29  1.56  1.72  1.06  1.70  1 
.90  1.98  1.03  1.01  1.72  1 
.05  G.O?  0.03  3.0J  0.01-0 
.10  0..S  5.11  0.3?  0.6?  : 
.66  C . 31  0.10  C.CO  0.20  c 
. 16-C.  e l-o. 2 3- 0.6 1-0. 5 2-C 
.66  C.60  0.51  3.59  0.53  ’ 
.11-0.08-0.  1 0-1. 16-0. 63-0 
.97-1.97-1,97-1 .97-2.26-2 
.19-2.10-?.  25-2. 16-2. 36-? 
.60-2.09-3.20-3. 51-3. 76-3 
.0C-S. 00-5.  0 0-5. 00 -5. 00 -5 
.00-5.00-5.00-5.00-5.00-5 
.06 -5.00-5, 0C-5, 00-5. 00-5 
.80-5.OC-5.OC-5.OC-5.CO-5 
.98 -2. 76-2.60-2.66-2.32-? 
.66-1.66-1.69-1 .69-1,25-1 
.79-0.62-0.56-8.30-0.62-0 
.20  :.#6  0.10  0.06  0.16-1 
.06-0.06-0,21-0 .39-0.6X-: 
.79  0.6?  8.26-0.31-0.67-0 
.13  9.32  0.35  0.91  t.60  0 
.13-0.1 1-a. I 1-0. 31 -C. 61 -J 
.65-1 .69-1.67-1.67-1 .51-1 
.30-1. 30-1.20-1 ,39-1. 33-t 
.61-1,62-1.60-1.61-1.31-1 
.90-1.92-1.78-1 .57-1.69-1 
.60-1.50-1,66-1.53-1.60-1 
.06-2,39-2. 76-3.09-3.66-3 
.31-2.15-2.01-1.99-2.16-2 
,52-1. 30-1. 29-1.22-0. 9i-3 
.*6-0.73-0.95-3 .89-0.61-; 
.76-0.19-0.10  3.  0.19-0 

.37-0.10  9.0?  0.16  0.20  3 
.21  C.20  0.37  (.26  0.0? 
.19  3.23  6.21  1.29  0.20  3 
.27  C. 06-0. 33-8. 81-1. 17-1 
,»8  2.00-2.23- 2 . 31-2.31-2 
.56-1.06-1.91-1 .75-1.03-1 
.36-1.66-1.60-1 .60-1.60-1 


.91-1.02-1. 73-1. **-i. §5  0375 
.56-0. 67-;. 39-0. 29-0.20  6650 
.65  0.6?  0.60  i. 61  6.6 6 0525 
.12  l.?3  1.26  1.29  1.36  060G 
•9?  1.96  1.09  1.01  1.65  0675 
.02  2.23  2.61  2.51  2.23  6750 
.12  2.13  2.61  2.18  1.99  0025 
.91  1.92  1.93  1.76  l.il  0903 
.31  0.13  9. 03-3. 07-0. 21  0975 
.91-2.09-2.27-2.65-2.il  9;50 
.66 -6. 83-6. 99-5. 16-5. 20  9125 
.09-3.57-3.32-3.11-2.91  9875 
.32-2.12-2. 00-2. 07-2. )7  9950 
.39-1.13-1. 29-1. 1 9- J .3010325 
.39-0.22-0. 16-0.06-0.3110100 
.12  0.15  0.20  ;.6]  C. 5913175 
.10  8.18  0.27  0.39  ..5310280 
.80  1.1b  1.16  1.30  1.6110325 
.69  1.63  1.6?  1.62  1.7310603 
.SO  l«Dc  x » 6 c i .61  1.5110675 
.6?  1.56  1.69  1.62  1 . 30 1 u 55  3 
.02  1.03  1.0?  1.09  1.3310625 
.69  1.59  1.50  1.36  1.20167C3 
.OB -0.10- 0.20-0. 11 -0.061 0 775 
.66  0.30  8.20  0.6?  C, 6110050 
.21  0.36  ..36  0.29  0.3510925 
.60-0. 5 »-0. 6 1-0. 6 S-l.J 311600 
.69  0.57  0.66  0.52  0.6211C75 
.62-0.00-1.09-1.16-1.3111150 
. 20 -2. C 1-1. 99-2.0; -2. 0611 225 
.66-2.6ii-2.69-2, 66-2.  6317.  3CC 
, 97-6. 2C -6. 6 3-6. 66-6. 5 911375 
.80-5.00-5.00-5.00-5.0011650 
. Bfl -5. C 3-5=00-5, 00-5. 051152' 
.CC-5.0; -5. CD -5, 0 3-5.  00 11600 
.0 0-5.00-5.00-5.02-5.0011675 
.17-2.0 J-l.  07-1 .79-1.7  611750 
.26-1. 3 0-0. 90 -1.06-3,9  HI  025 
.68-0.36-0.27-0.17-;. 2011900 
.12-0.02-8. 62-0.13-0,1111975 
.60-0.31-0.62-0 .55 -C. 5 712350 
.00-8. 00-0. 56-0. 39-0.1 01*125 
.51  8.51  0.60  0.60  C, 6312200 
.61-0. 39-0. 5 3-S. 69- C. 1612275 
.66-1  .60-1. 69-1. 03*1. 5112 350 
.60-1.35-1. 37-1.39-1.6112625 
.15-1.1 J-l. 20-1. 61-1. 001 ?*S0 
. 78- x. 73 -1.66-1. 06- 1.5 0125 75 
.65-1 .58-1. 58-1. 69-1. 791**50 
. 79-6. l%-6. 69-6. 86-5. 19 12 725 
,61-2,12-1.99-1.86-1.7913608 
. 90-1. C l-O. 76-0.98-0.9013675 
.61-0.69-0. 51-0.97-6. 8313 550 
.10  0 .20  0.20  0.02  .2013625 

0.09  8.09  9.99  0.0713708 
. J9-8. 29-9. *9-J. 69-0. 76iJ775 
.29  8.5?  9.56  3.51  .6011058 

.11-1.37-1. 52-1 . 56- 1,96 i 3925 
.53-2,31-2.31-?. 31-2.2816081 
.76-1 .56-1. 9»-i  .80 -l.it  16.5 75 
.36-1 .65-1.69-1.05-1.3916150 


S?ECTRAL  CAT*: 


T 


e A2. 


Listing  of  Data  for  IjOWTHAN  3 (Conti 


-1.23-1.16-1. 13-1. 34-1. 36-1 .23-1. ?3*1. 3 7-1 
-l.??-l,2«-l.'?5-1.69~2.0  7-2.42-2.56-2.  5S-? 
-I « 23-1. 10- l . 23-i . 1C -0.81 -fl .80-: . 60-0. 60-0 
-1. 50-1.00-2. 16-2. 3? -2, 6 3-3.91  -4,22-4.4  9-4 
-6. *5 -3.70-3. 20 -2.75-1. 90-1. 73-1.51-1. 26-1 
8.49  0.76  1.06  1.23  1.56  1.76  1.91  2.08  2 

3.66  3.66  3.74  3.66  3.88  3.66  1.73  3.58  S 

2.01  1.64  1.77  1.63  1.67  1.21  ’.92  C.5J  * 

-1.61-1. 18-1. 62-1. 61 -1.66-2. 10-2. ’9-2.51-2 
-5. 00 -5. 00-5.03 -5 .00-5.98-5.00-2, 68-2.5  7-2 
-1.09- 1.11 -1.13 -1.09-1. 01 -1.0 1-1 .11-1. 33-1 
-1. ’8-1. 59-1. 33-1. 13 -1.0 1-8. 96-". 9! -C. 90 -C 
-2. 11-2. 75-1. #9-3. 50-5. 03-2. 60-2.23-1. 86-1 
-1. 21-1.17-1, \2-l. 15-1.19-1. 2G-1. 17-1.32-0 
0.67  0.77  0.46  1.07  1.11  1.11  1.06  1.15  1 

0.5%  0.66  0.76  0.61  0.66  6.6?  :.6S  C.47  0 

-2.61 -5. 60-6. 6C-5, 63-3. l%-2.%7-2. CO-1.71-1 
-2.0% -2. 10-2. 23-?. 3?-?. %6-2.7i-?. 66-3. 39-2 
-1.7* -2. 01-2.%1 -2. 65-2.63-2.59 36-2. 27-2 
-l./Z-l. 81-1. 92-2.0 1-2. 27-2. 61 -1.2 5 -5. 3 1-5 
-5*66-5. 00-5, 90 -5.0* -5. #|-S.0C-5.0C-6. 00-5 
-5. 80-5. C6-5. 60-5. 03-5.00-5. 00-5. 00-5. 10-5 
-2.83-2,71-2.67-2.67-2.68-2.58-2. 13-2.01-1 
-0.63-0. 69-0. 67- 1.0 8-1.26 -1.5 3-1 .87-1.91-1 
-3.09-2.93-2.76-2. 39 -2. 01 -1.69-1. 36-0.99-3 
0.0?  0.0 1-0.06-0. 2 '-0.60 -0.51-;. 5 3-4. 5 7-C 
-0.91-8. 58-9. 51 -0.01  0.19  0.51  ..76  1,00  l 

1.91  2.02  2.10  2.3?  2.5|  2.61  2.69  2,01  2 

3.72  3.98  5.03  5.27  5.5?  5.61  5.71  5.73  5 

3. ?8  2.79  ?,30  1.85  1.35  0 .62 - ? . 2 5- 1. 69 -2 

-1. 1?- 1.1 5- 1.25- 1.26 -1.20-1. 17-1 . 28-1. 32-1 
-1.34-1.06-0.81-0.61-0.59-0.55-3.57-0.59-3 
-1. 43 -2. 01-?. 60 “2. 89 -2. 67-2. 74 -2. M-2. 42-2 
—2-63  — 2.  ,6  — 2.53- 2 . v * - 7 . ? — — * • o < — . . a . — c . 6a  — 2 

-2. ?fe -2. 2 3- 2. 20 -2, 03 -2. C2-1. 96 -1.85-1. 55-1 
-1. 65 -1.69-1. 75 -1.79-1. 87  1.78-1. 63-1. 5C-1 
-0.38-0.19-0.B9-0.05  0.02  0.10  0.1b  0.18  0 
0.12  0.97-0.01-8,07-0.09  0.32  C.72  0.91  1 
-0.17-0.08-0.88  0.09  0.13  i.18  ?.25  0,27  3 
0.89  0.02-0, 85 -8. 13-6. 32-C.51-C.72-C. 98-1 
--2. 62-2.  8 7-3. 03-3. 21 -5. 0 8 -5. 80 -5.  0 0-5.80-5 
-5.88-5. 0 1-3, 38-3.01 -2. 63-2 • 12-2, 0 9- 1. 98-1 
-1. 69-1, 51-1, 38 -1.55 -1.65 -1,90 -2. 99-2, 55-2 
-3.37-3.38-3.16-3.01-2.76-2.51-2. ’0-1.80-1 
0.20  8.36  8.51  0.61  0.67  0.83  1.30  1.22  1 

2.57  2.61  2.76  2.9?  3.81  3.85  J.O?  2-98  2 

1.83  1.59  1.59  1.50  1.67  1.96  2.22  2.5C  2 

3.26  3.19  2.98  2.59  2.15  1.7C  1.22  C.55-0 

-2. 6 W. 68-2. 63-2. 60 -.'.57-?. 5 3-2. 5 7-2.65-2 
-5.88-5.88-5.88-5.88-5.08-5.80-5. 00- 5.80 -i 
-5. 88-5. 88-3.73-3. 62-3. 59-3. 53-3, 56-3,57-1 
-3.2  7-3. $6-3.68 -3.96 -5. 08 -5. 80 -5.00-5.  0 0 -5 
-5 ,88-5.88-5. 88-5. 98 -5. 08-5. 8C -5.  0 8-5.0  0-5 
-5.8  7-3. 89-3. 76-3.67-3.56-3.52-3.  3 5-3.26-3 
-3.38-3, 38-3.37-3.24-3.15-3.88-3.88-2.93-2 
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Appendix  B 

Basic  Fla.  Chart  tor  LOWTRAN  2 


A general  flow  chart  for  l.OWTRAN  3 is  given  in  Figure  B1  which  shows 
the  overall  mode  of  operation  of  the  program.  More  detailed  flow  charts  are 
also  given  for  the  two  main  blocks  in  the  program,  that  is,  where  the  equivalent 
absorber  amounts  and  refraction  calculations  are  made  (Figure  B2)  and  for  the 
transmittance  calculation  (Figure  B3). 

The  notation  used  in  the  flow  charts  is  as  follows: 

(1)  If  a condition  stated  within  a given  block  is  fulfilled,  then  the  direction 
of  flow  is  sideways  as  indicated  by  the  direction  in  which  the  block  points  (for 
example,  — • for  the  following  block  f > ), 

(2)  If  the  condition  stated  within  a block  is  not  fulfilled,  the  flow  is  down- 
wards. 

The  numbers  appearing  on  the  flow  charts  correspond  to  the  statement  numbers 
given  in  the  main  program  (see  Table  Al). 
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Figure  Bl.  General  Flow  Chart  for  LOWTRAN  3 
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Figure  B2.  Flow  Chart  for  Calculation  of  Equivalent  Path  Quantities 
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Figure  B3.  Flow  Chart  for  Transmittance  Calculations 
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Appendix  C 

Iterative  Refraction  Scheme  Uiad  in  Subroutine  ANGL 


I.  INTRODUCTION 

An  iterative  scheme  is  presented  here  for  determining  the  initial  zenith 
angle  0 ^ required  for  a path  trajectory  from  altitude  HI  to  H2  for  a given  earth 
center  angle  3Q  (see  Figure  Cl),  taking  into  account  refraction.  The  theoretical 
background  for  the  scheme  can  he  conveniently  divided  into  upward-  and  downward- 
looking path  trajectories. 


2.  IIPfARD  TRAJECTORIES  (0q  < 90°) 

If  we  let  the  geometrical  angle  9 in  Figure  Cl  be  our  initial  guess  for  & and 
use  this  In  the  LOWTRAN  3 program,  we  can  calculate  the  corresponding  3, 
taking  refraction  into  account  (fi  >.3Q  for  upward  looking  paths).  It  is  apparent 
that  3 is  a function  of  9 for  a given  HI,  H2  and  model  atmosphere.  That  is: 

3 * f<9)  . (Cl) 

We  can  then  differentiate  3 with  respect  to  6 and  write  the  differential  d3  = t'(9  )d0 
in  the  form 

dfl  = «!3/f»<0)  . <C2) 
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Figure  Cl.  Schematic  for  Upward  Looking  Paths 
Showing:  (a)  The  first  guess  at  the  trajectory  using  the 
geometrical  angle  8 from  Hi  to  H2,  and  fb)  the  final 

iteration  where  6-0  and  B - & 
o o 


Replacing  the  differentials  by  differences,  we  have 

8'  - 0 - <0o  -0)/f<8)  (C3) 

where  is  the  next  guess  at  0 . That  is: 

o 

d*  - 0 +<0o -0)/r>(0)  . (C4) 

Thun  we  can  substitute  6 ' in  Eq.  (Cl)  and  go  through  the  same  procedure  described 
above  to  ottain  successive  iterations,  until  0 > converges  finally  to  0o  when 

is  negligibly  small  (for  example,  < 10  ^ radians).  The  major  unknowns  in 
the  above  equation  are  1(6)  and  f’(8)  . From  Eqo.  <9>  and  (10)  of  Section  4.  2,  it 

was  seen  that 

m-l 

? = t <«.  -o.)  = f(8)  . (C5) 

i 1 1 
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Thus 


m- 1 

f’<0)  * Z 
i 


/ d6|  doA 

\ de  ' de  / 


(C6) 


■where  m is  the  number  of  levels  between  HI  and  H2.  Also  from  Eq.  (7) 

sin  8 = n (R  + HI)  sin  8 /n.(R  + z.)  . 

i o o '10  1 

Differentiating  the  above  equation  with  respect  to  8,  we  have 

de 

coa  ei~d r nomo  + H1>  cos  6 'ni(Ro  + 2i)  • 

Therefore 


(C7) 


d0  . tan  8 . 

i _ i. 

d<3  ' tan  8 


(C8) 


assuming  that  m and  z i are  independent  of 


Similarly,  it  can  be  shown  that 


doj  tan  o j 

d 0 tan  fl 


' (C9) 


so  that 


f'(fl) 


1 

tan  9 


m -1 

£ (tan  6 ^ * tan  o^)  . 
i 


(CIO) 


Equations  (C4>  and  (CIO)  form  the  basis  for  the  iteration  scheme  provided  in 
the  subroutine  ANGL  which  contains  all  the  angular  calculations  given  in  the  main 
LOWTRAN  2 program  for  both  upward-  and  downward -looking  paths.  For  an 
initial  guess  0,  the  quantities  0 and  f'(®)  are  calculated  by  summing  3 j and 
tan  dj  - tan  Oj  (defined  in  Figures  2 and  Cl)  for  each  layer  along  the  total  extent 
of  the  path  trajectory.  Using  Eq.  (C4),  a second  guess  8 ' is  obtained  and  is  used 

to  obtain  successive  9 's  until  the  final  value  of  0 is  sufficiently  close  to  (3  (for 

-7  ° 

example,  to  within  10  radians).  When  the  latter  condition  is  satisfied,  the  final 

iteration  gives  the  required  initial  zenith  angle  gQ. 

t This  assumption  holds  for  upward  looking  trajectories  but  not  for  some 
downward  looking  cases  where  the  trajectory  passes  through  the  minimum  height 
HM1N  (see  Figure  3).  The  latter  case  will  be  discussed  In  the  following  section. 
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3,  DO»N*ARD  TRAJECTORIES  (0c  > W) 


The  iteration  scheme  for  determining  for  downward  looking  trajectories  is 
Identical  to  that  described  above,  provided  that  the  path  does  not  pass  through  the 
minimum  height  HMIN  [that  is,  see  Figure  3(d)).  If  the  path  does  pass  through 
HMIN,  however,  it  is  apparent  that  the  refractive  index  between  HMIN  and  the 
level  above  (say,  n ) and  HMIN  are  both  dependent  on  successive  guesses 
of  $. 

Consider  the  case  shown  in  Figure  C2  where  >?2  < Hi.  For  a given  initial 
guess  S,  let  us  divide  the  contribution  to  the  total  angle  0 into  three  parts: 


Figure  C2.  Schematic  for  Downward  Looking  Paths  Showing  the 
First  Guess  at  the  Trajectory  Using  the  Geometrical  Angle  9. 
HMIN  is  the  minimum  height  for  the  first  iteration,  which  will  be 

lower  than  the  final  value  when  9 - 0 „ and  0=0 

o o 
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Bj,  3 2 and  fiy  Here  £5^  is  the  angle  between  HI  and  H2  (on  the  short  side  of  the 
path),  is  the  angle  between  H^  and  the  level  * above  HMIN,  and  £3  is  the 


angle  from  zf  to  HMIN.  That  is.  0 = + 202  + 203  = f(0). 


f'(fl) 


40,  2d0„  2d0- 

de»  + ~59~+  ~dS~- 


Then 


(Cli) 


Following  the  same  procedure  ou.Tined  in  the  previous  section,  and  defining 

sin  6 . s n <R  + Ht)  sin  8 in.  (H  +2.)  (C12) 

i o o 1 ' 1 o i 

where  n now  refers  to  the  mean  refraction  from  HI  to  the  level  below  HI.  we  can 
show  that 


k-1 

1 £ (tan  0 . - tan  a ■)  - F. 

d®  ' tan  6 j 1 11 


d(3. 


/-I 


d ^Zl^re  l (tan  0 . - tan  aj  - 


d(3  d8 

-J.  - - — HL  - F 

53^  ' " d»  *2 


(where  6 90  . p ).  Now 

m j 


sin  6 n (R  + H. ) sin  0 /n  (R  + z J 
moo  1 'mol 


'"•.“Tv  R io  the  refi  active  iiiuc-x  between  KjVIIN  and  z. 

m 1 

with  respect  to  0 and  dividing  by  cos  0 we  have 

m' 


(C13) 

(C14) 


(CIS) 


(C 16) 


r\  t . ..  ...  »■»«  jo  ■ «»\ 

L/UieX'eiK&MVUl^  £«Q.  \c,iu/ 


d8  tan  6 tan  6 / dn  \ 

in  _ rn  _ m I m ] 

d6  ~ tan  9 ~ n \ d®  / ' 

m 


(C17) 


We  will  now  digress  for  a moment  to  show  some  relationships  between  HMIN, 
nm  anil  8,  which  will  assist  us  in  defining  the  quantity  dn^/dO  . If  we  let 
X =•  R0  + HMIN  and  write 


dn 


m 


dX 


dX 

d8 


(C18) 
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j 


then  clearly  dnm/dX  is  the  mean  refractive  index  gradient  above  HM1N,  which  we 
can  specify  for  a given  model  atmosphere.  The  other  parameter  dX/dfl  defines 
the  variation  of  HMIN’  with  8 . 

Applying  Eq.  (C12)  and  replacing  n ar  and  ® . by  n HMIN  and  v,  2 
respectively,  we  have 

n X = n (R  + HI)  sin  8 . <C19) 

m o o 

Differentiating  Eq.  (C19)  with  respect  to  ® and  using  Eq,  (18),  it  will  be  seen 
that 


(C20) 


Also,  if  we  assume  that  n - 1 varies  exponentially  with  altitude  where  n is  the 
mean  refractive  index,  then  the  refractive  index  gradient  is  given  by 


dn  n- 1 

ax  ' it 

where  H is  a scale  height  parameter  which  can  be  defined  as  follows: 
H = (*i  - *i_1>/loge  j (n^  - X)/(n.-l)  } . 


Therefore 


dn 
— to 
IdT 


■-(^brih )■»«.( 


(C21) 


Using  Eqs.  (C18)  through  (C20)  and  substituting  for  dn^/d®  in  Eq.  (C17),  we  have 


d®  tan  0 


— "Jjffl  [,  . _L  an  h\_  _1_  <in\ 

tan  6 1 n dX/  \X  n dX  ) 

m / ' m ' . 


(C  22) 


Thun  wc  can  now  write  Eq.  (C15)  in  terms  of  known  quantities: 
tan  8 


F3 


tan  0 


/ / dn  x*1 

_l-(1+nn/XlDf  ) J* 


<C23) 
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The  final  expression  for  the  iteration  scheme  for  the  general  case  where  H2  < HI 
and  6 > 90°  is 

O'  - 6 i (p  - ji)/(Fl  + 2F2  + 2F3>  (f'24) 

where  Fj.  and  are  defined  in  Eqs.  <Cl3),  (04)  and  \C23/.  For  t e case 
where  0 > 90°  and  H2  > Hi,  the  above  expression  would  become 

0 1 = e + (do  - i3)/(2Fl  + 2F3  + F) 

where  F refers  to  the  upward  looking  contribution  from  Hi  to  H 2 and  is  giver,  by 


F = 


tan  6 


m = 1 

£ (tan  0 . - tan  a 
i -r 


In  the  above  expression,  0 is  the  angle  of  refraction  ax  III. 


0 = sin" 1 (n  sin  9/n1 ) 

r*  n ' n 


That  is 


W'here  n'  and  n are  the  mean  retractive  indices  of  the  layers  above  and  below  HI 
o o 

respectively,  and  9 is  the  current  initial  zenith  angle  guess. 
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Appendix  D 

Symbols  art!  Definitions 


A HAZE 

AHZZ 

AJ 

A LAM 

ALP 

ANGLE 

AO 

AVW 

BET 

BETA 

B.I 

CA 

CO 

cvv 

Cl 

C2 

C3 

C4 

C5 

ce 


Absorption  at  frequency  \ , aiso  average  transmittance 
Aerosol  number  density  for  MOPEI.  - 7 
Aerosol  number  density  for  MODEL  - 7 
Equivalent  absorber  amount  per  km  at  level  .1 

Wavelength  ( ^ m) 

Angle  of  arrival  at  adjacent  level 

Input  zenith  angle  (degrees'  (compare  with  (1  in  the  text] 

Constant  A defined  in  Kq.  (10':  that  is.  (11  ■*  H i hi  sin  <! 

1 o o o 

Average  wavelength  used  in  refractive  index  expression 

Angle  subtended  at  the  earth's  center  as  path  traverses  adjacent  levels 
jcf/L  in  Eq.  ini] 

Total  angle  subtended  by  path  at  earth's  center  (compare  £ in  Eq.H'l] 
Equivalent  absorber  amount  per  km  at  level  J + 1 
Conversion  factor  from  degrees  to  radians 

Wavelength  dependent  coefficient  used  m refractive  index  expression 

Wavelength  dependent  coefficient  used  in  refractive  index  expression 

Log  absorption  coefficient  for  water  vapor 

Log  absorption  coefficient  for  uniformly  mixed  gases 

Log  absorption  coefficient  for  ozone 

Absorption  coefficient  for  nitrogen  ('4  pm' 

Absorption  coeff. merit  lor  water  vapor  continuum  (-  lowin' 

Extinction  coeffu  ient  for  molecular  scattering 


Preceding  page  blank 


10b 


C7 

CIA 

C8 

D 

DP 
DS 
DV 
DZ 
E(K) 
EH(1.  I) 
EH( 2,  It 
EH(3,  n 
EH(4, 1) 
EH(5,I> 
EH(6. 1) 

EH(7. 1) 
EH(8, 1) 

EH(9, 1) 

EV 

F 

FO 

FW 

HI 
II 2 

HAZE 

HM 

HMIN 

HZ  1 

HZ2 

I 

IATM 

IDV 

IFIND 

IHAZE 

IM 

IP 


Extinction  coefficient  for  aerosol  models 
Aerosol  absorption  coefficient 

Absorption  coefficient  for  ozone  (UV  and  visible  regions) 

Water  vapor  amount  (pr.  cm/km)  at  level  J 
Dew  point  temperature  ("C) 

Path  length  from  level  I to  Level  1 + 1 

Wavenumber  increment  at  which  transmittance  is  calculated 
Height  increment  from  level  I to  level  I + 1 

Equivalent  absorber  amounts  per  km  at  height  HI  [see  in  Eq,  (2)] 
Equivalent  absorber  amount  per  km  for  HgO  at  level  Z(I) 

Equivalent  absorber  amount  per  km  for  COg  + NgO  etc  at  level  Z(I) 
Equivalent  absorber  amount  per  km  for  Og  at  level  Z(I) 

Equivalent  absorber  amount  per  km  for  Ng  at  level  Z(I) 

Equivalent  absorber  amount  per  km  for  HgO  continuum  at  level  Z(I) 

Equivalent  absorber  amount  per  km  for  molecular  scattering  at 
level  Z([) 

Equivalent  absorber  amount  per  km  for  aerosol  extinction  at  level  Z(I) 

Equivalent  absorber  amount  per  km  for  ozone  (UV  and  visible)  at 
level  Z(I) 

Mean  refractive  index  of  layer  above  level  Z(I) 

Integrated  absorber  amount  from  level  I to  level  I + 1 [cf  W.  defined  in  Eq.  (lb)] 

^ 3 

Function  for  determining  saturation  vapor  density  of  water  (gin  m ) 

Transmission  function  logarithmic  absorber  amount  scale  for 

Transmission  function  logarithmic  absorber  amount  scale  for  HgO 
and  the  uniformly  mixed  gases 

Initial  altitude  (km) 

Final  altitude  (km) 

-3 

Aerosol  number  density  (no.  cm  ) 

Estimated  tangent  height  (km) 

Minimum  altitude  of  path  trajectory  (km) 

_ 3 

Aerosol  number  den?  ity  (no.  cm  ) for  23  km  visual  range 

.3 

Aerosol  number  density  (no.  cm  ) for  5 km  visual  range 

Running  integer  used  as  altitude  (level)  indicator  and  frequency 
indicator 

Number  of  levels  ir.  model  atmosphere 
Frequency  increment  (cm  *) 

Indicator  for  using  subroutine  ANGL 
Aerosol  model  indicator 

Parameter  used  when  reading  in  a new  atmospheric  model 
(see  Section  5.  2.  1) 

Indicator  for  using  subroutine  POINT  to  calculate  refractive  index 
only  (IP  = 0)  or  equivalent  absorber  amounts  also  (IP  4 0). 
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I 

1 


ITYPE 
IV 
IV I 
IV2 
DfY 

J 

JMIN 

JP 

J1 

J2 

K 

K1 

K2 

K4 

L 

LEN 

LI 


L2 


M 

ML 

MODEL 

Ml 

M2 

M3 

N 

NH 

NL 

NP 

NP1 

NP2 

P(M,I) 


Indicator  for  type  of  atmospheric  path  (see  Section  5.1) 

Frequency  at  which  transmittance  is  calculated 

Starting  frequency  (equivalent  to  VI  — see  Section  5.1) 

Last  frequency  (equivalent  to  V2— see  Section  5.1) 

Parameter  for  terminating  program  and  cycling  indicator 
(see  Section  5.1) 

Running  integer  for  altitude  identification 
Altitude  indicator  for  minimum  heigh  of  path 
Print  option  parameter 
Level  indicator  for  altitude  HI 
Level  indicator  for  altitude  H2 

Absorber  indicator,  K = 1.  2,  3,  etc.,  corresponds  to  HgO, 
uniformly  mixed  gases,  Og  etc,  respectively 

Integer  used  in  reading  two  model  atmospheres  on  one  card 

Integer  used  in  reading  two  model  atmospheres  on  one  card  and 
cycling  parameter  for  downward  looking  paths 

Frequency  indicator  for  nitrogen  continuum  transmittance  calculation 

Frequency  indicator  for  ozone  transmittance  calculation 

Parameter  used  for  defining  longest  of  two  paths  (see  Section  5.1) 

Frequency  identifier  for  UV  and  visible  ozone  transmittance 
calculation 

Frequency  identifier  for  UV  and  visible  ozone  transmittance 
calculation 

Integer  used  to  identify  required  model  atmosphere 
Number  of  levels  in  radiosonde  data  input  (MODEL  = 7) 

Integer  used  to  identify  required  model  atmosphere  (ree  Section  5,1) 
Integer  for  selecting  HgO  altitude  profile  for  (M=M1) 

Integer  for  selecting  temperature  altitude  profile  for  (M=M2) 

Integer  /or  selecting  Og  altitude  profile  for  (M=M3) 

Indicator  for  level  below  given  input  altitude  used  in  rOINT  sub- 
routine; also  as  frequency  indicator  in  UV  and  visible  ozone 
transmittance  calculation 

Frequency  indicator  for  water  vapor  continuum  transmittance 
calculation 

Number  of  levels  in  model  atmosphere  data 

Indicator  for  determining  whether  HI  or  H2  coincide  with  levels 
in  the  model  atmospher 

Value  of  NP  for  altitude  HI 
Value  of  NP  foi  altitude  H2 

Pressure  (mb)  at  level  I for  model  atmosphere  M 
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I 

i 

i 


PHI 

PI 

PPW 

PS 

PSI 

PT 

RANGE 

RE 

REF 

RH 

RN 

RO 

RX 

R1 

SALP 

SPHI 

SR 

SUM 

SUMA 

TIM.  I) 

THET 

THETA 

TMP 

TR 

TS 

TT 

TX(K> 

TXI9) 

TX(10) 

TX1 

TX2 

TX3 

V 

VH(K) 

VIS 

VX 

VI 
V2 


Angle  of  arrival  at  H2 
3.141592654  that  is  <TH 

Partial  pressure  of  water  vapor  (in  atmospheres) 

Total  pressure  in  atmospheres 

Angular  deviation  of  path  from  initial  direction  [cf  <l>  in  Eq.  (12)] 
Product  of  total  pressure  (atm)  and  the  square  root  of  273/T(M,  I) 
PaU.  length  (km) 

Earth  radius  (km)  (cf  R in  text) 

o 

Refractive  index  of  air  at  level  I 
Relative  humidity  (%) 

Ratio  of  refractive  indices  of  air  above  and  below  a given  level 
Earth  radius  (km)  read  in  as  input  (=  RE) 

Ratio  of  earth  center  distances  between  adjacent  levels 

The  product  of  the  sine  of  the  initial  zenith  angle  and  the  earth  center 
distance  to  starting  altitude 

Sine  of  angle  of  arrival  at  adjacent  level  (cf  sin  a) 

Sine  of  the  local  zenith  angle  at  a given  level  (cf  sin  0 ) 

Slant  range  (km) 

Sum  of  the  optical  thicknesses  of  absorbers  4 through  H 
Accumulated  integrated  absorption 
Temperature  (°K)  for  model  atmosphere  M at  level  i 
Zenith  angle  at  a given  level  (in  radians) 

Zenith  angle  at  a given  level  (in  degrees) 

Ambient  temperature  (°C) 

Transmittance  scales  for  transmission  functions 

Ratio  of  standard  temperature  (273. 0°K)  to  temperature  at  level  I 

Ratio  273. 15/ (TMP  + 273.15) 

Equivalent  absorber  amounts  per  km  at  a given  altitude  obtained 
from  POINT;  also  transmittance  values  at  a given  wavelength  for 
each  absorber  type  (K  = 1,8) 

Total  transmittance  at  frequency  V 
Absorption  due  to  aerosol  only  at  frequency  V 
Refractive  index  of  layer  above  initial  altitude  HI 
Refractive  index  of  layer  above  final  altitude  H2 
Refractive  index  of  layer  above  minimum  altitude  HMIN 
Running  frequency  (cm  *) 

Integral  of  the  equivalent  absorber  amounts  from  HI  to  level  I 
Visual  range  (km)  at  sea  level 

Wavelength  at  which  aerosol  coefficients  are  read  in  ( pm) 

Initial  frequency  for  transmittance  calculation,  cm 
Final  frequency  for  transmittance  calculation,  cm  * 
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Total  equivalent  absorber  amount  for  entire  path 

_ 3 

Water  vapor  density  for  atmospheric  model  M at  level  1 (gm  m 1 

_3 

Ozone  density  for  atmospheric  model  M at  level  I (gm  m 1 
Transmission  function  scaling  factor  for  H^O  at  given  wavelength 

Transmission  function  scaling  factor  for  COj  , etc.,  at  given 
wavelength 

Transmission  function  scaling  factor  for  O-  at  given  wavelength 

J -3 

Water  vapor  density  for  atmospheric  model  M at  level  1 + 1 (gm  m ) 

Input  height  to  POINT  subroutine 

Wavenumber  interpolation  parameter  in  UV  ozone  transmittance 
calculation 

Wavenumber  interpolation  parameter  in  HjO  continuum  calculation 
Wavenumber  interpolation  parameter 

Wavenumber  identification  parameter  for  UV  ozone  transmittance 
calculation 

Earth  center  distance  of  level  1 
Eartv  center  distance  of  level  1+1 
Input  zenith  angle  in  radians 

Refractive  index  of  layer  below  input  height  from  POINT  subroutine 
Refractive  index  of  layer  below  initial  altitude  If  - 
Refractive  index  of  layer  below  final  altitude  Hi- 
Aerosol  absorption  coefficient  at  frequency  V 
Altitude  at  level  f in  km 
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